1/2 


x  AD-A178  745 


MULTIDIMENSIONAL  FILTERING  INVESTIGATIONS^)  STEVENS 
INST  OF  TECH  HOBOKEN  N  J  S  BHSU  FEB  87  RADC-TR-8E-231 
F3B682-81-C-B185 


UNCLASSIFIED 


F/Q  12/1 


electe 

MAR  3  teW 


*E 


ROOM  AM  DEVELOPMENT  CENTER 
Air  Tore*  Systems  Command 
Orlffiss  Air  Fore*  Rasa,  NY  13441-9700 


1 V  |J 


UNCLASSIFIED 


AJ>-  A  / 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBS R(S) 
N/A 


(a.  NAME  OF  PERFORMING  ORGANIZATION 
Stevens  Institute  of  Technology 


AOORESS  (City,  State,  sad  ZIP  Cod*; 

Department  of  Electrical  Engineering 
&  Computer  Science 
Hoboken  NJ  07030 


Bb.  OFFICE  SYMBOL 
(If  spptksbh) 


Sa.  NAME  OF  FUNOING/  SPONSORING 
ORGANIZATION 

AFOSR 


Be  AOORESS  (City,  Stan,  and  ZiPCodt) 
Washington  DC  20332 


11.  TITLE  (Muds  security  Osaifkstion) 
MULTIDIMENSIONAL  FILTERING  INVESTIGATIONS 


lb  RESTRICTIVE  MARKINGS 
N/A 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 
RADC-TR-86-231 


7a.  NAME  OF  MONITORING  ORGANIZATION 
Rome  Air  Development  Center  (DCCD) 


7b.  AOORESS  (City,  State,  and  ZfPCode) 
Grlfflss  AFB  NY  13441-5700 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F30602-81-C-0185 


10.  SOURCE  Of  FUNOING  NUMBERS 


PROiECT 

NO. 


12.  PERSONAL  AUTHOR(S) 
Dr.  Sankar  Basu 


13a.  TYPE  OF  REPORT 
Final 


IE.  SUPPLEMENTARY  NOTATION 


14.  DATE  Of  REPORT  (Yssr,  Month,  Dsy)  IIS.  PAGE  COUNT 

February  1987  I  142 


17.  COSATl  COOES  |  H4.  RIBJECT  TERMS  (Continue  on  mono  if  naceoary  and  identify  by  block  number)  , 

group  I  sub-group  |  Filters j  f  Passive  Filters 

Multi-Dimensional  Filtering^  Lossless  Filters 

Signal  Processing  . 


19.  ABSTRACT  (Continue  on  rovarte  if  nocemary  and  identify  by  block  numbsr) 

New  classes  of  multivariable  polynomials  arising  in  studies  of  passive  multidimensional 
systems  have  been  identified  and  their  properties  have  been  studied.  The  problem  of 
structurally  passive  synthesis  of  multidimensional  digital  filters  as  a  cascade  inter¬ 
connection  of  more  elementary  building  blocks  has  been  addressed  via  the  factorization 
®*trix .  The  problem  of  sampling  rate  alteration  of  deterministic  multidimensional  signals 
is  addressed  on  the  basis  of  frequency  domain  description  of  the  signal.  Experiments  with 
digitized  images  are  also  reported  to  demonstrate  the  performance  of  the  designed  inter¬ 
polation  and  decimation  schemes.  - - 

7  (  f  f  /  t 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
Cfl  UNCLASSIFIED/UNLIMITEO  □  SAME  AS  RPT 


22a  NAME  OF  RESPONSIBLE  INOIVIOUAL 
Richard  N.  Smith 


DO  form  1473.  JUN  K  i 


121  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22b  TELEPHONE  (fncfuS*  Ar*«  Cod*)  22c  office  SYMBOL 
(315)  330-3224  RADC 


Previous  editions  are  obsoltt*. 


JRITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


INDEX 


1.  Introduction 


page 

1 


2.  Simple  proofs  of  stability  results 

on  multidimensional  polynomials  .  6 

3.  New  results  on  stable 

multidimensional  polynomials  . .  38 

4.  Realization  of  structurally 

passive  multidimensional  digital  filters  .  81 

5.  Multidimensional  interpolation  and 

decimation  schemes  for  fast  implementation  .  100 

6.  Conclusions .  121 


Figures. 


130 


CHAPTER  X 


INTRODUCTION 


The  necessity  of  dealing  with  signals  and  systems  having  two 
or  more  independent  parameters  is  now  well  recognized. 
Specific  areas  where  multidimensional  signal  processing  find 
applications  include  image  processing,  remote  sensing,  target 
tracking,  robotics,  geophysical  and  astronomical  data 
processing.  A  wide  variety  of  problems  belonging  to  these 
categories  can  be  dealt  with  either  by  the  methods  of 
multidimensional  frequency  selective  filtering  or  by 
employing  multidimensional  modelling  schemes.  For  most 
practical  applications  two  specific  types  of  multidimensional 
filters  are  of  interest.  To  exemplify  the  situation  it  may  be 
recalled  that  in  image  processing  type  applications  the 
processing  scheme  is  mandated  to  be  rotation  invariant  which 
requires  that  pass/stop  regions  of  the  frequency  response  of 
the  filter  be  spherically  symmetric  in  the  multidimensional 
frequency  hyperspace.  On  the  other  hand,  in  the  problem  of 
discriminating  image  regions  moving  with  different  velocities 
(e.g.,  a  target  against  a  background)  the  dynamic  scene  can 
be  modelled  as  a  time  dependent  intensity  distribution 
s(x,y,t)  or  its  Fourier  transform  S(u,v,f)  where  u,  v  are  the 
spatial  frequencies  and  f  is  the  temporal  frequency.  We  then 
have: 


S(x,y, t)  -  m(x-x0-cxt,y-yQ-cyt)  + 

(l-X(x-x0-cxt,y-y0-cyt) J . b( x-dxt , y-dyt ) 

(1.1) 


where  m(x,y)  and  b(x,y)  are  textural  functions  of  the  object 
and  the  background;  X(x,y )  is  a  mask  function  whose  value  is 
one  in  regions  where  the  object  is  present  and  zero  in 


regions  where  no  object  is  present;  (xQ,y0)  is  the  coordinate 
of  the  sensor.  Considering  the  fourier  transform  of  (1.1)  we 
then  have: 

S(u,v,f)-  M(u,v,f)  +  B(u, v, f )  -  A(u,v, f ) *B(u, v, f )  (1.2) 

where  *  denotes  3-D  convolution.  Analysis  of  this  equation 
reveals  that  the  spectral  energies  of  image  regions  moving 
with  different  velocities  are  concentrated  in  different 
planes  in  the  3-D  (u,v,f)  frequency  space.  Therefore,  they 
can  be  discriminated  by  linear  spatio-temporal  filters  the 
frequency  responses  of  which  are  fitted  to  different 
velocity  planes. 

Furthermore,  due  to  possible  interrelations  and  redundancies 
in  the  large  amount  of  image  data  one  needs  a  spatio-temporal 
3-D  stochastic  model  of  the  data  field,  where  again  two  of 
the  dimensions  indicate  the  space  coordinates  and  the  third 
dimension  represents  the  time  coordinate.  The  need  for 
estimation,  prediction,  noise  filtering  etc.  associated  with 
problems  such  as  displacement  estimation,  movement 
compensated  prediction  of  time  varying  images  thus  become 
apparent. 

Considerations  other  than  those  mentioned  above  are  also 
highly  important  in  the  successful  operation  of  a 
multidimensional  recursive  digital  filter  when  realization 
in  hardware  is  sought.  These  are:  (i)  insensitivity  to 
coefficient  pertubation  i.e.,  numerical  stability  to  counter 
the  effects  of  inevitable  rounding  and  truncation  of  signal 
values  (ii)  fault  tolerance  i.e.,  insensitivity  of  the 
overall  filter  performance  to  sudden  and  unexpected  faults  in 
digital  circuitry  (iii)  feasibility  of  hardware 
implementation  in  currently  emerging  high  speed  architectures 
e.g.,  the  systolic  architecture  and  finally,  (iv)  the 
property  that  it  can  be  conveniently  made  adaptive  when  the 


signal  characteristics  vary  in  the  spatial  or  in  the  temporal 
domain. 

The  class  of  one-dimensional  digital  filters  which  has  been 
found  to  satisfy  all  of  the  desirable  characteristics 
mentioned  above  can,  from  a  fundamental  point  of  view,  be 
broadly  categorized  as  the  structurally  passive  filters. 
Structurally  passive  digital  filters  are  those  which  are  not 
only  passive  (i.e,  nonenergy  generating  in  a  discrete  sense) 
from  the  input  output  point  of  view,  but  the  most  elementary 
building  blocks  which  constitute  the  filter  structure  are 
also.  The  fact  that  structurally  passivity*  in  addition  to 
being  responsible  for  the  properties  of  numerical  stability, 
fault  tolerant  and  adaptivity,  can  also  be  exploited  towards 
the  goal  of  implementing  the  filter  in  highspeed  VLSI 
structures  is  now  known  in  the  one-dimensional  context.  The 
well  known  lattice  filters,  wave  digital  filters  and  the 
orthogonal  filters  belong  -to  this  class  of  filters. 

Steps  towards  the  analysis  and  design  of  multidimensional 
structurally  passive  digital  filters  have  been  taken  in  the 
investigation  reported  here.  Since  an  appropriate  transform 
domain  description  of  multidimensional  passive  (or  lossless) 
linear  shift  invariant  quarter  plane  filters  requires  a 
proper  notion  of  stable  multidimensional  polynomials 
previously  not  considered  in  the  literature  our  investigation 
starts  from  a  reexamination  of  multidimensional  stability 
concepts  from  a  fundamental  standpoint.  In  chapter  2  it  is 
shown  that  all  of  the  known  results  on  stable 
multidimensional  polynomials  can  be  derived  via  the 
elementary  artifice  of  a  continuity  property  of  zeros  of  a 
polynomial  as  a  function  of  its  coefficients.  Several 
previously  unrecognized  classes  of  multidimensional  stable 
polynomials  crucial  to  the  transform  domain  description  of 
multidimensional  passive  continuous  systems  are  identified 
for  the  first  time  in  chapter  3  by  utilizing  analytical 


techniques  which  are  largely  similar  to  those  used  in  chapter 
2.  Certain  theoretical  results  on  the  analysis  of  passive 
multidimensional  continuous  systems  are  also  derived  here. 
Synthesis  of  multidimensional  structurally  passive  recursive 
digital  filters  directly  in  the  discrete  domain  form  the 
contents  of  chapter  4.  Motivated  by  one-dimensional  examples 
such  as  the  digital  lattice  filters  and  its  potential 
desirability  in  pipelined  implementation,  the  further 
topological  constraint  that  the  filter  consist  of  cascade 
interconnection  of  elementary  passive  building  blocks  is 
imposed  in  the  problem  formulation.  Necessary  and  sufficient 
conditions  for  the  feasibility  of  such  synthesis  are  derived. 
Since  the  design  problems  of  many  of  the  known 
one-dimensional  structurally  passive  filters  can  be  viewed  as 
special  cases  of  the  results  developed,  new  algorithms  for 
one-dimensional  design  fall  out  as  a  byproduct  of  this 
discussion. 

More  specific  practical  questions  of  image  interpolation  and 
decimation  are  examined  in  chapter  S  by  formulating  the 
problems  in  terms  of  frequency  domain  digital  filtering. 
Certain  computational  FIR-type  structures  are  derived  which 
exploit  full  advantages  of  combined  pipelineability  and 
parallelism  when  high  speed  hardware  implementation  of  the 
resulting  filters  are  sought.  Experiments  with  real  as  well 
as  synthetic  two-dimensional  data  are  reported. 

Conclusions  and  recommendations  for  further  research  are 
chalked  out  in  chapter  6. 

Each  of  the  following  chapters  are  self  contained  and  can  be 
read  independently.  For  similar  discussions  in  the  open 
literature  we  refer  to  the  publications  in  [1],  [2],  [31  and 
(41  in  the  following. 
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CHAPTER  2 


SIMPLE  PROOFS  OF  STABILITY  RELATED 
PROPERTIES  OF  MULTIDIMENSIONAL  POLYNOMIALS 


2.1. Introduction 


The  criterion  for  bounded-input-bounded-output  property  of 
multidimensional  (k-D)  linear  shift  invariant  systems,  when 
the  rational  transfer  function  associated  with  the  system 
does  not  have  non-essential  singularities  of  the  second  kind 
on  the  distinguished  boundary  of  the  polydomain  under 
consideration,  is  well  established  [1-5] .  More  specifically, 
if  the  transfer  function  of  a  k-dimensional  discrete  time 
system  is  given  by  H(z)  as  in  (2.1),  where 

H(z)  -  A(z)/B(z)  (2.1) 

A(  j! )  ■  A  (  z^ ,  ,  z^  )  and  B(^)  *  B(z^,z^,...,z^)  are 

relatively  prime  polynomials  in  the  k  variables  z  « 
( z^ , Zj * • . • » 2^ )  then  under  the  restrictive  hypothesis  that 
A(  3! )  and  B(z)  do  not  have  any  common  zero  on  the 
distinguished  boundary  Iz^  -  1,  i  -  1  to  k  (also  to  be 

denoted  as  |z|  •  1  in  the  forthcomming  discussion)  of  the 
polydomain  1 z ^ |  <  1,  i  -  1,2, ...,k  (also  to  be  expressed  as 

1^1  £  1),  the  system  produces  a  bounded  output  in  response  to 
a  bounded  input  if  and  only  if  B(z)  is  a  strict  sense 
Hurwitz  polynomial,  i.e.,  (2.2)  is  satisfied. 

B(z)  jt  0  for  |z|  <  1  (2.2) 

Since  for  a  given  polynomial  B(z)  it  is  not  possible  to  test 
for  condition  (2.2)  directly,  a  number  of  alternative  but 
equivalent  conditions,  which  are  easier  to  test  for,  have 
been  derived  by  Strintzis,  DeCarlo,  et.  al.  [7]  and  others. A 


simaary  of  all  these  results  are  available  in  the  works  of 
Jury  (8J  and  Bose  (1J.  A  variety  of  different  methods  of 
proofs  of  these  results  have  appeared  so  far  in  the 
literature.  Strintzis  in  (6]  uses  analytic  function  theory, 
DeCarlo  et.  al.  uses  homotopy  theory  expounded  in  Rudins  book 
(91.  It  has  been  pointed  out  (1,101  that  the  results  just 
mentioned  can  be  derived  as  special  cases  of  Rudin's  theorem 
(Theorem  4.7.2,  pp.87  in  (9]).  Delsarte,  Genin  and  Kamp  [10] 
have  shown  that  all  these  results  including  Rudin's  theorem 
can  be  proved  via  a  number  of  elementary  one-dimensional 
( 1-D )  steps.  However,  the  proofs  given  in  [10]  still  require 
the  use  of  some  function  theoretic  results,  which  may  be 
inaccessible  to  an  engineering  reader.  The  present  report 
deals  with  proofs  of  the  above  results  ,  which  are  very 
simple  and  highly  intuitive  as  well.  The  technique  dwells  on 
the  fact  that  the  zeros  of  a  polynomial  can  be  viewed  as 
continuous  functions  of  its  coefficients.  A  correct  and 
complete  statement  of  this  continuity  property,  which 
includes  all  possible  degenerate  cases  (although  restricted 
to  polynomials  in  one  variable),  is  available  e.g.  in  [11] 
(Theorem  4,  p.19,  including  footnote),  [14]  (§4.4),  and  [14a] 
(p.200).  The  technique  has  already  proven  to  be  a  very 
powerful  tool  in  recent  studies  on  passive  multivariable 
network  theory  [12],  A  similar  effort  in  this  direction  is 
noted  in  [13].  Howev.r,  [13]  deals  with  the  two-variable  case 
(k  -  2)  only,  and  a  fully  appropriate  discussion  of  results 
including  all  possible  degenerate  cases  are  not  given. 

A  complete  statement  of  the  continuity  property  of  zeros  of  a 
multivariable  polynomial  as  a  function  of  its  coefficients, 
which  includes  all  the  degenerate  cases,  is  given  in  Section 
2.2.  Section  2.3  contains  proofs  of  the  main  theorems  of 
Anderson,  Jury  [16]  and  those  due  to  Strintzis  [6],  DeCarlo 
et  al.  [71.  Another  result  originally  proved  by  DeCarlo, 
Murray  and  Saeks  [7]  forms  the  main  topic  of  discussion  in 
Section  2.4,  where  it  is  again  shown  that  all  related  results 


rasanaFasrsnanra 


can  be  derived  from  the  continuity  property  of  zeros  of  a 
polynomial  as  a  function  of  its  coefficients.  In  section  2.5 
it  is  shown  how  the  results  of  section  2.3  and  2.4,  when  used 
with  further  continuity  type  arguments  lead  up  to  Rudin's 
theorem  mentioned  earlier.  At  this  point  a  few 
generalizations  of  Rudin's  theorem  are  also  proven  as 
consequences  of  discussions  of  earlier  sections  of  the 
report.  Finally,  the  report  is  summarized  and  conclusions  are 
drawn  in  Section  2.6. 


2.2.  Continuity  property  of  the  zeros  of  a  polynomial  as  a 
function  of  its  coefficients. 


Due  to  the  fact  that  the  utilization  of  the  aforementioned 
continuity  property  in  stability  related  problems  in  the 
context  of  multi-dimensional  polynomials  is  the  main 
contribution  of  this  report  and  that  a  proper  formulation  of 
the  property  needed  for  our  purposes  has  not  until  recently 
been  known  to  be  available ,  at  least  not  in  the  engineering 
literature,  we  undertake  to  give  a  brief  exposition  of  the 
result  without  giving  any  details  of  proofs. 


The  one  variable  version  of  the  following  result  is  well 
known  (11]  in  mathematical  literature.  For  a  proof  see,  for 


example,  [14].  Let  g(z)  -  I  At-  be  a  polynomial  in  z  - 


(z. ,z,, . . . ,zu) ,  where  N  deXS^es  the  set  of  k-tuples  \>  such 

1  *  K  2  2  2 
that  is  not  zero.  Also  let  ||z||  •  |ZjJ  ♦  |z2|  +  —  + 

| z^ | 2 .“Then,  if  the  coefficient  Av  are  assumed  to  be  variable 


quantities  with  certain  initial  values  Av  ,  and  gQ(z)  be  the 
corresponding  expression  for  g(z),  the  following  two  mutually 
exclusive  cases  can  arise: 


(1)  gQ(z)  is  identically  zero,  i.e.,  A  ■  0  for  all 
Vq  c  N . 


(2)  g  ( z )  0  0,  i.e.,  there  exists  at  least  one  v  e  N 


such  that  A  00. 

-0 


If  in  this  latter  case  there  exists  a  zQ  with  finite  IIZqII 
such  that  9q(!q)  ■  0,  then  to  any  c  >  0  we  can  make 

<  S  ,  for  all  v  c 


correspond  a  b  >  0  such  that  for  |A^  -  A^ 

N,  there  exists  a  value  of  z  for  which  we^have  ||z  - 


-0 


c  as  well  as  g(z)  »  0.  The  proof  of  the  above  result  follows 
from  its  one-variable  counterpart,  and  is  available  in  [15]. 


t-.e  u  .  £ 


For  the  present  purpose,  however,  a  different  formulation  of 
the  above  principle  proves  to  be  more  useful  and  is 
elaborated  upon  next.  Assume  the  A^  for  each  v  e  N  to  be 
continuous  functions  of  some  vectoT  t,  say,  Av  «  Av(t). 
Consequently,  we  can  write  g  -  q(z,t) .  Consider~f irst~a  fixed 
value  t  ■  tg  of  t  such  that  g(£Q#tQ)  -  0  for  some  z  -  ssq. 
Then  one  of  the  following  two  mutually  exclusive  cases  must 
hold  true: 


(1)  Av(tQ)  -  0  for  all  v  c  N,  i.e.,  g(z#t0)  -  0  for  all 

£  ■  (  2  ^  »  Z 

(2)  A ( Jt q )  /  0  for  at  least  one  v  e  N. 


In  the  latter  case,  if  we  move  t  along  a  continuous  curve 
from  its  initial  value  tQ,  it  is  possible  to  move  also  each 
2i#  i  -  1  to  k  along  certain  continuous  curves  in  their 
respective  z^  -  planes  in  such  a  way  that  g(z,j:)  ■  0 
continues  to  hold  true,  until  a  value  t  ■  tj  of  t  is  reached 
sjch  that  either  A^tt^)  -  0  for  all  v  e  N  or,  for  t 
approaching  tf,  Iz^  *♦  •  for  at  least  one  i  in  1  <  i  <  k. 
For  k»l,  this  result  amounts  simply  to  a  reformulation  of  a 
corresponding  result  in  (14a]  ($39),  and  it  can  be  easily 
extended  to  the  case  k>l  by  using  the  same  simple  approach  as 
in  Appexdix  1  of  (IS). 
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2.3.  Simplified  proofs  of  results  of  Strintzis  [6)  and 
DeCarlo  et.  al.  17]  and  Anderson  and  Jury  C 16 ] 


We  will  first  need  the  following  lemma. 

Lemma  2.1:  Let  f(z)  ■  f ( z^ , , . . • * zk )  be  a  polynomial  in  z  « 
(z1,z2, . . . ,2^) .  Assume  that  there  exists  a  zQ  ■ 

^ZX0'Z20' *  *  * ,2k0^  in  1-0*-  1  such  that  £(£q)  -  0,  and  |ziQ|  - 
1  holds  for  v  of  the  k  variables  z where  v  is  any  integer 

0  <  v  <  k-1.  Then  there  exists  a  z£  -  ( z{q • z20 •  *  *  *  * zk0 *  in 

\z'q  |  <  1  such  that  f(z£)  -  0  and  |  z£Q  |  -  1  holds  for  at  least 
(v+1)  of  the  z£q.  In  particular,  such  a  zero  then  exists 
with  z£q  -  z^q  for  all  those  i  for  which  |ziQ|  ■  1* 

Proof :  No  loss  of  generality  occurs  in  assuming  that  | z i Q  |  - 
1  for  i  -  1  to  v  .  Let  us  freeze  the  variables  z^  at  z^  -  z^q 
for  i  -  1  to  v.  If  f(jz)  is  independent  of  at  least  one  of  the 
variables  zv+^  to  z^,  say  z^,  then  the  proof  is  immediate  by 
choosing  zfQ  ■  ziQ  for  i  ■  1  to  k-1  and  any  z£Q  with  | z£Q I  ■ 
1. 

Next,  assume  that  f(z)  involves  {i.e.,  actually  depends  on) 
zk-l  *nd  zk*  Let  us  now  freeze  the  variables  z ^  at  z ^  »  ziQ 
for  i  ■  1  to  (k-2),  where  thus  |z^q|  ■  1  for  i  -  1  to  v  and 
I z ^ q I  <  1  for  i-  (v+1)  to  (k-2),  and  consider  a  continuous 
path  from  zJc_1  Q  leading  up  to  the  unit  circle  in  the  zk 

plane.  As  z^^^  is  moved  continuously  along  the  variable 

z^  can  be  moved  along  continuous  path  starting  from  z^q 
such  that  f(z)  -  0  remains  satisfied.  The  only  exception  that 
could  arise  is  that  for  some  zJj_1  -  on.  with 

|«k_il<l#  the  polynomial  f^  defined  by 
f (z10, . .  .Zjt_2  o,zk-l'zk*  *■*  2ero  f°r  *11  Zjj?  the  proof  is 
then  completed  by  choosing  zfQ  -  ziQ  for  1  <  i  <  k-2,  z£^  0 

-  *^.1*  and  any  z£q  such  that  [z£q|  -  1.  If  however  no  such 

exists,  then  by  invoking  the  continuity  property  of 
zeros  of  a  polynomial,  mentioned  in  Section  2.2,  it  follows 
that  either  z^^  or  z^  (or  both)  will  reach  the  unit  circle 


in  the  corresponding  z  plane.  We  would  then  have  constructed 
a  zero  z£  of  f(z)  such  that  l*iol  -  1  for  at  least  (v+1)  of 
the  k  variables  z^. 

We  are  now  in  a  position  to  prove  the  theorem  due  to 
Strintzis[6]  and  DeCarlo  et.  al.(7]  mentioned  earlier. 

Theorem  2.2;  If  f(z)  is  a  polynomial  in  z  -  (z^,z2, . . -z^ ) 
then  f(z)  /  0  in  |z|  <  1  if  and  only  if  the  following 

conditions  simultaneously  hold  true: 

(a)  f(z)  /  0  for  |zj  -  1,  i .  e .  # 1 z^ 1  -  1  for  i  ■  1  to  k. 

(b)  if  j  is  any  integer  in  1  <  j  <  k,  then  f(z)  j*  0  for 

zi"Yij'  *  "  1'2'  *  •  * » (  *  •  •  •Jc*  and  l*jl  *  1*  where  the 

r^j's  are  some  complex  numbers  with  |r^j|al. 

Proof:  Necessity  of  the  theorem  is  obvious.  To  prove 
sufficiency  let  us  assume  that  there  exists  azQ- 
( Zi0,z20, . . . /Zjjq)  with  f(zQ)  -  0,  l*i0l  S  1  for  1-1  to  In 
view  of  lemma  2.1  we  may  assume  that  (k-1)  of  the  k  variables 

(z^Zj . zk)  are  located  on  the  unit  circles  in  the 

corresponding  z^-planes.  A  renumbering  of  the  variables,  if 
necessary,  will  show  that  no  loss  of  generality  occurs  in 
assuming  |ziQ|  •  1  for  i  -  1  to  (k-1).  Due  to  condition  (a) 
of  the  theorem,  we  then  have  <  1-  This  latter 

conclusion  ,  however,  implies  in  particular  that  f(z)  is  not 
independent  of  z^.  Next  we  successively  move  the  variables 
z^,  i-1  to  (k-1)  from  ziQ  to  YiJt  continuously  along  the  arcs 
of  the  corresponding  unit  circles.  The  variable  z^  may 
then  be  moved  along  a  path  starting  from  z^q  such  that 
f(z)  -  0  remains  satisfied.  In  this  process  it  is  impossible 
to  have  f(z)  -  0  for  some  z^  -  ri0  on  r.,  1  <  i  <  k-1,  and 
all  z^  because  otherwise  a  contradiction  with  condition  (a) 
of  the  theorem  is  easily  arrived  at  by  choosing  zk  at  any 
point  on  the  unit  circle  Iz^t  -  1.  Invoking  the  continuity 
property  of  the  zeros  of  a  polynomial,  it  then  follows  that 
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rk  can  be  chosen  to  be  a  continuous  path  in  the  z^-plane  and 
the  process  described  above  can  only  have  one  of  the 
following  two  outcomes.  Either  we  reach  2i  "  Yik'  1  "  1  to 
(k-1)  with  zk  remaining  such  that  |zk|  <  1  or  the  path  rk 
described  by  zk  has  at  least  one  point  in  common  with  the 
unit  circle  | zk |  •  1.  In  the  former  case,  the  condition  (b) 
of  the  theorem  is  violated,  whereas  in  the  latter  case  the 
condition  (a)  of  the  theorem  is  violated.  The  proof  of  the 
theorem  is  thus  complete. 

Remark:  The  proof  does  not  assume  that  f(z)  is  a  polynomial 
or  that  it  is  a  holomorphic  function.  All  that  is  required  is 
that  the  continuity  property  presented  in  Section  2.2  holds. 

A  similar  remark  holds  for  other  results  in  this  paper. 


We  next  need  to  prove  a  result  due  to  Anderson  and  Jury  [16] 
via  the  continuity  argument.  We  need,  however,  the  following 
lemma,  which  will  also  prove  to  be  useful  in  other 
developments  to  follow. 

Lemma  2.3:  Let  f(z)  be  a  polynomial  in  z  -  (z2,z2, . .  .z^)  and 
let  the  set  of  indices  i  -  1  to  k  be  the  union  of  two 
disjoint  subsets  1^  and  Ij.  Then  f(z)  *  0  in  |z|  <1  if  and 
only  if  the  following  conditions  hold  true  sumultaneously : 


(a)  f(z)  #  0  if  z^  -  a^  for  i  e  1^;  |z^|  <  1  for  i  c  I2  where 
the  a^,  for  all  i  t  1^,  are  some  complex  numbers  such  that 
1^1  i  1. 

(b)  f(z)  *  0  if  IZjJ  <  1  for  i  c  1^ ;  |  z  ^  |  -  1  for  i  c  Ij. 

Proof :  Necessity  is  obvious.  To  prove  sufficiency,  let  it  be 
assumed  for  contradiction  that  f(z)  has  a  zero  at  z  ■  Zg  » 

(  z10 '  z20 '  *  ’  '  zk0  )  w*-th  I£qI  S  !•  By  virtue  of  Lemma  2.1,  no 
loss  of  generality  occurs  in  assuming  that  there  exists  a 
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fixed  integer  p  such  that  |  r  ,;j  |  £  1  and  that  |z^g|  *  1  for 
all  other  i  t  (1^  U  ij).  In  view  of  condition  (b)  of  the 
present  lemma,  we  must  have  in  fact  | I < I  and  p  e  I 2 •  We 
move  the  variables  z^  for  i  e  1^  from  their  initial  locations 
*i  -  zig  on  1 2^1  “  1  along  continuous  curves  1*^  lying  in  the 
respective  unit  disc  |z^|  <1  and  leading  up  to  the  terminal 
points  z^  -  a^,  while  the  variables  z^  for  i  e  I2,  i  *  u  are 
held  frozen  at  their  corresponding  values  z^  •  ziQ.  The 
variable  z  can  then  trace  out  a  contour  T  in  the  z  -plane 

p  p  p  r 

so  that  f ( je )  •  0  is  satisfied.  Note  first  that  as  long  as  z^ 
is  on  Tl  for  all  i  c  and  z^  -  ziQ  for  i  c  l2\{;/},  the 
polynomial  f(z)  cannot  be  zero  for  all  z^  because  otherwise 
in  arbitrary  choice  of  z^  on  |z^|  >1  yields  a  value  of  the 
k-tuple  z  with  |zi|  <1  for  i  t  and  | * i |  •  1  for  i  el2 
such  that  f(z)  «  0,  which  contradicts  condition  (b)  of  the 
present  lemma.  However,  with  this  alternative  excluded,  it 
would  be  possible,  due  to  condition  (a)  and  the  continuity 
argument,  to  arrive  at  the  same  unpermitted  situation  in  a 
way  in  which  z^  reaches  a  point  with  |z^|»l  by  moving 
continuously  along  r while  the  z^,  for  i  t  1^,  remain  on 
their  respective  r^. 

Remark :  If  k-2  and  1^(1}  and  I2«{2]  then  lemma  2.3  coincides 
with  a  result  well  known  (see  e.g.,  [2]  and  (13b])  in  the 
literature . 

Corollary  2.3.1:  Let  f(z)  be  a  polynomial  in  z  » 

(z^Zj#  •  •  •  #*k*  *  Tben  f(z)  f  0  for  |z|  <  1  if  and  only  if  the 
following  two  conditions  hold  true  simultaneously: 

(a)  f (a,z2»2j, .  •  •  #zJl)  ^  0  for  |z^|  <1,  i  •  2  to  k,  where  a 
is  some  complex  number  with  |a|  <  1. 

(b)  f  ( Zj ,  z, ,  —  ,  z^ )  p  0  for  |Zj |  <  1,  and  iz^  »  1,  i-2  to  k. 
Proof :  Follows  from  lemma  2.3  via  the  choice  of  1^  •  {1}  and 


-14- 


mv’.v’av, 


I 2  "  {2»3,...k}. 

The  proof  of  the  theorem  due  to  Anderson  and  Jury  [16] 
mentioned  earlier  is  given  next. 

Theorem  2.4:  If  Hz)  is  a  polynomial  in  z  -  ( z^ , z2 * • • • « zk } 
then  t(z)  ft  0  in  |z|  <  1  if  and  only  if  for  some  complex 
number  a  with  |a|£l  the  following  conditions  simultaneously 
hold  true: 

(1)  f(z)  ft  0  for  IZjJ  £  1,  and  |z^|  -  1,  j  -2  to  k. 

(2)  f(z)  ft  0  for  z2  -  a,  |z2l  £  1#  and  |  |  -  1,  j  • 

3  to  k. 

(3)  f(z)  ft  0  for  zi  -  a,  i  -  1,2,  | z 3 1  <  1,  and  | z ^ |  - 
1,  j«  4  to  k. 


(k)  f(z)  ft  0  for  zi  »  a,  i  -  1  to  (k-1),  |  zk  |  <  1. 

Proof :  Necessity  is  obvious.  To  prove  sufficiency,  let  us 
define  the  n-variable  (1  <  n  <  k)  polynomial  from  f(z)  by 

freezing  each  of  the  first  (k-n)  variables  equal  to  a,  i.e., 
z^  ■  a  for  i  »  1  to  (k-n).  Note  that  Bk  «  f(z)  and  B^^  - 
B, ( Zfc )  •  f (a, . . . ,a,zk) .  We  first  claim  that  Bn  t  0  in  |zil  < 
1,  for  i  -  (k-n+1)  to  k  for  each  n  in  1  <  n  <  k.  The  proof  of 
this  assertion  is  via  induction  on  n. 

Obviously,  due  to  condition  (k)  of  the  theorem,  B^  -  B^z^)  ft 

0  in  |zkl  <  1.  The  assertion  is,  therefore,  true  for  n  -  1 . 

Assume  now  that  our  assertion  is  correct  for  B_  i.e.,  B 

n  n 

-  *k-n+l '  *  ’  ’ zk *  **  ®  *or  I  * i  I  £  1#  i  ■  (k-n+1)  to  k.  Note 
that  this  implies  n+l,***,^k^  m  ^n^^k— n+l,**’*^k^  ^ 

0  for  | z ^ |  <  1,  i  -  (k-n+1)  to  k.  However,  condition  (k-n)  of 
the  theorem  states  that  2k-n' * ' ’ ' zk >  f  0  for  'zk-n'  -  1 

and  I z ^ |  ■  1,  i  »  (k-n+1)  to  k.  The  last  two  conditions,  in 
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view  of  corollary  2.3.1,  imply  that  Bn4^  /  0  for|z^|  <  1,  i 
(k-n)  to  k.  The  proof  of  the  theorem  by  induction  is. thus 

complete . 


2.4.  Proof  of  a  theorem  of  DeCarlo,  Murray,  and  Saeks  [71  and 
it*  extension 

In  this  section  we  undertake  the  proof  of  the  theorem  of 
DeCarlo, et.al.  [7]  based  on  continuity  property  of  zeros  of 
polynomials.  We  note  that  this  theorem  is  not  mentioned  by 
Srintzis  [6],  and  its  proof  in  [7]  makes  use  of  homotopy 
theoretic  arguments.  Subsequently  it  was  pointed  out  [1,10] 
that  the  result  can  be  considered  to  be  a  special  case  of 
Rudin's  theorem  [9].  We  will  need  to  have  the  following  lemma 
as  a  preperation  for  our  proof  of  the  theorem. 

Lemma  2.5:  Let  f(z)  be  a  polynomial  in  z  »  (z^,z2 , ...,zk) 
such  that  f(z)  /  0  for  |z|  -  1,  and  the  set  of  indices  i  »  1 
to  k  be  the  union  of  disjoint  subsets  of  indices  1^,  I2  and 
Ij.  If  f(z)  /  0  in  |ZjJ  <  1  for  i  c  I,,  |z^|  -  1  for  i  e  I2, 
and  z^  »  Yi#  for  i  c  I3#  where  the  r^'s  are  some  complex 
numbers  such  that  | I  •  1,  then  we  also  have  f(z)  y  0  in 
JZjJ  <  1  for  i  c  1^ ;  |z.J  -  1  for  i  e  I£,  and  zi  -  for 

i  c  1^,  where  j  is  any  integer  belonging  to  l3,  I£  ■  I2  u 
{j},  and  1^  -  I3\{j}. 

Note  that  the  lemma  is  meaningful  only  if  I3  is  nonempty  and 
that  we  may  also  assume  1^  to  be  nonempty  since  otherwise  the 
result  is  trivial.  However,  I2  may  be  empty. 

Proof :  Let  it  be  assumed  for  the  purpose  of  contradiction, 

that  f(zQ)  ■  0,  where  zQ  -  ( *j.q * z20 ' *  ’  ’ zk0 *  is  such  that 
I  ziQ  1  i  1  for  i  e  'I1,  |zi(J|  -  1  for  i  c  I£,  and  ziQ  -  Yi  for 
i  c  I 3 .  In  view  of  Lemma  2.1  and  the  fact  that  f(z)  /  0  for 
|z|  -  1  we  may  assume  that  for  one  of  the  i  e  I,,  say  for  i  - 
//,  we  have  |z^Q|<l  and  that  |ziQ|  ■  1  for  all  i  t  I£,  where 

1[  -  Ij\{/i}. 

Consider  the  polynomial  f.(z  ,z.)  obtained  by  freezing  in 

•  ^  J 

f(z)  the  rest  of  the  variables  as  follows:  z ^  -  ziQ  for  i  t 
I£  U  l2  and  zi  -  ri  for  i  c  1$.  Clearly,  z.Q)  -  0, 
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but  by  the  hypothesis  of  the  lemma  we  have  f^(z^,Zj)  /  0  for 
|z^|  <  1,  Zj  ■  Yj  as  well  as  for  |z^|-|Zj|«l.  Consider  then 
a  directed  arc  Tj  of  the  circle  |Zj|  -  1  originating  from  z ^ 

«  ZjQ  and  terminating  at  Zj  -  y j .  Obviously,  for  any  fixed 
Zj  e  I\  the  polynomial  f^(z^,Zj)  cannot  be  zero  for  all  z  ^ 
because  otherwise  a  contradiction  with  the  fact  that  f(z)  * 

0  in  |z|  -  1  would  immediately  be  arrived  at  by  choosing  zfj 
do  be  arbitrarily  located  on  |z^|  -  1.  By  invoking  the 
continuity  argument  we  can  thus  move  Zj  continuously  on 
from  Zjg  towards  Yj  and  simultaneously  move  z y  continuously 
from  z^q  such  that  t^(z  f/,z^)  -  0  remains  satisfied.  We  will 
then  reach  a  situation  either  with  |z^|  -  |Zj|  -  1  or  in 
which  |z^|  <  1,  Zj  »  y j #  which  both  have  been  seen  to  be 

impossible. 

The  following  result  is  immediately  obtained  by  repeated 
application  of  Lemma  2.5. 

Lemma  2.6:  Let  f(z)  be  a  polynomial  in  z  -  ( z^ , z2 , . . . , z^ ) 
such  that  f(z)  /  0  for  |z|  -  1,  and  the  set  of  indi  as  i  «  1 
to  k  be  the  disjoint  union  of  two  subsets  and  I2.  Then  if 
f(z)  *  0  in  |z^|  <  1  for  iel^  and  z^^  -  y^  for  iel2  where  the 
y^'s  are  some  complex  numbers  such  that  |  YjJ  -  1,  then  f(z)  * 
0  in  |ZjJ  <  1  for  i  c  1^,  and  | z^ |  *  1  for  iel2. 

Note  again  that  the  set  result  is  trivial  if  I2  is  empty. 

Lemma  2.7:  The  transformation 

Zj^  »  (u-v)/(  l-v*u) ,  z2  •  (u+v)/(  l+v*u)  (2.3a,b) 

has  the  following  properties: 

Property  2.1:  If  v  -  0  then  z^  -  z2  -  u. 

Property  2.2:  If  |u|  -  1  and  v  is  such  that  v*  ?  +l/u,  then 

I  I  ■  I z2 1  *  1  * 
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Property  2.3:  If  |u|  <  1  and  |v|  <  1,  then  |  |  <  1  and  |z2l<l. 

Property  2.4;  For  any  z^  and  z2  with  |z^|<l  and  |z2(<l  there 
exist  u  and  v  with  |u|  <  1,  |v|  <  1  such  that  (2.3)  is 
satisfied.  (Note:  (2.3)  does  not  represent  a  simple 
one-to-one  transformation,  but  Property  2.4  is  in  a  sense  the 
converse  of  Property  2.3.) 

Proof :  Properties  2.1  and  2.2  are  easily  verified.  For 

proving  Property  2.3,  consider  any  fixed  v  in  |v|  <  1.  Then 
z^  is  an  analytic  function  of  u  in  the  domain  |u|  <1. 
Furthermore,  from  property  2.2,  )z^|  -  1  for  all  u  on  the 
boundary  |u|  -  1  of  the  domain  |u|  <1.  Thus,  maximum  modulus 
theorem  implies  that  |z^|  <  1  for  |u|  <  1.  Similar  arguments 
hold  for  z2< 

For  proving  Property  2.4,  note  first  that  if  z^^  »  z2  then  the 
proof  immediately  follows  by  choosing  u  «z^  -  z2  and  v  -  0. 
Henceforth  zj/z2  will  be  assumed.  By  eliminating  u  from 
(2.3a)  and  (2.3b)  it  follows  that: 

(l-|v| )2/2|v|-(z1+exp( ja) ) ( z2«xp(-ja)-l )/( z2-z2 )  (2.4) 

where  v- | v | exp( ja) . 

We  next  claim  that  by  choosing  a  properly  the  right  hand  side 
of  (2.4)  can  be  made  to  be  equal  to  a  finite  real  and 
positive  number.  To  substantiate  this  claim  consider  the 
angle  y  defined  by  (  z^expf  jo) )  (  z2exp( -ja)-l )«  |c|exp(jr). 
Since  Iz^l  <  1  and  |z2|  <  1,  it  is  obvious  that  |c|  /  0,  and 
r  is  thus  well  defined.  Furthermore,  y  satisfies  (2.5). 

exp(  2  j  r )  -  [(^z1)/(l+z1*5)U(C-*2)/(l-z2*OJ  (2.5) 

where  K  -  exp( ja) 

Since  the  right  hand  side  of  (2.5)  is  the  product  of  two 
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allpass  functions  of  the  variable  £  it  is  well  known  that  as 
a  increases  continuously  by  2it,  the  angle  associated  with 
each  of  the  two  factors  also  increases  continuously  by  an 
amount  2n.  Consequently,  r  then  also  increases  continuously 
by  an  amount  2it,  i.e.,  by  proper  choice  of  a,y  can  be  given 
any  arbitrary  value.  The  proof  of  our  claim  then  follows 
from  the  fact  that  r  is  the  angle  associated  with  the 
numerator  of  the  right  hand  side  of  (2.4),  and  that  the 
corresponding  denominator  depends  only  on  the  given 
quantities  z ^  and  z2<  By  solving  (2.4)  for  |v|  it  then 
easily  follows  that  one  of  the  two  solutions  satisfies 
0< | v | <1 . 

Finally,  since  (2.6)  follows  from  (2.3),  by  using  arguments 
similar  to 

u  -  { z1+v)/( l+v*zx )  -  (z2-v)/(l-v*z2)  (2.6) 

those  used  in  the  proof  of  property  2.3,  it  follows  that  |v| 
<  1  together  with  either  |z.|  <1  or  |z,|  <  1,  implies  that 

lul  <  1. 


We  prove  the  theorem  due  to  DeCarlo  et.  al.  (7]  next. 

Theorem  2.7:  Let  f(z)  be  a  polynomial  in  z  -  ( z^ , z2 , . . . , z^  )  . 
Then  f(z)  /  0  in  |z|  <  1  if  and  only  if  the  following 

conditions  hold  true  simultaneously: 


(a)  f(z)  *  0  in  |*|  -  1 

(b)  f(z)  ^  0  in  |z|  <  1,  where  z^  -  z  for  each  i-1  to  k 


First  proof  of  Theorem  2.7 


We  first  provide  a  particlarly  simple  proof  of  theorem  2.7 
for  the  two-variable  case  i.e.,  when  k-2.  Subsequently,  it 
will  be  shown  via  induction  on  k  that  if  the  theorem  is  true 


for  k-2,  then  it  is  true  for  any  k.  For  all  this,  no  use  will 
be  made  of  the  results  obtained  in  Section  2.3. 


Only  sufficiency  needs  to  be  proved,  the  necessity  being 
obvious . 

(i)  k«2. 

Consider  the  function  h(u,v)  in  (2.7)  obtained  from 

via  the  transformation  (2.3),  where  n^ ,  and  n2  are  degrees  of 

f(z^,Z2)  in  z ^  and  Z2  respectively. 

h(u,v)  -  (l-v*u)  1(l+v*u)  2  ftZj^Zj)  (2.7) 

The  function  h(u,v)  is  not  a  polynomial  in  u  and  v,  but  may 
be  considered  as  a  polynomial  in  u  whose  coefficients  are 
polynomials  in  v  and  v*.  The  coefficients  just  mentioned  are, 
therefore,  continuous  functions  of  v. 

Property  2.1  of  Lemma  2.7  along  with  condition  (b)  of  theorem 
2.7  then  yields: 

h(u,0)»  f(u,u)/  0  for  |u|  <  1  (2.8) 

Furthermore,  it  follows  from  property  2.2  and  condition  (a) 
of  tneorem  2.7  that: 

h(u,v)  /  0  for  |u|  -  1,  | v |  <  1  (2.9) 

We  next  claim  that  condition  (2.8)  and  (2.9)  imply  that 
h(u,v)  0  for  |u|  <  1  and  |v|  <  1.  To  show  this  we  assume 

for  contradiction  that  h(uQ,v0)  -  0  with  |uQ|  <  1,  |vQ|  <  1. 
Consider  a  continuous  directed  arc  r  in  the  complex  v-plane 
originating  from  v-v^  and  terminating  in  v-0 .  For  any  fixed  v 
on  ry,  h(u,v)  cannot  be  zero  for  all  u  because  otherwise  a 
contradiction  with  (2.9)  would  be  arrived  at  by  choosing  u  to 
be  arbitrarily  located  on  |u|-l.  Furthermore,  since  as  v 
describes  u  describes  a  continuous  path  in  the  u-plane, 

it  follows  from  (2.8),  by  invoking  the  continuity  argument, 
that  there  must  exist  a  v-Vg  on  I"v  and  a  corresponding  u-Ug 


with  |u£|-l  such  that  h(u£,v£)-0,  which  again  contradicts 
(2.9)/  thus  proving  the  claim  that  h(u,v)^0  in  |u|<  l,|v|<l. 

If  f(z^,Z2)“0  for  some  z^-z^q  and  z2mZ20  *n  lzil  <  *  an<*  lz2^ 

<  1,  then  it  folLows  that  from  (2.7)  and  property  2.4  of 
Lemma  2.7  that  there  exists  Uq  and  Vq  in  |u|  <  1  and  |v|  <  1 
such  that  h{u0,Vg)-O,  which  however  contradicts  the 
conclusions  of  the  last  paragraph.  Thus  f(z.,z,)/0  in  |z.|  < 
1.  |Z2I  <1. 

It  only  remains  to  show  that  ftz^/ZjMO  if  one  of  the  two 
variables  z^  and  z2  is  strictly  inside  and  the  other  on  the 
corresponding  unit  circle.  Assume  e.g.,  ^zi0'z20^“°  f°r 

| z^o I “1 » I z2q I  <  Then  f(z^Q»z2)  cannot  be  zero  for  all  z2 
because  otherwise  an  arbitrary  choice  of  z2  on  |  z2 1  ’■■1  would 
violate  condition  (a)  of  theorem  (2.8).  The  continuity 
property  of  zeros  then  implies  that  by  moving  the  variable  z^ 
from  z^-z^g  to  inside  the  unit  circle  by  an  arbitrarily  small 
amount  it  would  be  possible  to  construct  a  zero  of  t(z^,z2) 
in  IZjJ  <  1 , I z  2 1  <  1#  the  impossibility  of  which  has  been 
demonstrated  in  the  previous  paragraph. 

(ii)  k  >  2. 

Assume  for  the  purpose  of  induction  that  the  theorem  is  true 
for  (k-1)  variables,  with  k  >  2,  and  consider  the  polynomial 
fl(Z,zk)  of  two  variables  zi>z2  def*ned  as  ^(Z/Z^)  - 

f ( z , . . . z , zk ) .  Then  it  obviously  follows  from  condition  (a)  of 
Theorem  2.7  that  0  for  |z|  -  1  and  |  zk  |  -  1, 

whereas  condition  (b)  implies  that  f^(z,z)^0  for  |zf  <  1.  The 
last  two  conditions,  due  to  the  proof  already  given  for  case 
(i),  imply  that  f^z.Zj^O  for  |z|  <  1  and  |  z^  |  <  1.  In 

particular,  f ^ ( z , zk ) -f ( z , . . . z , zk )/0  for  |z|-l  and  |zkl  <  1. 
The  latter  conclusion  along  with  condition  (a)  and  Lemma  2.6 
imply  that  f ( z^ , z2 , . . . zk )/0  for  |  z i I  — 1 ,  i-i  to  (k-1)  and  |zkl 

<  1. 


-22- 


■>  - - 


Consider  next  the  ( k— 1 )  variable  polynoaial  defined  as 

f2(z^,z2*  •  •  m  f(*^#*2'**,xk-l'*k0^'  SfcQ  ^  ®  any 

fixed  complex  number  in  |zk|  <  1.  Zt  then  follows  that 

f2(z,z,...z)  /  0  for  |zj  <  1,  because  otherwise  f^tz.z^g)  - 

f ( z , . . . z , z^g  )  would  have  a  zero  in  |z|<  1,  the  impossibility 
of  which  has  already  been  demonstrated.  Furthermore,  from  the 
concluding  sentence  of  the  last  paragraph  it  follows  that 

*2*zl'z2'  *  *  *zk-l^°  for  I*!!"1*  t0  t***1)*  Therefore,  by 
invoking  the  induction  hypothesis  we  can  assert  that 

£2^*l'z2# ’ * ,zk-l*  *n  *2i*  t0  which  in  view 

of  the  fact  that  z^q  is  arbitrarily  located  in  (z^l  <  1,  in 
turn  implies  that  f(z)^0  for  |z|  £  1. 

Corollary  2.7.1:  The  polynomial  f(z)  in  z  »  (  z^ , z2 , . . . zk )  is 
devoid  of  zeros  in  |£|  <  1  if  and  only  if  the  following 
conditions  simultaneously  hold  true: 

(a)  f(z)  «t  0  for  |  z  |  »1 

(b)  f (£( z ) )  *  0  in  |  2 | <1, 
and  each  g^z),  for  i 
analytic  in  | z | £1  with 
unit  circles  into  unit 
discs . 

Proof :  Let  us  note  that  any  function  g^(z)  which  maps  unit 

disc  into  unit  disc,  unit  circle  into  unit  circle,  and  is 
also  analytic  in  | 2 | <1  is  a  rational  function  which  can  be 
written  [19,  p.12]  as  in  (2.10),  where  the  constants  and 

«l>(  satisfy  IrJ-l,  l«ivl<l. 

zi"^i*ui*"ri  ®  *  ui”*i  v^//^“aivui  * '  vi -1  (2.10) 

Furthermore,  we  also  have  that: 

Iglu^l  «  1  fot  |u.|  •  1  (2.11) 


where  e[(z)-(gL(z),g2(z),...gk(z)), 
*  1  to  k  are  functions  which  are 
the  further  property  that  they  map 
circles  and  unit  discs  into  unit 
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Consider  the  polynomial  g(u)  in  u-(u1#u2, . . .u^) ,  defined  as 
in  (2.12),  where  n,  is  the  partial  degree  of  f(z)  in  z^. 


*  ».  *  ni 

h(u)-f(z)  n  n1  (l-<*.  u.)  1 
”  i-lv-1  1  1 


(2.12) 


Due  to  condition  (b)  of  the  present  theorem  it  then  follows 
that  h(u)*0  in  |uf^l,  where  each  u^»u  for  i-1  to  k. 
Furthermore,  it  follows  from  equation  (2.11)  and  condition 
(a)  of  the  present  theorem  that  h(u)*0  for  | u | —1 .  Invoking 
theorem  2.7  it  then  follows  that  h(u)#0  in  |u|<l. 


Next,  consider  any  z-(z^,z2l 


. . .zk) 


with  | z ) <1 .  For  each 


corresponding  z^  we  can  compute  a  u^  by  means  of  equation 
(2.10)  i.e.,  by  solving  an  algebraic  equation  of  degree  voO. 
Hence,  we  can  find  a  u-(u. ,u2, . . .u^)  such  that  (2.10)  is 
satisfied,  and  due  to  (2.11),  we  then  have  |u|<l.  Hence,  it 
follows  from  (2.12)  that  f(z)*0  for  |z|£l. 


Remark : 


He  note  that  Theorem  2.7  is  a  special  case  of 


Corollary  2.7.1  when  the  choice  ( z )  -  z  for  all  i  -  1  to  k 

is  made. 

The  following  multidimensional  (k>2)  versions  of  results 
stated  in  [xO]  and  (1],  respectively,  follow  immediately  from 
Corollary  2.7.1. 

Corollary  2.7.2:  If  f(z)  is  a  polynomial  in  z  -  ( z^ , 
*2,...'*k)  then  f(z)  /  0  in  |z|  £  1  if  and  only  if  the 
following  conditions  simultaneously  hold  true: 


(a) 

(b) 

where  the 


f  (  z)  i*  0  in  |z|  -  1 
f  ( z )  p  0  in  jz 


<  1  with  z^  -  c^z 


c  v  s  are  some  unimodular  constants 


\  i . 

(i.e., 


•1  to  k 

!cil  “ 


1),  and  the  p^s  are  some  positive  integers  (i.e.,  px  >  0 


Corollary  2.7.3 


Let  f(z)  be  a  polynomial  in  z  - 
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(z1#22, . . . ,ak) .  Then  f{z)  ^  0  in  |z|  <1  if  and  only  if  f(z) 
yi  0  for  |  z  |  <1  with  z^  -  z  exp  (j(5^),  i  -  1  to  k,  where  the 
0^'s  are  arbitrary  real  numbers. 

Remark:  An  alternative  proof  of  theorem  2.7,  which  heavily 
makes  use  of  the  results  developed  earlier  in  this  report, 
but  does  not  necessitate  .the  technique  of  induction  on  the 
number  of  variables,  can  be  formulated  along  lines  explained 
later.  We  first  consider  the  following  theorem. 

Theorem  2.8:  Let  f(z)  be  a  polynomial  in  z  -  ( z., z«, . • . ) • 
Then  f(z)  *  0  in  |z|  <  1  if  and  only  if  the  following 
conditions  simultaneously  hold  true: 

(a)  f(z)  ^  0  in  |z|-l 

(b)  t(z)  *  0  for  Zj  -  a,  and  | z ^ |  <  1  for  i  /  j 
where  a  is  some  complex  constant  such  that  |a|  <1  and  where 
j  is  some  integer  1  <  j  <  k. 

(c)  f(z)  *  0  if  |Zj|  <  1  and  z ^  for  the  indices  i  *  j, 
where  the  r^’s  are  some  complex  numbers  such  that  j ri I  •  1. 

Proof:  Necessity  of  the  theorem  is  obvious.  To  prove 
sufficiency  we  first  note  that  no  loss  of  generality  occurs 
if  we  assume  j  -  1.  In  view  of  Lemma  2.6,  with  1^  •  {1}, 
conditions  (a)  and  (c)  of  the  theorem  together  imply  that 
f (z)  /  0  if  | z ^ |  <  1  and  |z^|  -  1  for  i  -  2  to  k.  However, 
this  latter  conclusion  along  with  condition  (b)  of  the 
theorem,  due  to  Corollary  2.3.1,  implies  that  f(z)  ^  0  in  |z| 
<  1- 

Second  proof  of  Theorem  2.7: 

Consider  a  polynomial  g(u)  in  u  •  ( u^ , u2 , . • • , u^ )  obtained  by 
making  in  the  polynomial  f(z)  the  substitutions:  z ^  -  uk  and 
z.  -  u.u.  for  each  i  -  1  to  (k-1).  Conditions  ( 2 . 1 3 ) , ( 2 . 14  ) 
then  immediately  follow  from  condition  (b)  of  Theorem  2.7, 
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whereas  (2.15)  follows  from  condition  (a)  of  Theorem  2.7  and 
the  fact  that  |u|«l  implies  |z|-l. 

g(u1# . . .  *uJt_1,0)«f  (0,0, . .  .0)/0 

for  all  finite  u^,  i-1  to  k-1  (2.13) 

g(l, — ,l,uk)  "  f^uk' — /Uk^  *  0  for  *uk'  -  1  (2.14) 

g(u1# . . . ,uk)  jt  0  for  |u|  -  1  (2.15) 

Invoking  Theorem  2.8  on  g(u),  with  j-k,  a-0 ,  and  i-1  to 

k-1,  it  then  follows  from  (2.13),  (2.14),  and  (2.15)  that 
g(u)  i*  0  for  |u|  <  1. 

We  next  note  that  |u|<l,  thus  f ( z )»g(u)/0,  if  any  of  the 
following  conditions  hold:  (a)  |  z^  |  <^1  and  |  z i  |  — 1  for  i-2  to 
k,  (b)  Zj^-0  for  i-1  to  j-1,  |z.|£l,  and  |z^|-l  for  i-j+1  to 

k,  where  j  is  any  integer  such  that  2£j£k-l,  (c)  z^O  for  i-1 
to  k-1  and  |zk|<l.  Hence,  due  to  Theorem  2.4,  f(z)  /  0  for 
|z|  <  1,  thus  providing  an  alternative  proof  of  Theorem  2.7. 


2.5.  Proof  of  Rudin's  theorem  and  an  Extension  of  Rudin's 
Theorem 


In  this  section  we  undertake  the  proof  of  Rudin's  theorem 
using  continuity  argument.  A  generalization  of  Rudin's 
theorem  is  also  reported.  Our  proof  of  Rudin's  theorem  is 
based  entirely  on  continuity  arguments  and  an  application  of 
Corollary  2.7.2.  Proofs  previously  reported  in  the  literature 
have  been  obtained  by  appealing  to  homotopy  theory  [9].  Our 
extension  of  Rudin's  theorem  is  a  further  generalization  of 
the  results  discussed  in  (17]  and  involves  a  simple 
application  of  our  lemma  2.3  and  lemma  2.6,  along  with  the 
Rudin's  theorem  itself.  We  first  state  and  prove  the 
conventional  form  of  Rudin's  theorem  for  covenience  of 
exposition. 

Theorem  2.9  (Rudin  [1],[9]):  The  polynomial  f(z)  in  :  • 
( z1 , z2 , . . . zk )  is  devoid  of  zeros  in  |z[  <  1  if  and  only  if 
the  following  conditions  simultaneously  hold  true: 

(a)  f(z)  *  0  for  |z|  -  1 

(b)  f(2(z))  *  0  in  |  z  |  <  1,  where  g[(z)  - 
(g^tz) ,g2(z) , . • .gk(z) ) ,  and  each  g^Cz),  for  i  »  1  to  k  are 
continuous  functions  of  z  with  nonnegative  winding 
numbers  with  respect  to  the  origin  with  the  further  property 
that  they  map  unit  circles  into  unit  circles,  and  unit  discs 
into  unit  discs. 

Note  that  this  theorem  reduces  to  corollary  2.7.1  if  the 
g^(z)'s  are  assumed  to  be  analytic  in  |z|<l;  indeed  this 
added  assumption  ensures  that  g^(z)'s  are  all-pass  functions 
(i.e.,  unit  functions  (19]).  A  proof  such  as  in  [10], 
however,  is  not  fully  correct  because  of  the  fact  that  if  the 
functions  g^z)  are  merely  continuous  the  principle  cf 
argument  may  not  be  invoked.  An  approach  modified  wherever 
required  will  therefore  be  used. 


Proof:  Necessity  is  obvious.  To  prove  sufficiency,  Let  c^  - 
g^l),  and  p.  >  0  b«  the  winding  nuaber  of  g^(z)  with  respect 
to  the  origin  for  all  i  ■  1  to  k.  Also,  let  Tq  and  be  the 


contours  described  in  the  complex 
by  the  functions  f(c^z 


plane  Refined  r 
#  #  •  •  •  C 


f(g1(z) ,g2(z) , . . .gk(z) ) ,  as  the  variable  z  describes  the  unit 


circle  z  -  exp(je)  beginning  froa  0-0  and  ending  in  0  -  2it 


in  the  anticlockwise  direction.  Clearly,  rQ  and  are  closed 
contours  lying  in  the  finite  complex  plane,  each  with  the 


same  initial  and  terminal  points  at  f (c1#c2, . . .ck) . 


He  first  claim  that  the  nuaber  of  encirclements  of  the  origin 
of  the  complex  plane  by  the  contours  rQ  and  are  the  same. 
To  prove  this  assertion  we  define  the  functions  ei(0)  for  0  < 
0  i  2s  as:  ^  -  argtg^expt  j0) )  for  each  i  -  1  to  k.  Note 
that  since  each  g^(z)  is  a  continuous  mapping  of  the  unit 
circle  into  the  unit  circle,  *^(0)'s  are  continous 
functions  of  0,  and  we  can  also  write: 


gi(exp(j0))  -  exp( 0) )  for  i  -  1  to  k  (2.16) 

In  particular,  we  have  (2.17a),  whereas  (2.17b)  follows  froa 
the  fact  that  the  winding  nuaber  of  g^ ( z )  with  respect  to  the 
origin  is  p^  for  each  i  -  1  to  k.  Consider  next  a  complex 
valued  function  h(t,0)  defined  for  0  <  t  <  1,  0  <  0  <  2re  as 
in  (2.18),  where  the  ♦i(0)'s  in  (2.18)  are  as  defined  in 
(2.19) . 


cl  -  gi(l)  -  exp(jei(0)),  <^(2*)  -  2*pt  ♦  ♦i(0)  (2.17a,b) 

h(t,0)  -  f(4l(0),02(0),...0k(0))  (2.18) 

♦  ^0)  -  exp  (jtei(9)  ♦  j  ( 1-t )  (  ♦i  (  0  )  ♦  pi0))  (2.19) 

for  each  i  ■  1  to  k. 


Subtituting  t-1  in  (2.19)  and  subsequently  making  use  of 

(2.18)  and  (2.19)  we  have  h(l,G)  -  f<91(z)  ,g2(z) , . .  .gk(z) ) 

with  z  -  exp(jd).  Similarly,  Ly  substituting  t  -  0  in  (2.19) 

and  subsequently  making  use  gf  ( 2 . 18 )^( 2 . 19 )  and  (2.17a)  we 

obtain  h(O,0)  -  f (c^z  ^CjZ  ^,...cJtz  k)  with  z  -  exp(j0). 

Also,  it  follows  from  ( 2 . 17 ) , ( 2 . 18 )  and  (2.19)  that  h(t,0)  ■ 

h(t,2n)  -  f(C]/c2, . .  .c^)  for  all  t  in  0  <  t  <  1.  Furthermore, 

since  the  <f^(0)'s  are  continuous  functions  of  6,  and  f(z)  is 

a  polynomial  and  hence  a  continuous  function  of  £  • 

(z^,z2,  •  •  .Zjj) ,  it  follows  in  view  of  (2.18)  and  (2.19)  that 

h(t,0)  is  a  continous  function  of  t  and  0.  The  function 

h(t,0),  when  viewed  as  a  function  of  0  only,  can  therefore  be 

thought  of  as  representing  a  family  of  closed  contours  with 

their  initial  and  terminal  points  fixed  at  f ( c^ ,c2 , • . .ck ) , 

which  are  continously  parametrized  by  the  variable  t  in  such 

a  way  that  we  obtain  Tq  when  t  -  0  and  when  t  -  1.  In 

addition,  since  it  follows  from  (2.19)  that  |+,(0)|  -  1  for 

each  i,  we  have  from  condition  (a)  of  the  present  theorem 

that  h(t,0)  #  0  for  all  0  <  t  <  1,  0  <  0  <  2it.  Therefore,  the 

function  h(t,0)  can  be  taken  to  represent  a  continuous 

deformation  of  the  contour  Tq  into  1*^  with  the  initial  and 

terminal  points  fixed  at  f (c^,c2, . . .c^)  such  that  at  any 

intermediate  stage  of  the  continous  deformation  process  the 

contour  may  never  pass  through  the  origin  (0,j0)  of  the 

complex  plane.  It,  therefore,  follows  that  the  number  of 

encirclements  of  the  origin  (0,j0)  by  Tn  and  r.  are  the 
1  u  1 

same . 

Next,  by  using  arguments  similar  to  the  one  used  above  we 
show  that  the  number  of  encirclements  of  the  origin  by  the 

*This  intuitive  notion  is  more  formally  expressed  by  saying 
that  if  Tq  and  are  Q— homotopy  of  each  other  then  they  are 
Q-homologous  (theorem  13.15  in  (18j).  Here,  Q  is  the  set  of 
all  complex  numbers  except  the  origin  (0,j0). 
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contour  is  zero.  Consider  the  family  of  circles  described 
by  the  complex  function  k(r,s)  of  real  variables  r,s  e  [0,1] 
defined  as  k( r,s)«rexp{2iijs) .  Also,  consider  the  function 


K( r,s) 


f ( 2 ( k ( r,s) ) ) .  Then  K(l,s)  describes 


K(0,s)«f (g(0) )  for  all  s.  Thus,  if  r  changes  from  1  to  0,  the 
contour  r  described  by  K(r,s)  for  s  e  [0,1]  contracts 
continuously  from  into  the  point  f(g(Q))*0.  Hence  the 
number  of  encirclements  of  the  origin  must  be  zero  for  r 
sufficiently  close  to  0,  and  it  must  therefore  be  equal  to 
zero  also  for  r-1,  because  due  to  condition  (b)  of  the 
present  theorem  r  never  goes  through  the  origin  for  r  e 
(0,1]. 2 

Thus  the  number  of  encirclement  of  the  origin  by  the  contour 
rn  is  zero.  Since  rn  is  the  image  of  the  unit  circle  |z|«l 


gapping 


defined 


function 


f(c^z  ,c2z  ♦•••CjJz  ),  it  follows  gy  using  the  principle  of 
argument  that  1,c2*  2,...ckz  *)  *  0  in  |z|  <  1.  Since 
Pj^  >  0,  and  from  (2.17a)  we  have  that  (c^J  ■  Ig^DI  -  1  this 
latter  statement  along  with  condition  (a),  due  to  Corollary 
2.7.2,  completes  the  proof  the  present  theorem. 

Before  we  can  undertake  the  proof  of  the  extended  version  of 
Rudin's  theorem  reported  in  [17],  we  need  the  following 
lemma . 


Lemma  2.8:  Let  f(z)  be  a  polynomial  in  :  ■  ( z^ , z2 , . . . z^ )  such 
that  f(z)  t  0  in  |z|  -  1.  Also,  let  the  indices  i  -  1  to  k  be 
the  union  of  disjoint  sets  of  indices  J^,  J2  and  J^.  If  f(z) 
satisfies:  (a)  f(z)  #  0  for  | |  $  1,  i  c  ;  zi  -  «i ,  i  c 
(J2  U  J3),  and  (b)  f(z)  t  0  for  all  z  with  zi  -  0i  for  i  c 
J1  zi  "  9*(*)  for  i  e  J2  and  |z|<l,  and  zi  -  ri  for  i  c  J3, 
where  the  goal's  are  functions  satisfying  the  same 

2This  step  does  not  indeed  follow  from  the  principle  of 
argument . 
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conditions  as  those  in  Rudin's  theorem,  and  the  ,  0^,y^'s 
are  some  complex  numbers  such  that  |  |  ■  |  I  **  1  and  |  0.J  < 
1,  then  f(z)  *  0  for  Iz-J  <  1,  i  e  (Jj^  U  J2);  i  e 

J3. 

Proof:  Since  f(z)  *  0  in  \z\  -  1,  invoking  Lemma  2.6  along 
with  condition  (a)  of  the  present  lemma  yields  that: 

f(z)j*0  for  |zi|<l,  ieJ^;  |*i|-l,  ie(J2  U  j3)  (2.20) 

Condition  (2.20)  implies,  in  particular,  that  f(jz)  *  0  for  z^ 
■  3^,  i  e  J^;  |z^|  »  1,  i  e  J2  and  Zj^  ■  y^,  i  e  J3,  which 

along  with  condition  (b),  due  to  Rudin's  theorem,  implies 
that  f(z)  *  0  for  z^^  -  i  e  J^,  |  z  i  |  <1,  i  e  J2  and  z^^  « 
r^,  i  e  J3.  This  latter  conclusion,  along  with  a 

particularization  of  (2.20)  via  the  choice  of  z^  -  ri >  i  e 
J3,  due  to  Lemma  2.3,  yields  the  desired  result. 

We  can  now  prove  the  generalized  version  of  Rudin's  theorem 
stated  as  follows. 

Theorem  2.11:  Let  the  set  of  indices  i  -  1  to  k  be  the  union 
of  n  disjoint  subsets  of  indices  Ij,  j  ■  l,2,...n.  Then  the 
polynomial  f(z)  in  z  -  (z^,z2, . .z^)  is  devoid  of  zeros  in  |:z| 
<  1  if  and  only  if  the  following  hold  true  simultaneously: 

(a)  f(z)  *  0  in  |z|  -  1 

(b)  for  each  y,  y  -  l,2,...n,  the  polynomials  obtained  by 

setting  zi  -  0^,  i  c  Ij,  j  ■  1  to  (y-l);  zi  -  g^fz),  i  z  1IJ 

and  zi  ■  r ^fJ,  i  e  1^,  j  -  (*/+l)  to  n  in  f(z)  are  devoid  of 

zeros  in  |z|  <  1,  where  the  g^(z)'s  are  functions  satisfying 

the  same  constraints  as  those  satisfied  by  the  g^(z)'s  in 

Rudin's  theorem,  and  the  0.  's  and  y.  's  are  some  complex 

iy  l  y  r 

numbers  such  that  |0.  |  <1  and  ly.  I  ■  1. 

iy  1  i y ' 


Proof:  Only  sufficiency  needs  to  be  proved.  We  claim  that  for 

all  m  in  1  <  m  <  n  tjje  polynomial  Bm  obtained  by  freezing 

at  z.  -  y.  ,  is  U  I.  satisfies  *  0  in  |z. |  <  1,  i  e 
mi  lm  i-m+1  3  mi 

U  1^.  The  proof  of  this  assertion  is  via  induction  on 

ilAce  B.  is  obtained  from  f(z)  by  setting  z.  -  y..,  i  s  U 

1  X  11  •  ^ 

Ij  and  | ri ^ I  ■  1#  it  follows  from  condition  (a)  that  B^  *  0 
for  | z ^ |  -  1,  i  e  I ^ .  Furthermore,  condition  (b)  with  u  -  1 
implies  that  B^  *  0  in  |z|  <  1  with  z^  -  g^(z),  i  c  1^.  The 
last  two  conditions,  in  view  of  Rudin's  theorem,  prove  that 
B^  *  0  in  |ZjJ  <  1,  i  e  1^.  Therefore,  our  assertion  is  true 
for  m  -  1.  We  next  assume  that  the  assertion  is  true  fgr  m. 


with  1  < 


<  n-1 .  i.e.,  f ( z )  t  0  for  |z.|  <  1,  i  e  U 
n  "  1 


and  z,  -  y.  ,  i  e  U  I..  Applying  Lemma  2.8  on  f(z)  4ith 
1  im  j-m+1  3 


XJ 


m 


I.  along  with  condition 


*  m+T+1 '  ifc 


■  U  I  j  #  2  *  9  ^ 3  m  ^  ^  j  dlon^  w 

(a)  and^cindition  (b)  of  the  preSeRt^lemma,  with 
directly  follow^  that  f(z)  *  0  in  |z^|  <1  for  i  e  U  Ij*  zi 
y ;  ;  i  c  U  I.,.  Therefore,  Bffl+1  ^  0  in  |zi|371l, 


j-m+2 


'j 


i ,m+l ' 

m+1 

i  c  U  Ij.  The  proof  of  our  assertion  via  the  induction  is 
thus  c^m^lete.  The  theorem  then  follows  by  noting  that  br  and 
f(z)  are  identical  in  all  the  variables  z^,  i  ■  1  to  k. 


Remark :  We  note  that  Rudin's  theorem  is  a  special  case  of 
Theorem  2.11  when  n-1.  Also  if  0.  -0.  for  all  u  then  we 

ifj  i 

obtain  the  extension  of  Rudin's  theorem  reported  in  [17].  in 
this  sense  Theorem  2.11  can  be  considered  as  a  slightly 
generalized  version  than  that  reported  in  [17]. 

As  a  corollary  to  the  above  theorem,  we  have  the  following 
result  reported  in  [20]. 

Corollary  2.11.1:  The  polynomial  f(z)  in  z  -  ( , z2 , . .  •  z^  )  is 
devoid  of  zeros  in  |z|  <  1  if  and  only  if  the  following  hold 

true  simultaneously. 


(b)  for  each  u ,  fj  -  l,2,...k,  the  polynomials  obtained  by 
setting  z i  -  8^,  i  »  1  to  (n-l),  z^  -  z,  i  -  n,  and  zi  -  y^  , 
i  -  (/j+1)  to  k  in  f{z)  are  devoid  of  zeros  in  |ft|  <  1,  where 
0^ ' s  and  r^'s  are  some  constants  such  that  |  fT  |  <  1  and 


Proof :  The  proof  of  the  above  corollary  clearly  follows  from 
Theorem  2.11  by  choosing  n  -  k,  Ij  “  { j }  for  j  -  l,2,...k, 
0^-0^  and  g^(z)  ■  2  for  each  i  and  fj . 

Remark :  It  is  possible  to  formulate  a  more  direct  proof  of 
Corollary  2.11.1  without  making  use  of  the  Rudin's  theorem 
via  the  use  of  Lemma  2.3,  Lemma  2.6  and  a  strategy  similar  to 
the  one  adopted  for  the  proof  of  Theorem  2.11.  However,  the 
details  of  such  a  proof  is  omitted  from  the  present 
discussion  for  the  sake  of  brevity. 

Also  note  that  if  0^  «  0  for  all  i,  and  r ^  »  1  for  all  i  and 
//  then  Corollary  2.11.1  coincides  with  a  result  stated  and 
proved  in  Theorem  2.4  in  [7J. 
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2.6.  Conclusion 


A  correct  formulation  of  the  continuity  property  of  zeros  of 
a  multivariable  polynomial  as  a  function  of  its  coefficients 
has  been  discussed.  This  property  has  been  used  to  derive  all 
the  conditions  for  a  polynomial  f(z)  in  z  -  ( z^ , Z2 , . • - , )  to 
be  devoid  of  zeros  in  jz|  <  1  previously  known  in  the 
literature.  The  proofs  presented  here  are  believed  to  be 
simple,  rigorous  and  more  intuitive  than  those  published 
earlier.  The  present  report  deals  with  results  arising  from 
studies  of  bounded-input-bounded-output  property  of 
multidimensional  discrete  time  systems  only.  It  is  well  known 
[1,10]  that  unlike  in ' one-dimension,  the  obvious  continuous 


domain  analogs  of  certain  discrete  time  domain  results  do  not 
hold  true  unless  proper  modifications  necessitated  by  the 
non-compactness  of  the  right  half  polydomain  are  duely  made. 
The  proofs  of  these  results  in  the  multidimensional  (k>2) 
context  based  on  the  continuity  property  of  the  zeros  of  a 
polynomial  as  a  function  of  its  coefficients  can  also  be 
worked  out  along  lines  similar  to  those  discussed  here.  This 
discussion  will  be  the  content  of  a  forthcomming  report  [21]. 
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CHAPTER  3 


NEW  RESULTS  ON  STABLE  MULTIDIMENSIONAL  POLYNOMIALS 

3.1.  Introduction 

Recent  studies  on  the  scattering  parameter  description  of 
passive  multidimensional  systems  have  given  rise  to  a  new 
class  of  multidimensional  Hurvitz  polynomials,  called 
scattering  Hurvitz  polynomials,  by  allowing  zeros  of 
restricted  nature  to  occur  on  the  boundary  of  the  domain 
under  consideration  (1].  Originally,  two  different 
definitions  of  these  polynomials  were  introduced,  the 
equivalence  of  which  has  been  demonstrated  via  an  extension 
of  the  maximum  modulus  theorem  for  analytic  functions  of 
several  complex  variables  [2],  and  some  of  their  properties 
have  been  discussed  in  (3].  Nontrivial  properties  of 
scattering  Hurvitz  polynomials  occuring  in  discrete  time 
domain  applications  are  considered  in  [4]  and  their  testing 
procedures  are  elaborated  in  the  two-dimensional  context  in 
[5].  An  alternative  approach  to  deriving  the  basic 
properties  of  scattering  Hurvitz  polynomials  has  also  been 
offered  (6).  The  present  report,  however,  will  be  organized 
such  that  all  proofs  are  complete  without  requiring  any 
knowledge  of  ( 6  ]  . 

The  aforementioned  investigation  is  carried  further  in  the 
present  report  by  classifying  a  wider  variety  of 
multidimensional  polynomials  occurring  in  transfer  function 
descriptions  of  passive  systems.  Whereas  the  scattering 
Hurvitz  polynomials  occur  as  the  denominators  of  bounded 
functions,  reactance  Hurvitz  and  immittance  Hurvitz 
polynomials,  as  defined  later,  are  characterized  as  the 
denominators  (and  hence  the  numerators)  of  reactance 
functions  and  immittance  functions  respectively.  Related 
other  results  on  properties  of  multidimensional  polynomials 


and  rational  functions  are  discussed  in  this  context.  Some 
related  investigations  in  [7]  have  been  brought  to  the 
attention  of  the  authors. 

Notations,  terminologies  and  definitions  and  some  general 
properties  of  multivariable  polynomials  to  be  used  in  the 
rest  of  the  report  are  introduced  in  Section  3.2.  Properties 
of  widest-sense  Hurwitz  polynomials  and  self-paraconjugate 
Hurwitz  polynomials  are  discussed  in  Section  3.3.  Properties 
of  scattering  Hurwitz  polynomials  previously  unpublished  in 
the  literature  are  discussed  in  Section  3.4.  Section  3.5 
deals  with  properties  of  positive  functions,  whereas  the 
reactance  Hurwitz  and  immittance  Hurwitz  polynomials  are 
defined  and  their  properties  studied  in  Sections  3.6  and  3.7 
respectively.  In  Section  3.8  a  few  results  potentially 
useful  for  testing  the  positivity  property  of 
multidimensional  rational  functions  are  derived.  A 
nontrivial  result  concerning  the  property  of  nonnegativity  of 
the  real  part  of  a  rational  function,  analytic  in  the  right 
half-polyplane,  in  terms  of  the  behavior  of  the  function  on 
the  distinguished  boundary  of  the  domain  of  holomorphy  is  to 
be  noted  in  this  connection.  Finally  results  obtained  are 
summarized,  and  conclusions  are  drawn  in  Section  3.9. 


3.2.  Notation,  terminology,  definitions  and  general 
properties  of  multivariable  polynomials. 


The  following  is  a  partial  glossary  of  notations  and 
terminologies  to  be  used  in  succeeding  discussions. 


A  polynomial  g  in  k  variables  Pi*P2'***Pk  will  be  denoted 
simply  by  g  or  g(p2,p2*  •  •  .p^)  or  by-  9(E)*  where  g  denotes  the 
k-tuple  of  variables  g  -  (p1 rp2> . . .p^) .  We  will  also  write  g 
-  g(g)  or  g  ■  9(p^*P2* • • -p^)  to  indicate  that  g  is  a 
polynomial  in  k-variables.  We  take  for  granted  that  a 
polynomial  g  in  k  variables  may  be  independent  of  one  or  more 
of  the  variables  P^*P2* . * .p^.  The  k-variable  polynomial  g 
will  be  said  to  involve  a  variable  p4  if  the  indeterminate  p. 
actually  exists  in 
the  polynomial  g. 
variable  polynomials  obtained  by  assigning  arbitrary  values 
to  p^  cannot  all  be  identical.  A  polynomial  g  -  g(g)  will  be 
called  nonconstant  if  g  involves  at  least  one  of  the 

The  set  of  integers  1 


at  least  one  of  the  monomials  composing 
If  g  ■  g(g)  involves  pi  then  the  (k-1) 


I  m  { 1 , 2 ,  .  .  .  k }  .  By  i1  ,  i, ,  .  .  .if.  we 


to  k  will  be 
k 


variables  p1 ,p2 , . . .pk . 

de s i gna ted  by  I  i •  e . ,  *  —  ( • .a; •  bj  ^  l ^  2 

designate  a  permutation  of  the  integers  l,2,...k.  A 

nonconstant  factor  of  the  form  dv  where  v  is  an  integer  2 

is  said  to  be  a  multiple  factor,  and  a  nonconstant  factor 
that  is  irreducible  is  called  an  irreducible  factor. 


If  g  ■  g(g)  is  written  as  a  polynomial  in  p^,  id,  as 
“i  v 

9  -  £  *  (3.1) 

v-0 

where  the  coefficients  A^  are  polynomials  in  the  remaining 
variables,  with  A.  M  0,  then  n,  is  called  the  partial  decree 
of  g  in  the  variable  p^  and  is  to  be  denoted  by  deg^.  Two 
polynomials  will  be  said  to  be  relatively  prime  if  they  do 
not  have  a  nonconstant  common  factor.  The  terms  factor 
coprime  and  proper  factor  are  also  to  be  used  respectively 


for  the  terms  relatively  prime  and  nonconstant  factor. 


The  asterisk  *,  when  used  as  a  superscript  along  with  any 
scalar  expression  (or  a  constant),  will  indicate  complex 
conjugation.  The  paraconjugate  of  a  polynomial  (also  of  a 
rational  function)  g  is  defined  as:  g*  »  g*(p)  - 
g*(-p*,-P2, . • .-p£) •  A  polynomial  g  is  said  to  be  self- 
paracon jugate  if  g*-  Cg,  where  C  is  a  constant  (necessarily 
unimodular ,  i.e.,  |C|-1).  g  is  said  to  be  paraeven  or 


paraodd  if  c  ■  1  or  C  ■  -1  respectively.  Sometimes  it  will 
be  appropriate  to  write  the  k-tuples  £  as:  £  ■  (Pj/E')* 
where  p'  indicates  the  (k-l)-tuple  £'  -  (p2#p3, . . .pk) . 
Correspondingly  we  will  also  write:  g  »  g  (p^,£').  A 
polynomial  f  in  the  variables  p2,p3,...pk  will  be  expressed 


as  f  -  f(p2»P3/ 
specific  iel, 


.p^)  or  equivalently  as  f  -  f(p').  For  any 
we  will  use  the  phrase  "( k-1 ) -variable 


polynomial  obtained  by  freezing  the  variable  p^  in  g(p)"  to 
mean  the  (k-1)  variable  polynomial  obtained  from  g(£)  by 


assigning  a  fixed  value  to  p^  A  second  subscript  to  the 
variables  pi#  1  <  i.  <,  k  such  as  piQ,  will  usually  mean  a 
fixed  value  of  pi.  Correspondingly  the  notations  pQ  - 

( P^O 'P20 '  *  *  *^k0 ^  an<*  £0'  "  ^P20'£30' *  *  *£k0^  etc*  wiH  be 
used.  The  notation  Rep  >  0  (or  Rep'  >  0)  will  be  taken  to 

mean  Rep^  >  0  for  all  iel  (or,  ie{2, 3, . . . ,k) 
correspondingly),  etc.  (all  p^  involved  being  obviously 
assumed  finite).  The  symbol  u  will  be  used  exclusively  for 
designating  real  numbers.  Thus,  notations  such  as  p^  -  jw. 
and  piQ  -  j«iQ  imply  that  Rep.^  -  0,  Repi0  ■  0  respectively. 
Similarly  u  is  used  for  the  k-tuple  of  real  numbers 
( w^,w2, . • ) ,  i.e.,  £  ■  jw  implies  Rep  -  0  etc.  Wherever 
appropriate,  definitions  and  notations  discussed  so  far  also 
apply  if  g  is  more  generally  a  rational  function  in  p.  The 
notation  I  |  .  I  I  denotes  any  suitable  norm,  say,  the  Euclidean 
norm,  k 


Let  P  be  a  set  of  k-tuples  £  -  (p^pj,  — pk) ,  where  all  p^ 
belong  to  the  same  number  field  K  (hereafter  always  the  field 
of  real  numbers  or  the  field  of  complex  numbers).  We  will 
say  a  certain  property  holds  for  almost  all  values  if  a 
variable  may  be  equal  to  any  element  of  the  field  except, 
possibly,  finitely  many  of  them.  The  set  of  all  values  that 
the  variable  may  then  take  is  said  to  be  almost  complete. 
The  symbol  S  will  be  reserved  to  denote  the  set  P  when  the 
variables  are  restricted  to  be  real. 


Definition  3.2.1a:  We  say  that  P  is  a  sequentially  almost 

complete  set  of  order  m  >  1,  with  m  <  k,  if  there  exists  a 
permutation  , i2 , . . . , i^  of  the  integers  1,2, ...,k  such  that 
all  £  c  P  can  be  generated  in  the  following  way:  There 
exists  an  almost  complete  set  K.ck  such  that  any  p.  c  K.  may 
be  chosen.  For  any  choice  thus  made,  assuming  m  ?  2,  there 
exists  an  almost  complete  set  K2ck  (possibly  depending  on  the 
particular  p^  c  selected)  such  that  any  p^  c  may  be 
chosen.  Again  for  any  choice  thus  made,  assuming  m  >  3, 
there  exists  an  almost  complete  set  K^ck  (possibly  depending 
on  the  particular  p^  and  p^  selected)  such  that  any  p^  c 
Kj  may  be  chosen,  e&c.  If  2  »  k  this  process  is  continued 

until  we  have  reached  p.  . 

1  k 

If  m  <  k,  once  we  have  reached  i  there  exists  at  least  one 


(k-m)-tuple  (p. 


)  (possibly  depending  on  the 


particular  p^  toa$i  selected)  that  may  be  chosen.  Finally, 
we  may  extend  the  aBove  definition  to  the  situation  m  -  0  by 
saying  that  in  this  case  the  set  P  is  not  empty. 


Definition  3.2.1b:  P  is  sequentially  infinite  of  order  m, 

1  i  m  £  k,  if  it  can  be  generated  as  in  Definition  3.2.1a 
except  for  replacing  everywhere  the  term  "almost  complete 
set"  by  the  term  "infinite  set". 


Note : 


In  Definitions  3.2.1a  and  3.2.1b,  the  permutation 


( i  1 ,  i 2 , . .  .i^)  will  be  called  the  ordering  of  P.  If  i^  -  a*, 
H  -1  to  k,  then  P  is  said  to  be  ordered  naturally. 


Definition  3.2.1c:  P  is 

(n^ ,n2, . . .nk )  with  order  a, 


sequentially  exceeding  n 
where  a  <  k,  if  it  can  be 


generated  as  in  Definition  3.2.1a  except  that  the  terms  an 
almost  complete  set  K^ck  etc.  are  replaced  by  "a  set  K^ck 


comprising  at  least  n^  ♦  1  elements”  etc.,  the  n^,  i  ■  1  to 


k,  being  finite,  nonnegative  integers. 


An  obvious  interrelationship  between  the  sets  defined  above 
is  also  summarized  in  the  following  theorem. 


Theorem  3.2.1:  A  sequentially  almost  complete  set  of  order  m 
is  also  sequentially  infinite  of  order  m,  and  a  set  of  the 
latter  type  is  sequentially  exceeding  n  »  ( n^ ,n2 , . • .nk )  with 
order  a,  and  this  for  any  choice  of  the  n. 


Definition  3. 2. Id:  P  is  almost  complete  of  order  m,  infinite 
of  order  m,  or  exceeding  n  ■  ( n1 , n^ , . . . n^ )  with  order  m  if 
P  -  Px  x  P2  x...Pk,  the  sets  Pj  ,  i  ■  1  to  k,  being  non-empty 
and  such  that  at  least  m  of  them  are  almost  complete, 
infinite,  or  contain  a  number  of  elements  larger  than  the 
corresponding  integer  n^,  respectively,  where  m  <  k. 

Some  properties  of  the  type  of  sets  just  defined,  as  they 
relate  to  the  sets  defined  earlier  in  Definitions  3.2.1a, 
3.2.1b  and  3.2.1c  are  mentioned  in  the  following. 


Theorem  3.2.2:  A  set,  P,of  k-tuples,  £,  that  is 


1.  almost 'complete  of  order  m  is  sequentially 
almost  complete  of  order  m. 


2. 


infinite  of  order  m  is  sequentially  infinite  of 
order  m, 


3. 


exceeding  n-(n1>n2, _ n^)  with  order  m  is  sequentially 

exceeding  n  »  (n^,n2» . • .n^)  with  order  m. 

Proofs  of  these  results  clearly  follow  from  a  close 
examination  of  the  definitions  of  the  sets  involved. 

For  the  purpose  of  the  present  report  we  will  also  adopt  the 
following  definitions.  More  definitions  *nd  terminologies 
will  be  introduced  as  they  occur  in  the  main  body  of  the 
text. 


Definition  3.2.2:  A  polynomial  g  -  9(E)  is  widest-sense 

Hurwitz  if  g(£)  *  0  foe  Re  p  >  0. 

Definition  3.2.3;  A  polynomial  g  ■  g(p)  is  strict-sense 

Hurwitz  if  g(£)  *  0  for  Re  g  >  0. 

Definition  3.2.4;  A  polynomial  g  »  g(£)  is  scattering 

Hurwitz  if  the  following  conditions  simultaneously  hold  true: 

(i)  g(£)  *  0  for  Re  g  >  0,  i.e.,  g  is  widest-sense 

Hurwitz 

(3.2) 

(ii)  g  and  g#  are  relatively  prime  polynomials. 

(3.3) 

Definition  3.2.5:  A  polynomial  g  -  9(£)  is  a 

self-paraconjugate  Hurwitz  polynomial  if  it  is  a  widest-sense 
Hurwitz  polynomial  and  is  self-paraconjugate. 

Definition  3.2.6:  A  polynomial  g  -  g(£)  is  reactance  Hurwitz 
if  it  can  be  written  as  a  constant  (possibly  complex) 
multiplied  by  the  paraeven  or  paraodd  part  of  a  scattering 
Hurwitz  polynomial. 


Definition  3.2.7:  An  immittance  Hurwitz  polynomial  is  the 

product  of  a  scattering  Hurwitz  and  a  reactance  Hurwitz 
polynomial . 

Definition  3.2.8:  A  function  F(£)  is  called  a  positive 

function  if  Re  F(g)  >  0  everywhere  in  Re  £  >  £,  where  F  is 
holomorphic. 

The  positive  function  F  ■  jC,  where  C  is  a  real  constant  is 
said  to  be  trivial .  All  other  positive  functions  are  non¬ 
trivial  . 

Definition  3.2.9:  A  paraodd  rational  positive  function  is 

called  a  reactance  function. 

Note  that  the  Definitions  3.2.8  and  3.2.9  do  not  assume  the 
function  under  consideration  to  be  a  real  function. 

Theorem  3.2.3:  If  g(£)  is  a  polynomial  in  k-variables  such 

that  the  set  of  zeros  of  g  comprises  a  sequentially  infinite 
set  of  order  k,  then  g  is  identically  equal  to  zero.  More 
generally,  let  ni  -  deg^  ,  i  ■  1  to  k.  If  g  -  0  for  all 
£  e  P,  where  P  is  sequentially  exceeding  n  -(n^nj,  . .  .n^) 
with  order  k,  then  g  is  identically  equal  to  zero. 

Proof :  Assume  that  a  set  P  of  the  type  mentioned  exists,  but 
that  g  M  0.  We  may  assume  P  to  be  ordered  naturally. 
Consider  g  as  a  polynomial  in  £' ,  the  coefficients  A^tp^)  of 
which  are  polynomials  in  p^  only;  the  Av  are  of  degree  <  n1 . 
Obviously,  there  exists  a  v'  such  that  A^fp^)  jt  0.  Hence, 
among  the  values  of  p^  to  be  considered  for  forming  P,  there 
must  exist  at  least  one,  say  p1Q,  with  Av,  (p1Q)  *  0.  The 
polynomial  g^,  defined  by  g^(£')  ■  g(p^Q,£')  is  t^ien  not 

identically  zero. 

Next,  we  proceed  with  g^  as  before  with  g  i.e.,  consider  g^ 
as  a  polynomial  in  ( p^ , p4 , • • • p^ )  the  coefficients  B v ( p 2 )  of 


which  are  polynomials  in  p 2  etc....  Finally,  we  arrive  at 
g(£0)  i*  0  where  £Q  e  P.  This,  however,  is  a  contradiction. 

Lemma  3.2.4:  For  any  non-constant  polynomial  g,  there  exists 
a  sequentially  almost  complete  set,  P,  of  order  (k-1)  such 
that  g{£)  -  0  for  all  £  e  P. 

Proof :  We  may  assume  that  g  involves  p^.  Freeze  p1  at  p1Q 
and  consider  the  polynomial  g^  defined  by  g1 ( £ ' )  -  g(P10#£'). 
Since  the  coefficients  of  g^^  are  polynomials  in  p1Q,  there 
are  at  most  finitely  many  values  of  p^g  for  which  g^(£')  is 
independent  of  p^.  For  any  other  choice  of  p1Q  we  may  apply 
the  same  argument  to  g^  as  formerly  to  g,  with  p2  taking  the 
role  of  pj^  etc.  Finally,  we  find  a  polynomial  g^^P^)  that 
still  involves  p^  and  thus  has  at  least  one  zero. 

Theorem  3.2.5:  If  f  and  g  are  polynomials  in  k  variables 
then  f  and  g  have  a  proper  common  factor  if  and  only  if  the 
set  of  zeros  that  are  common  to  f  and  g  is  sequentially 
infinite  of  order  (k-1). 

Proof :  Necessity  is  obvious  in  view  of  Lemma  3.2.4.  To 
prove  sufficiency,  let  P  be  the  set  mentioned  in  the 
statement  of  the  theorem  and  assume  that  is  the  last 
variable  selected  in  forming  P,  i.e.,  in  the  terminology  of 
Definition  3.2.1a,  that  we  have  i^  -  1.  Let  P'  be  the  set  of 
(k-l)-tuples  £'  involved  in  forming  P.  There  exist 
polynomials  u,v  and  w  such  that  uf  +  vg  -  w,  where  w-w(£')  is 
a  polynomial  independent  of  p^,  deg^  u  <  deg^  g  and  deg^  v  < 
deg^  f;  a  polynomial  w  thus  defined  has  the  property  that  w  > 
0  holds  if  and  only  if  f  and  g  have  a  common  proper  factor 
involving  p^  (8J.  Furthermore,  since  P'  is  sequentially 
infinite  of  order  (k-1)  it  is  clear  from  Theorem  3.2.3  that 
w(£' )  is  actually  identically  zero,  i.e.,  that  f  and  g  have 
indeed  a  common  factor. 


Lemma  3.2.6;  Let  g  be  a  polynomial  in  £.  Let  us  select  one 
of  the  p^,  say  p^,,  and  let  us  freeze  p^,  at  say,  p^,g. 
There  exist  at  most  finitely  many  choices  of  P^,q  such  that 
any  partial  degree  of  g  is  lowered. 

Proof :  We  may  assume  i'  -  1.  Write  g  as  a  polynomial  in  £' 
whose  coefficients  are  polynomials  in  p^.  Clearly,  there  are 
at  most  finitely  many  values  of  p^g  such  that  the  leading  one 
of  these  coefficients  becomes  zero. 


Theorem  3.2.7:  Let  f  and  g  be  two  relatively  prime 

polynomials.  For  any  m  such  that  1  <  m  <  k  let  us  freeze  m 

of  the  variables  p. ,  say  for  i  ■  i,  to  i  ,  at  corresponding 

i  i  n& 

values  PiQ«  Let  ^  and  <3±  be  the  resulting  polynomials  in 
the  remaining  variables.  Then  there  exists  a  sequentially 


almost  complete  set  Pm,  of  m-tuples  of  order  m  such  that  for 


(P|  /Pj  »•••?< 
reli?ivei^  prime?0 


c  Pm,  the  polynomials  f^  and  g^  are 
Furthermore,  for  any  ordering  chosen 


still 

there 


exists  a  set  Pm  with  the  property  given. 


Proof :  It  is  enough  to  prove  the  theorem  for  m  -  1  and  to 
assume- that  the  variable  to  be  frozen  is  p^.  Define  u,  v  and 
w  as  for  the  proof  of  Theorem  3.2.5,  and  let  u^,  v^  and  w^^  be 
the  polynomials  resulting  from  the  former  set  of  polynomials 
by  freezing  p^  at  p^g.  Clearly,  w  is  not  identically  zero, 
and  there  are  at  most  finitely  many  values  of  p^g  for  which 
at  least  one  of  the  relations  deg1f1  <  deg^,  deg1g1  <  degxg 
(cf.  Lemma  3.2.6)  or  w^pj, — *  0  could  hold.  For  all 
other  choices  of  p^g  the  conditions  for  ensuring  that  and 
gj^  have  no  common  factor  involving  p^^  are  fulfilled.  In  a 
similar  way  one  can  show  that  there  are  at  most  finitely  many 
p^g  for  which  and  g^  have  a  common  factor  involving  any  of 
the  variables  p2  to  p^^. 


Theorem  3.2.8:  For  any  polynomial  g  that  is  not  identically 


p 


zero  there  exists  a  set,  P,  of  k-tuples  p,  that  is 
sequentially  almost  complete  of  order  k  and  such  that  g(p)/0 
for  peP.  Furthermore,  a  set  P  with  the  given  property  exists 
for  any  choice  of  the  ordering  of  this  set. 

Proof:  It  is  sufficient  to  assume  natural  ordering.  Let  g^ 
be  the  polynomial  in  p'  obtained  by  freezing  p^  at  i.e., 
gx(p' )  -  g(p^Q,p' ) .  The  coefficients  of  g^  are  polynomials 
in  p^Q .  Hence,  g1(p')j<0  for  almost  all  choices  for  p^.  For 
any  of  these,  we  may  apply  the  same  argument  to  g^  as  before 
to  g  except  that  P2  now  plays  the  same  role  as  p^  had  etc. 
Finally,  we  arrive  at  a  polynomial  g^^  that  depends  on  pk 
alone  and  is  not  identically  zero,  which  is  thus  different 
from  zero  for  almost  all  p^. 


Theorem  3.2.9:  If  a  polynomial  g  is  devoid  of  zeros  in  the 
region  Reg  >  0  as  well  as  in  the  region  Reg  <  0  then  g  is  a 
nonzero  constant. 

Proof :  If  g  is  not  a  constant  let  us  assume  that  it 
involves,  say,  the  variable  p^ .  Let  A  be  the  leading 
coefficient  of  g  when  writing  it  as  a  polynomial  in  p^ 
clearly,  A  is  a  polynomial  in  g'  and  is  not  identically  zero. 
According  to  Theorem  3.2.8  there  exists  a  sequentially  almost 
complete  set,  P',  of  (k-l)-tuples  such  that  A(g')  ?  0  for 
g'cP'.  Hence,  we  may  consider  a  fixed  value  g£  -  jw£  such 
that  gg(p^)  ■  g(p^»j<Sg)  still  involves  p^.  Then  there  exists 
a  P^  “  P10  such  that  g(p^g,w£)  -  0.  This,  however,  leads  to 
a  contradiction,  because  both  Rep^g  >  0  and  Rep10  -  0  are 
impossible  due  to  our  hypothesis. 


Theorem  3.2.10:  Let  g  be  a  polynomial  in  g  and  assume  that 
there  exists  an  i'e  I  and  a  fixed  value  Pg  such  that  g(g)»0 
if  Pi*“Po  and  the  remaining  take  any  arbitrary  value 
belonging  to  a  sequentially  infinite  set  of  order  k-1.  Then 
g  contains  (p-pn)  as  factor. 


Proof : 


We  may  assume  i'«l.  Write  g  as  a  polynomial  in  £' 
whose  coefficients  are  polynomials  in  p^  .  In  view  of  Theorem 
3.2.3,  with  k  replaced  by  k-1,  all  these  coefficients  are 
zero  for  pi-pg. 

Theorem  3.2.11:  If  g(£)  is  a  polynomial  and  g(£g)*.  0  for 

some  £q  then  there  exists  an  n  >  0  such  that  g(£)  *  0  for  all 
£  in  the  neighbourhood  I|£-£qII  <  h*  More  generally,  for  any 
c  <0  there  exists  an  n  >  0  such  that  | g(£) I > I g( £g ) I  -  e  for 

l£-E0l<n* 

Proof :  Since  g  is  a  polynomial  and  hence  a  continuous 

function  of  £,  for  any  given  e  >  0  there  exists  an  n  >  0  such 
that  for  all  £  satisfying  II£-£qII  <  h  we  have  that 
|g(£)-g(£0)|  <  e,  i.e.,  |g(£) I > |g(£0) I-e,  thus  I g ( p ) I > 0  if  we 
choose  e< |g(£fl) | . 

Theorem  3.2.12:  Let  g  be  a  polynomial  in  £  having  a  zero  at 
£g.  Let  U  be  any  neighborhood  of  £Q .  There  exists  a 
sequentially  infinite  set  Pcu  of  k-tuples  £  that  is  of  order 
k-1  and  such  that  g(£)-0  for  £EP. 

Proof ;  Since  any  polynomial  can  be  decomposed  into  a  product 
of  irreducible  factors,  it  is  sufficient  to  assume  that  g  is 
irreducible.  On  the  other  hand,  the  result  is  true  for  k«l. 
Assume  thus  that  it  holds  for  k-1,  we  will  show  that  it 

remains  valid  for  k.  For  this,  consider  g  as  a  polynomial  g^ 

in  £'  whose  coefficients  are  polynomials  in  p^  and  write 
£g*( Pig ,£q ' ) .  In  view  of  Theorem  3.2.10,  g^  cannot  be 
identically  zero  for  P]_“Pio  since  otherwise  g  would  not  be 
irreducible.  Hence,  we  may  move  p^  from  p^  to  a  position 
p^  close  to  P]_q>  apply  to  g^  the  continuity  property  of  the 
zeros  of  a  polynomial  in  several  variables  [2],  and  conclude 
that  there  exists  a  £^'  such  that  g^(£,')-0  for  P]_“Pj_^- 

assumption,  there  exists  thus  a  sequentially  infinite  set  P' 


of  (k-l)-tuples  £'  that  is  of  order  k-2  and  such  that 
g^E'J-O,  thus  gtPjyE'J-O,  for  £'eP'  and  Pi"Pn* 
Furthermore,  for  any  p^  sufficiently  close  to  p^Q  any 
corresponding  value  £'  then  resulting  in  the  way  just 
explained  is  such  that  (p^,£')eU.  Hence,  a  set  P  as 
required  does  indeed  exist. 

This  result  together  with  The.rem  3.2.5  yields  immediately 
the  following  corollary: 

Corollary  3.2.12.1:  Let  g  and  h  be  two  relatively  prime 

polynomials  in  £  and  let  £q  be  a  common  zero  of  g  and  h. 
Then  in  any  neighborhood  of  £q  there  exist  points  £  for 
which,  say,  g(£)-0  and  h ( £ ) #0 . 

The  following  theorem  is  known  in  a  more  general  form  from 
the  theory  of  functions  in  several  complex  variables  [10], 
but  is  included  for  the  sake  of  completeness  and  in  order  to 
point  out  a  simple  proof  based  on  the  above  results. 

Theorem  3.2.13:  Consider  the  function  F-h/g  where  g  and  h 

are  relatively  prime  polynomials  in  £.  If  F  is  known  to  be 
bounded  for  all  those  £  in  a  domain  D  where  g(£)*0,  then 
g(£)*0  for  all  £eO. 

Proof :  without  loss  of  generality,  we  may  assume  I F ( £ ) ! < 1 

for  £cD.  Assume  that  there  exists  a  point  £q£D  such  that 
In  v*ew  Corollary  3.2.12  it  is  sufficient  to 
assume  h(£Q)#0.  Due  to  continuity  of  the  polynomial  g  we  may 
state  that  for  any  c>0  there  exists  an  n1>0  such  that 
lg(E^)|<e  for  all  £^  satisfying  I  t Ei"Eo I  I < ni •  Similarly,  due 
to  Theorem  3.2.11,  for  any  c>0  there  exists  an  02  such  that 
| h(£1 ) | < | h( £0 ) | -c  if  I  I Ei”Eo I  I < n2 •  Clearly,  we  may  choose 
e<|h(£g)|/2,  in  which  case  the  expression  for  h(g^)  becomes 
I  h  ( £^  )  |  >  e .  Choose  then  n>0  such  that  n<  and  n<n-,.  Due  to 

Theorem  3.2.8,  we  can  choose  £^  such  that  |  ig^-ggi  !  <  n  and 
gtE^)**),  i.e.,  by  assumption,  that  |h(£1)/g(g1)|<l.  This  is 


in  contradiction  with  |g<£|)|<c  and  |h(£^)|>c. 


3.3.  Properties  of  widest  Sense  Hurwitz  Polynomials  and  Self- 
paraconjugatc  Hurwitz  Polynomials 

The  proof  of  the  following  Theorem  is  rather  trivial: 

Theorem  3.3.1:  Let  g  be  a  widest-sense  Hurwitz  polynomial  in 
£.  The  following  holds: 

(i)  Preeze  one  of  the  p^ ,  say  p^,,  at  Pi'O  with  R*  Pi'0>0* 
The  resulting  polynomial  in  the  remaining  k-1  variables  is 
also  widest-sense  Hurwitz. 

(ii)  Factors  and  products  of  widest-sense  Hurwitz  polynomials 
are  also  widest-sense  Hqrwitz. 

Lemma  3.3.2:  Let  g  be  a  widest-sense  Hurwitz  polynomial.  Let 
be  the  polynomial  in  k-1  variables  obtained  by  freezing 
any  one  of  the  variables  pi#  say  pi#,  at  p^-jw^g.  Then* 

(i)  there  exists  an  almost  complete  set,  Q ' ,  of  real  numbers, 

such  that  gL  is  widest-sense  Hurwitz  if  w^gCQ'  and  that  g^ 

is  identically  zero  if 

(ii)  There  exists  an  almost  complete  real  set  Q"cq'  such  that 

for  far  ,qCQ",  g^  is  widest-sense  Hurwitz  and  has,  in  its  k-1 

variables,  the  same  partial  degrees  as  g. 

Proof :  We  may  assume  i'-l.  Write  g  as  a  polynomial  in  g' 

whose  coefficients  are  polynomials  in  p^. 

(i)  There  are  at  most  finitely  many  values  of  p^ ,  thus  a 
fortiori  finitely  many  such  values  with  Re  p^-0,  for  which 
all  these  coefficients  can  become  zero.  Choosing 

different  from  any  of  these  values,  moving  slightly  from 

?ur  ,  into  Pe  p  ^  >  0  ,  and  applying  the  continuity  property  cf 
the  zeros  of  a  polynomial,  we  see  that  a  zero  of  g^(g')  in 
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Re  £'>0  would  imply  a  zero  of  g  in  Re  £>0.  The  proof  of  (ii) 
follows  from  (1)  and  Lemma  3.2.6. 

The  following  result  follows  by  applying  Lemma  3.3.2  m  times: 

Theorem  3.3.3:  Let  g  be  a  widest-sense  Hurwitz  polynomial  in 
£,  and  let  £a  denote  an  m-tuple  obtained  by  selecting  any  m<k 
of  the  variables  p^^  to  pk .  Let  us  freeze  £a  at  2mo"j“mo  and 
let  9^  be  the  resulting  polynomial  in  the  k-m  remaining  p^ . 
Then 

(i)  9^  is  either  widest-sense  Hurwitz  or  identically  zero; 

(ii)  more  precisely,  for  any  ordering  selected,  there  exists 
a  sequentially  almost  complete  set,  9'a  ,  of  real  m-tuples 
such  that  9^  is  widest-sense  Hurwitz  if  and  that  it 
is  identically  zero  lf  s.o*°V 

(iii)  there  exists  a  sequentially  almost  complete  set, 
fl"mcfl'a  *  of  feel  m-tuples  such  that,  for  “aoc3"m  '  91  is 
widest-sense  Hurwitz  and  has,  in  its  k-m  variables,  the  same 
partial  degrees  as  g. 

Theorem  3.3.4:  If  9  is  a  proper  sel f-paracon jugate  Hurwitz 
polynomial  then  there  exists  a  sequentially  almost  complete 
set,  Q,  of  order  (k-1),  composed  of  real  k-tuples  w  such  that 
9  ( jw)  ■  0  for  any  w  c  9. 

Proof :  We  may  assume  that  g  depends  on  p^.  Consider  the  one- 
variable  polynomial  g^p^)  ■  gtp^.jwg),  obtained  by  freezing 
£'  at  £'  -  where  is  a  real  (k-l)-tuple.  Due  to 
Theorem  3.3.3,  there  exists  a  sequentially  almost  complete 
set,  9',  of  order  (k-1)  such  that,  for  any  c  9',  g^  is  of 
degree  >_  1  and  is  widest-sense  Hurwitz.  Hence,  there  exists 
at  least  one  P^g*  necessarily  witn  Re  £  0,  such  that 
gl(p10}  "  °*  However,  since  gfc  •  Cg,  p^  -  -  p^j  is  also  a 
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zero  of  g1#  where  Re  p£0  <0,  i.e.,  Re  p1Q  >  0.  Hence, 

Re  P^g  “  0. 

Theorem  3.3.5:  Let  g  be  a  widest-sense  Hurwitz  polynomial 

and  let  d  be  the  greatest  common  divisor  of  g  and  gA.  Then  d 
is  self-paraconjugate  Hurwitz. 

Proof:  In  view  of  Lemma  1  in  [3],  we  have  d*  -  Cd  where  C  is 

a  constant.  The  rest  follows  from  Theorem  3.3.1. 

Oheorem  3.3.6:  If  g(g)  is  a  widest-sense  Hurwitz  polynomial 
in  £,  then  g(£)  and  g*(£)  have  a  proper  common  factor  if  and 
only  if  g( j«)  -  0  where  the  real  k-tuple  w  can  assume  any 

value  belonging  to  a  certain  sequentially  infinite  set  of 
order  (k-1). 

Proof :  Necessity  follows  Theorem  3.3.4  and  3.3.5, 

sufficiency  from  Theorem  3.2.5. 

Theorem  3.3.7:  A  widest-sense  Hurwitz  polynomial,  g,  can  De 
written  as  a  product  of  a  scattering  Hurwitz  polynomial  and  a 
self-paraconjugate  Hurwitz  polynomial. 

Proof :  Write  g  -  ad  where  d  is  the  greatest  common  divisor 

of  g  and  g#.  Then,  a  is  relatively  prime  with  a#.  The  rest 
follows  from  Theorem  3.3.1  (item  (ii))  and  3.3.5. 

The  following  corollary  follows  immediately: 

Corallary  3. 3. 7.1;  An  irreducible  widest-sense  Hurwitz 
polynomial  is  either  scattering  Hurwitz  or  self-paraconjugate 
Hurwitz . 

Theo cem  3.3.8:  A  polynomi a  1  g  is  self-paraconjugate  Hurwitz 
if  and  only  if  g  t  0  for  Rep  >  0  and  for  Reg  <  0. 


Proof:  Necessity:  Sine*  g#  -  Cg,  a  zero  for  Re£  <  0  would 

imply  a  zero  for  Re£  >  0  and  i*  hence  excluded. 

Sufficiency:  If  g(£)  *  0  for  Re£  >  0  then  due  to  Theorem 

3.3.3  g  is  the  product  of  a  self-paracon jugate  Hurwitz  factor 
and  a  scattering  Hurwitz  factor.  However,  factors  of  the 
latter  type  are  excluded  because  they  are  known  to  have  zeros 
in  rep  <  0  (Theorem  3  in  (3]). 

Theorem  3.3.9:  A  polynomial  g  is  self-paracon jugate  Hurwitz 
if  and  only  if  all  its  irreducible  factors  are 
self-paraconjugate  Hurwitz. 

Proof :  Sufficiency  is  quite  obvious.  To  show  necessity, 

observe  that,  due  to  Corollary  3. 3. 7.1  irreducible  Zactors  of 
g  are  either  self-paraconjugate  Hurwitz  or  scattering 
Hurwitz.  Presence  of  a  scattering  Hurwitz  factor,  in  view  of 
Theorem  3  in  (3],  would  imply  that  g  has  a  zero  in  Rep  <  0, 
which  is  ruled  out  by  Theorem  3.3.8. 

The  following  corollary  follows  immediately: 

Corollary  3. 3. 9.1:  Factors  and  products  of 

self-paraconjugate  Hurwitz  polynomials  are  self-paracon jugate 
Hurwi tz . 

Corollary  3. 3. 9. 2:  Let  g  be  a  scattering  Hurwitz  polynomial 
and  h  oe  a  self-paracon jugate  polynomial.  Then,  g  and  h  are 
relatively  prime. 

Proof :  Otherwise,  g  would  contain  a  proper  self-paracon ju¬ 
gate  Hurwitz  factor  and  thus  would  not  be  relatively  prime 
with  g#.  Alternatively,  apply  Theorem  3.3.8  as  well  as 
Theorem  3  in  (31. 

Lemma  3.3.10:  Let  g  be  a  polynomial,  and  let  us  select  one 
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of  the  pj,»  say  p,  ,  and  choose  m  £  n .  -  deg^g.  Then  the 
polynoaial  g' (»»,£')  -  h*.g(n”1#£#)  i*  widest-sense  Hurwitz  if 
and  only  if  g  is  widest-sense  Hurwitz. 

Proof :  Obviously  follows  from  the  fact  that  g' ( Hg,£' g )"0  for 
Jten0  >  0,  Re£g  >  0  if  and  only  if  p2  -  Hq1  #  £'  •  £'g  is  a 
zero  of  g  in  Re£  >  0. 

Theorem  3.3.11:  Let  g  be  a  widest-sense  Hurwitz  polynomial 
in  £  and  let  i  be  any  of  the  integers  1  to  k  such  that 
n^*deg^g  >  1.  Then  3vg/8p^  ,  for  \>  »  1,2,..., is  also 
widest-sense  Hurwitz. 


Proof ;  It  is  enough  to  prove  the  theorem  for  v  »  1.  Assume 
i  -  1.  Since  the  polynomial  g^lp^)  "  9(P2/£g)  obtained  by 
freezing  £'  at  £Q  in  Re£g  >  0  is  widest-sense  Hurwitz,  by 
invoking  a  classical  result  it  follows  that  dg^/dp^  is  also 
widest-sense  Hurwitz.  The  result  then  follows  by  noting  that 
ag/ap:  -  dg1/dp1  ,  when  £  -  Eg. 


Note,  however,  that  3g/)p^  is  not  even  scattering  Hurwitz 
when  g  is  strict  sense  Hurwitz.  As  an  example,  the 
polynoaial  g  -  PiP2*pl‘fl  i*  strict  sense  Hurwitz,  but  9g/Jp2 
-  p1  is  only  a  widest-sense  Hurwitz  polynoaial. 

Let  us  write  the  polynoaial  g(£)  with  deg^  g  »  n^  as 

nl 

g  -  £  av  (£'>Pj  (3.4) 

v-0 


where  the  Av(£' )  ace  polynomials  in  £' . 

Lemma  3.3.12:  Let  g  be  a  widest-sense  Hurwitz  polynomial 
involving  one  of  the  variables,  say  p^.  If  g  is  written  as 
in  *2.4»,  then  for  any  two  integers  a,fl  with  0  <  a  <  6  n. 
and  A^eO ,  there  exists  a  set  of  positive  integers  N^,v«a,a  * 
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1,...,0,  such  that  the  polynomial  in  (3.5)  is  widest-sense 
Hurwitz . 

3 

*  Nv  Av  (£')  pp* 
v-a  (3.5) 


Proof:  Due  to  Theorem  3.3.11,  ga  -  a^/apj  is  widest-sense 

Hurwitz,  while  deg,  ga  •  (n^-a).  Consequently,  due  to  Lemma 
3.3.10,  g^  -  p^  g^( Pi”^ )  is  also  widest  sense  Hurwitz, 
with  deg^g^  -  n^-a  .  Invoking  Theorem  3.3.11  again,  it 

n.-3  n, -3 

follows  that  g^  ■  a  9i/8Pl  i#  widest-sense  Hurwitz 

with  deg.  g-  -  (n.-«)  -  (n.-3)  -  (3-a).  Therefore,  from 

A  P  +  A  —1 

Lemma  3.3.10  we  have  that  g£  -  p^p  g3^pl  *  is  widest-sense 

Hurwitz.  The  proof  is  then  completed  by  noting  that  g£  has 

the  form  (3.5). 


Theorem  3.3.13:  Let  g  be  a  widest-sense  Hurwitz  polynomial, 
expressed  as  in  (3.4),  and  u  be  any  integer  0  <  u  <  n^. 
Then  the  following  hold  true:  (i)  If  0,  then  A^  is 
widest-  sense  Hurwitz;  in  particular,  A„  ( p ' )  *  0  for  Rep'>0. 
(ii)  If  0  <;  u  ,!/♦! ,..  .//+r<  n^  and  r  >  2  Chen  it  is  impossible 


to  have  A 


*♦1 


-  A  .... -  A  ^  . ■  0,  A  0  0, 

m+2  0+v- 1  v 


n^-1 , 


if  A  •  0, 
n 


*  Qm  For  *nY  v  satisfying  1  <  u  < 

A  .  00  and  A  .  0  0 

y- 1  y*l 

If  for  any  y  satisfying  0  <  y  <  n^,  K^0  0,  A^+1  0  0 


0  0  and  A^+1  0  0  then  A^_^/A^^  is  a  positive  constant. 

(iv) 

then  A  /A  .  is  a  positive  function. 
y  y+l  r 


Proof :  (i)  Follows  immediately  by  choosing  a-0-//  in  Lemma 

‘ .3.12. 


(ii)  Me  prove  this  for  r  •  3,  the  proof  being  similar  for 
r  >  3 
-  (N 


Lemma  3.3.12,  with  a  -  y  , 0  •  y  ♦  3 ,  implies  that  g^ 
A  .  ,pi*  4N^Ay)  i»  widest  sense  Hurwitz,  where  and 

However,  this  latter  conclusion  is 


y*  3  //♦  3F1 
are  positive  integers. 


impossible  due  to  the  fact  that  for  any  fixed  g'  in  Reg’  >  0 
(where,  due  to  (i),  A^  0  0  and  A^+j  0  0)  the  cubic  equati 


on 
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g3  -  0  in  p,  has  at  least  one  solution  in  Rep^  >  0. 

(iii)  Lemma  3.3.12  with  a«//-l,(3  -  p+1  implies  that  g2  ■ 

(N^+1  A^+i  P^  +  N^_iA^_i )  is  widest  sense  Hurwitz,  where 
N^+1  and  N^_i  are  positive  integers.  We  may  claim  that 

is  real  and  positive  for  each  in  Rep'>  £  , 
because  otherwise  there  would  exist  a  £'  *  E'g  *n  ReE'o  >  — 
for  which  the  equation  g^  0  would  have  a  solution  in  Rep^  > 
0.  Since  A^_^/A^+^  *s  thus  holomorphic  as  well  as  real 
everywhere  in  Reg'  >  0,  invoking  a  standard  result  in  the 

theory  of  functions  of  complex  variables  (e.g.  the  theorem  on 
the  minimum  of  the  imaginary  part),  we  have  that  A^_^/A^+1  is 
a  constant. 

(iv)  Lemma  3.3.12  with  a-/*  ,  0-//+1  yields  that  g^  -  N^+l 
A^+iPi  +  N^A^  **  widest-sense  Hurwitz.  Clearly,  if  A^/A^+1 
has  negative  real  part  for  some  £'  a  £'g  in  Rep'g  >  0  then 

9i<Pio'E'o)  "  where  Pl0““^Ni/A//^,0^/^N^+lA/;+l(^,,0)  1  * 

Since  Rep1Q  >  °»  the  latter  conclusion  is  impossible  due  to 
the  widest-sense  Hurwitz  property  of  g^. 

Theorem  3.3.14;  Let  g  be  a  widest-sense  Hurwitz  polynomial 
in  £  and  g^  be  the  ( k-1 )-variable  polynomial  obtained  by 

freezing  any  one  of  the  variables,  say  p^,  at  p^  -  p^  in 

Rep^  >  0.  Then  degig  -  deg^g^  for  the  remaining  variables 
p^,  i  -  2  to  k  except  possibly  for  finitely  many  values  of 
p1(J  on  Repx  -  0. 

Proof :  Write  g  as  a  polynomial  in  £'  whose  coefficients  are 

polynomials  in  p^.  Let  A  -  A(p^)  be  any  one  of  these 

coefficients  that  is  not  identically  zero  for  all  p^.  By 

repeated  application  of  Theorem  3.3.13,  item  (i),  it  follows 
that  A(p^)  is  widest-sense  Hurwitz.  Hence,  A(p1q)  ^  0  for 
Re  p^Q  >  0  and  there  are  at  most  finitely  many  p^Q  -  j^g  for 
which  A(jw^g)  -  0  (cf.  also  Lemma  3.3.2). 


3.4.  Properties  of  scattering  Hurwitz  polynomials. 


Theorem  3.4.1:  If  g  is  a  scattering  Hurwitz  polynomial  then 
g  cannot  have  a  zero  for  a  £g  with  Rep^g  ”  0  for  one  of  the  i 
and  Rep.^  >  0  for  the  remaining  i. 

Proof :  We  may  assume  Rep^g  ■  0  for  i  -  1  i.e.,  p-^g  -  3wio* 
In  view  of  Lemma  3.3.2,  a  zero  of  the  type  just  mentioned 
would  require  g(jw^g,  £' )  ■  0  for  all  £' .  Thus,  in  view  of 
Theorem  3.2.10,  (p^-ja^g)  would  be  a  factor  of  g.  Since 
(p^-j<a^o)  is  self-paracon jugate ,  this  is  excluded  due  to 
(3.3). 


It  is  important  to  note,  however,  that  a  scattering  Hurwitz 
polynomial  can  indeed  have  zeros  e.g.  for  p^  -  3wi'P2  "  3U2 
and  Rep£  >  0,  i  -  3  to  k.  Consider  the  polynomial  g  -  P1P2  + 
P2P3  +  P3P1  +  P1P2P3'  which  is  scattering  Hurwitz  11],  but 
g-0,  when  p1  -  p2  -  0. 


Theorem  3.4.2:  Let  g  be  a  scattering  Hurwitz  polynomial  and 
let  i'  be  any  specific  one  of  the  i  •  1  to  k.  Then  the 
( k-1 )-variable  polynomial  g^  obtained  by  freezing  pif  at  pif 
-  jfti^ , g  is  also  scattering  Hurwitz  and  has  the  same  partial 
degrees  in  the  remaining  variables  as  g,  with  the  possible 
exception  of  at  most  finitely  many  values  of  u^,g. 


Proof :  We  may  assume  i'  •  1.  Referring  to  Theorem  3.4.1, 
g^g')-  g(j«^Q,£')  k  0  in  Re£'  >  0.  Furthermore,  g(£)  and 
g#(£)  do  not  have  any  nonconstant  common  factor.  Therefore, 
due  to  Theorem  3.2.7  with  m-1,  the  polynomials  g^lp')  and 
9]>*(£')  •  g*(  j«10,  -  £'*)  -  9*(j<*>1g»£')  can  have  a  common 
nonconstant  factor  for  at  most  finitely  many  values  of  w^g. 
Furthermore,  in  view  of  Lemma  3.2.6,  a  lowering  of  a  partial 
degree  can  also  occur  at  most  for  finitely  many  values  of 


w 
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In  view  of  the  above  theorem  and  Definition  3.2.1a  it  is  also 
possible  to  state  the  more  general  result: 

Theorem  3.4.3:  Let  g  be  a  scattering  Hurwitz  polynomial. 
For  any  m  such  that  1  <  m  <  k  consider  the  ( k-m) -variable 
polynomial  g'  obtained  by  freezing  m  of  the  p^  at  p^^  -  jwiQ 
for,  say,  1  <  i  <  ml  Then  there  exists  a  sequentially  almost 
complete  set,  &m,  of  order  m  of  m-tuples  such  that  for 
< ***1 0 , ***2 O'... “mO )  c2m'  scattering  Hurwitz,  with 
the  same  partial  degrees  in  Pm+^  to  pk  as  g.  Furthermore, 
any  ordering  may  be  chosen  for  the  set  &m. 

Proof:  Proof  of  the  above  theorem  follows  by  sequentially 
freezing,  in  any  order,  the  variables  pi  to  pm  on  the 
imaginary  axis  and  observing  Theorem  3.4.2. 

Theorem  3.4.4:  A  polynomial  g  in  k  variables  is  scattering 
Hurwitz  if  and  only  if  (i)  g  is  widest-sense  Hurwitz,  and 
(ii)  the  set  of  real  k-tuples  w  such  that  g(j«)  -  0  does  not 
form  a  sequentially  infinite  set  of  order  ( k— 1 ) . 

Proof :  Follows  immediately  from  Definition  3.2.4,  and 
Theorem  3.3.6. 

Lemma  3.4.5:  Let  g  be  a  scattering  Hurwitz  polynomial. 
Freeze  one  of  the  pi#  say  p.,,  at  p^g  with  Repi,fl  >  0.  The 
resulting  polynomial,  g^,  in  the  remaining  p^  is  still  a 
scattering  Hurwitz  polynomial,  and  the  partial  degrees  of  g^ 
in  these  remaining  p^  are  the  same  as  for  g. 

Proof :  Without  loss  of  generality  we  may  assume  that  i'-l 
i.e.,  g^  »  g^(g').  Obviously  g^  is  widest  sense  Hurwitz 
(Theorem  3.3.1).  We  show  in  the  following  that  g^  and  g^ 
are  relatively  prime.  If  g^  and  g^*  have  a  nontrivial  common 
factor,  then  due  to  Theorem  3.3.6  there  exists  a 
sequentially  infinite  set,  2'  of  real  (k-l)-tuples  u>'  of 


order  (k-2)  such  that  g(p^g,  «' )  ■  0  for  all  w'  c  Q' .  Thus, 
applying  Theorem  3.3.3,  item  (i),  for  m  -  k-1,  we  conclude 
that  9(p^#jw')  ■  0  for  all  p^  and  all  u' c  S'.  We  therefore 
have,  in  particular,  gfjw^jw'  )  -  0  for  all  w'e2'  and  any 
arbitrary  Thus  g  would  be  zero  for  real  k-tuples  « 

belonging  to  a  sequentially  infinite  set  of  order  (k-1), 
which  is  impossible  in  view  of  Theorem  3.4.4.  Furthermore, 
the  preservation  of  the  partial  degrees  follows  from  Theorem 
3.3.14. 

Repeated  use  of  Lemma  3.4.5  yields  the  following  more 
general  result: 

i 

Theorem  3.4.6:  If  g  is  a  scattering  Hurwitz  polynomial  in  £ 
then  for  1  <  m  <  k  the  m-variable  polynomial,  g' ,  obtained  by 
freezing  (k-m)  of  the  pj^  at  p^^  -  p^  in  Rep^^  >  0  is 
scattering  Hurwitz  and  the  partial  degrees  of  g'  are  the  same 
as  the  corresponding  partial  degrees  of  g. 

Theorem  3.4.7:  If  g  is  a  scattering  Hurwitz  polynomial 

expressed  as  in  (1)  with  deg^g  -  n^,  then  Av  M  0  for  each 
v  ■  0,1,2, ...n^. 

Proof :  Assume  i  -  1.  For  any  £'  -  £'g  with  Re£'g  >  0,  the 

polynomial  g^p^)  "  9(Pi'Eq)»  due  t0  Theorem  3.4.6,  is  a 

scattering  Hurwitz  polynomial  in  the  single  variable  p^ 
(i.e.,  a  Hurwitz  polynomial  in  the  classical  sense  [3])  of 
degree  equal  to  n^  -  deg^g  and  thus  has  nonzero  coefficients 
for  v  -  0  to  n^. 
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3.5.  Rational  Positive  Functions  and  Related  Results 


Theorem  3.5.1:  Let  F-h/g  be  a  rational  positive  function  in 
irreducible  form  (g  and  h  thus  being  relatively  prime 
polynomials)  with  h(£)  jt  0  (and,  obviously,  g(£)  M  0)  .  Then, 

(i)  the  polynomial  g+h  is  scattering  Hurwitz, 

J 

I 

(ii)  both  g  and  h  are  widest-sense  Hurwitz. 

Proof :  (i)  Consider  the  function 

p  -  ( F-l )/( F+l )  (3.6) 

which  we  can  write  in  the  form  p«c/d,  c  and  d  being  defined 
by  c-h-g,  d-h+g.  Clearly,  c  and  d  are  relatively  prime. 
Then,  for  Re  £>0^  and  d(£)  *  0,  we  have  Re  F ( £ )  >^0  and  thus 
|  p ( £ )  |  <_1  if  g(£)*0,  and  p(£)-l  if  g(£)-0.  Hence,  by  Theorem 
1  in  [3]  (where  in  view  of  the  proof  of  Theorem  1  in  [1],  on 
which  the  proof  of  Theorem  1  in  (3]  is  based,  (3.4)  should  be 
interpreted  to  mean  that  the  first  inequality  is  known  to 
hold  for  all  those  £  in  Re£>0  for  which  g(£)#0),  d  is 
scattering  Hurwitz.  (Note  that  d(£)#0  for  Re£>0  follows  from 
Theorem  3.2.13. ) 

(ii)  Since  d(£)#0  in  Re  £>0,  we  cannot  have  g(£)-h(£)-0  in 
Re  £>0.  If  only  one  of  the  polynomials  g  and  h  has  a  zero 

for  a  £q  with  Re  £g  >  0  we  have  p-+l .  This,  however,  is 
excluded  since  due  to  the  maximum-modulus  theorem  [10]  we 
have  in  fact  |p(£)|<l  in  Re  £>0,  except  if  p  is  a  unimodular 
constant,  i.e.,  if  F-jC,  C  being  a  real  constant  (in  which 
case  g-1 ,  h- jC ) . 

Theorem  3.5.2:  If  F  is  a  rational  positive  function,  then 
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(i)  Re  F(£)>0  if  Re  £>0  and  F  is  nontrivial; 

(ii)  Re  F(£)>0  for  any  £  in  Re  £>0  where  F(£)  is  holomorphic; 

(iii)  assuming  F{£)^0,  1/F  is  also  a  rational  positive 
function . 

Proof :  (i)  Define  p  by  (3.6)  and  consider  a  closed 

polydomain  U  in  the  neighborhood  of  an  arbitrarily  selected 
point  £q  with  Re  £g>0.  In  view  of  Theorem  3.5.1,  F  and  p  are 
holomorphic  in  U.  Thus  |p(p)|^l  for  £  e  0,  while  application 
of  the  maximum-modulus  theorem  in  its  simplest  form  [10] 
yields  |p(Pq)I<1»  i.e..  Re  F(£q)>0.  The  proof  of  (ii) 

follows  by  simple  continuity  arguments,  and  that  of  (iii) 
follows  in  an  obvious  fashion  from  (i). 

Theorem  3.5.3:  If  g  and  h  are  relatively  prime  polynomials, 
then  F-h/g  is  a  positive  function  if  and  only  if  (i)  the 
polynomial  d-g+h  is  scattering  Hurwitz,  (ii)  Re  F(j«)^0  for 
all  £-jw  where  F(£)  is  holomorphic. 

Proof :  Necessity  of  (i)  follows  from  Theorem  3.5.1, 

necessity  of  (ii)  from  Theorem  3.5.2  item  (ii).  For  proving 
sufficiency,  observe  that  for  p  defined  again  by  (3.6)  and 
d(jw)#0,  we  have  Re  F(ju)>0  and  thus  |p(jw)|£l  if  g(j«)*  0, 
and  p(ju)-l  if  g(jw)»0,  altogether  thus  |p(j«)|<l  wherever 
p(jw)  is  holomorphic.  Hence  applying  Theorem  1  in  [2] 
(generalized  maximum-modulus  theorem),  we  conclude  that  in  Re 
£>0  we  have  | p ( £ ) | < 1  and  thus  also  Re  F(£)>0.  The  only 
exception  to  this  is  if  p  is  a  unimodular  constant,  in  which 
case  F  is  an  imaginary  constant,  thus  a  trivial  positive 
function . 

lemma  3.5.4:  Let  F  be  a  positive  function  (reactance 

function).  For  any  m  such  that  1  <  m  <  k,  consider  the 

( k-m ) -va r l able  rational  function  F'  obtained  by  freezing,  in 

F,  m  of  the  variables  pi  at  p^  -  ju^g,  3<ay  for  1  *  1  t0  m- 

Then  there  exists  a  sequentially  almost  complete  set,  Q  ,  of 

m 
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order  a  of  real  m-tuples  such  that  for  any  •  •  ••**(}*  c 

S  ,  F'  is  still  a  positive  function  (reactance  function).  If 

a 

in  addition,  F  is  given  in  irreducible  fora,  F-h/g,  then 
can  be  chosen  such  that  F'-  h'/g'  is  also  in  irreducible 
fora,  h'  and  g'  being  obtained  by  applying  the  corresponding 
freezing  operation  to  h  and  g,  respectively,  and  that  the 
partial  degrees  of  h'  and  g'  in  the  reaaining  variables  are 
the  same  as  for  h  and  g  respectively,  with  fia  reaaining 
otherwise  as  stated. 

Proof :  We  need  to  prove  the  leaaa  for  a  -  1  only.  Due  to 

Theorem  3.5.1,  g  is  videst-sense  Hurvitz.  Also,  let  g'(g') 
-g(  joi^Q .  Due  to  Leaaa  3.3.2,  there  exists  an  alaost 
complete  real  set,  21,  such  that  g'(g')  *  0  for  Reg'  >  C  and 
ui1  q  e  2^  and  that  the  partial  degrees  of  h'  and  g'  in  p2  to 
p^  are  the  same  as  those  of  h  and  g,  respectively.  In 
particular,  for  e  Q£,  F  is  regular  for  Reg'  >  0,  and  by 

invoking  Theorem  3.5.2,  it  follows  that  the  function  F'ig'»  • 
F(j«i0,  2'),  for  c  2J ,  satisfies  the  property  that 

ReF'(g')  >  0  for  Reg'>0,  i.e.,  it  is  a  positive  function. 

Furthermore,  if  F  is  a  reactance  function  then 
F  +  F#  -  0  implies  that  F'  ♦  F'  •  0.  Consequently,  F'  is  a 
reactance  function. 

Also,  due  to  Theorem  3.2.7,  there  exists  an  almost 
complete  real  set  2£  such  that  the  polynoaials  g'(g')  and 
h'(g')  -  h( j«10 ,g' )  are  relative  priae  for  c  . 

Therefore,  h'/g'  is  in  irreducible  fora  if  w  is  required  to 
be  in  the  almost  complete  set  2^  -  (Q^nQJ). 

Lemma  3.5.5a:  The  numerator  and  denominator  polynomials  of  a 
positive  function  in  irreducible  form  cannot  contain 
self-paracon jugate  factors  of  multiplicity  larger  than  one. 

Proof:  Let  F  ■  h/g  be  a  positive  function  in  irreducible 
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fora,  and  g  contain  a  propar  salf-paracon jugate  factor,  a,  of 
aultiplicity  v.  Lat  a  involve,  say,  p^ •  Than  it  follow*  by 
invoking  Leaaa  3.5.4  with  a  -  k-1  that  there  exists  a 
sequentially  almost  complete  real  set,  8',  of  order  'k-1; 
•uch  that  for  any  real  k-tuple  m' q  c  8’,  the  rational 
function  -  hj/gj.  with  h^Pj)  -  h(pl,3w'0.i  and  gl'pl;  - 
g(plf  3a'0 1 ,  i*  a  ona-veriable  positive  function  in 
irreducible  form  and  that  the  partial  degree  of  9,  in  is 
the  same  as  that  of  9.  Obviously,  -  a'p^.^w'g  is  a 
nonconatant  factor  of  g^'p^i  multiplicity  at  least 
equal  to  v.  Since  a  is  sel f -paraconjugate ,  a ^  is  also  self- 
paracon jugate .  Furthermore,'  taking  into  account  Theorem 
3.3.3,  a^  is  widest-sense  Hurwitt  and,  therefore,  its  teres 
are  restricted  tc  be  on  p,  »  Consequent./  9.  has  teres 

a  *  * 

of  multiplicity  at  least  «  on  the  p.  •  'm.  axis  Since  9.  is 
Che  denominator  of  a  positive  function  in  irreducible  form, 
the  latter  conclusion  dictstes  that  *  •  1  Similar  arguments 
hold  for  h  e  9  .  in  view  c?  Theorem  r  c  2  .'em  ... 


Lemma  355b  If  f  •  1  9  is  a  pcs.*.-e  furc'.cr 
polynomials  then  deg.h  -1  deg  9  '  deg,*’  ♦  .  for 
1.2,  x 


and  z 


Proof 


The  p 
d 
c 


1 

^  1  - 

f 


invoking  a  wall  known  ona-variabla  rasult  wa  assart  that 
(dg^/dp^/g,  is  a  positive  function.  However,  since  this 
latter  conclusion  is  true  for  any  g£  in  Reg'  >  0,  we  have 
Re( ( )/g ]  >  0  for  Reg  >  0. 

Finally,  for  the  sake  of  completeness ,  we  offer  a  sore 
complete,  but  obvious  version  of  Theorem  3.5.2,  item  (i). 

Theorem  3.5.7:  A  rational  function  F  is  a  nontrivial 

positive  function  positive  if  and  only  if  ReF(g)  >  0  for  Reg 

>  2  • 

Fr oof :  Sufficiency  follows  by  a  simple  continuity  argument. 

Necessity  has  been  shown  in  Theorem  3.5.2,  item  (i),  but  it 
also  follows  by  simple  application  of  the  theorem  on  the 
minimum  of  the  real  part  of  a  holomorphic  function. 
Alternatively,  we  may  freete  g'  at  a  g£  with  Reg£>0,  in  which 
case  r  reduces  to  a  function  r,  in  p^  alone.  Applying 
Theorems  3  3.1,  3.314,  and  3.5.1,  F,  is  found  to  be  a 

i. 

nontrivis*  positive  function,  i.e.,  the  proof  is  reduced  to 


3.6.  Reactance  Hurwitz  Polynomials 


The  following  lemma  follows  directly  from  Definition  3.2.6 
and  Theorem  8  in  (3]: 

Lemma  3.6.1:  A  reactance  Hurvitz  polynomial  is 
self-paracon jugate  Hurwitz. 

Definition  3.2.6  is  justified  by  the  following  theorem: 

Theorem  3.6.2:  1.  If  g  is  a  reactance  Hurwitz  polynomial, 
there  exists  a  polynomial,  h,  relatively  prime  with  g,  such 
that  h/g  is  a  reactance  function  in  irreducible  form. 

2.  Vice  versa,  if  h/g  is  a  reactance  function  in  irreducible 
form,  the  following  holds:  (i)  g  and  h  are  reactance  Hurwitz 
polynomials.  (ii)  Por  any  constant  C,  the  polynomial  d 
defined  by  d-gg+hg,  gg-Cg,  hg-Ch,  is  scattering  Hurwitz.  In 
particular,  it  is  always  possible  to  choose  C  in  such  a  way 
that  the  paraeven  and  paraodd  parts  of  d  are  equal  to  gg  and 
hg,  respectively,  or  that  these  parts  are  equal  to  hg  and  gQ, 
respectively . 

Proof :  Por  proving  the  first  statement,  observe  that  in  view 
of  Definition  3.2.6  there  exist  polynomials  gQ  and  hg  such 
that  g0.-i9g,  hg.-Jhg,  and  gQ  -  Cg  and  that  g0+hg  is 
scattering  Hurwitz,  the  two  upper  and  the  two  lower  signs 
corresponding  to  one  another  and  C  being  a  nonzero  constant, 
•y  Theorem  7  of  [3],  hg/g0  is  then  a  reactance  function  in 
irreducible  fora  and  the  sane  is  thus  true  for  h/g  where 
h-hg/C. 

Por  proving  the  second  statement,  observe  first  that  in  view 
of  Theorem  3.5.1,  g  and  h  are  widest-sense  Hurwitz  while 
9g*h0  is  scattering  Hurwitz  where  gg-Cg  and  hg-Ch,  C  being  an 
arbitrary  nonzero  constant,  and  that  in  view  of  Definition 
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3.2.9  we  have  h*/g#—h/g,  thus  g*-Y9#  h#»-rh,Y  b«in9  • 
constant  (nscsssarily  uniaodular).  If  in  particular  we 
choose  C  such  that  C2  -♦  r  we  have  gQ*-+g0,  ho*"*hO' 

• 

Theorem  3.6.3;  If  9  is  a  self-paracon jugate  Hurwitz 
polynoaial  then  for  any  i  -  1  to  k  with  -  deg^g  >  1,  f  • 
h/g  is  a  reactance  function  where  h  -  ig/ap^. 

Proof :  Proa  Theorea  3.5.6  it  follows  that  h/g  is  a  positive 
function.  Assuae  i  -  1.  If  g  is  written  as  in  (3.4),  then 

9*  -  I1  A  <E'>*  (-Pi>V- 

v-0 

Since  g  is  sel f-paracon jugate ,  g  -  Cg ,  for  soae  constant  C 
with  |C|  «  1.  Therefore,  -  (-1)VCA%1#  whence  it  can  be 
shown  that  Jg/ap^— C(  ag/ap^ )  # .  It  thus  follows  that  r#  —  r. 
Consequently,  r  is  a  reactance  function. 

Corollary _ 3 . 6 . 3 . 1 ;  If  g  is  an  irreducible, 

self-paracon jugate  Hurvitz  polynoaial  then  for  any  i  •  1  to  k 
with  n ^  -  deg^g  £  1,  (ag/ap^/g  is  a  reactance  function  in 
irreducible  fora. 

Proof :  Since  g  is  irreducible  and  the  partial  degree  of 
*9/3Pl  1*  smaller  than  that  of  g  in  the  variable  p^  the 
polynoaials  3g/3p^  and  g  are  relatively  priae.  The  proof  is 
then  completed  by  observing  Theorem  3.6.3. 

Theorem  3.6.4:  A  polynomial  g  is  reactance  Hurwitz  if  and 
only  if  all  its  irreducible  factors  arc  self-paracon jugate 
Hurwitz  and  it  contains  no  multiple  factors. 

Proof : _ Necessity :  If  g  is  a  reactance  Hurwitz  polynomial, 

invoking  Theorem  3.6.2  there  exists  a  reactance  function, 
thus  a  positive  function,  F»h'g,  such  that  h  and  g  are 
relatively  prime.  In  view  of  Lemma  3.6.1  and  Theorea  3.3.9, 
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irreducible  factors  of  9  are  necessarily  self-paracon jugate . 
Zt  then  iaaediately  follows  froa  Leaaa  3.5.5a  that  9  cannot 
contain  aultiple  factors. 

Sufficiency:  Let  g^,  i  -  1  to  m,  be  the  irreducible, 

distinct,  self-paracon jugate  factors  of  g.  Let  g^  involve 
the  variable  p[c(pL ,p2, • • • ,Pk)  end  consider  the  polynoaial  h 

v 

defined  by  h/g  -  l  (hl/gi),hi  -  *g*/*Pj • 

froa  Corailary  3.6. 3.1  we  conclude  on  the  one  hand  that  h/g 
is  a  reactance  function,  and  on  the  other  that  each  h^  is 
relatively  priae  with  the  corresponding  g^.  Thus,  since  by 
assuaption  the  g^  are  pairwise  mutually  priae,  h  is 
relatively  priae  with  g.  The  proof  is  then  coapleted  by 
observing  Theorea  3.6.2,  second  part,  itea  (i) 

Corollary  3. 6. 4.1:  Products  of  reactance  Hurwitz 

polynomials,  that  are  pairwise  relatively  prime  are  reactance 
Hurwitz  polynomials.  Conversely,  any  factor  of  a  reactance 
Hurwitz  polynomial  is  also  a  reactance  Hurwitz  polynomial. 

Proof :  T ollowa  immediately  froa  Theorem  3.6.4 

Corallory  3. 6. 4. 2:  A  polynomial,  g,  is  reactance  Hurwitz  if 
and  only  if  it  is  sel f-pa racon jugate  Hurwitz  and  contains  no 
aultiple  factor. 

Proof :  rollows  froa  Theorems  3.3.9  and  3.6.4. 

Corollary  3. 6. 4. 3:  Any  scattering  Hurwitz  polynomial  is 

relatively  priae  with  any  reactance  Hurwitz  polynoaial. 

Proof  :  Tollows  from  Corollary  3.  3.9.2  and  Lemma  3.6.1. 

Theorem  3.6.5;  Let  g  be  a  reactance  Hurwitz  polynomial.  For 
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•ay  integer  a  such  that  1  j  a  <  k,  freese  a  of  tho  variables 
at  p^  •  aay  for  i  •  1  to  a.  Consider  the  polynoaial 

9'  to  pk  daftnad  by  gMp^ . Pk>-  g<  j*10 . j*a0. 

P»» 1 * ‘ . • *  Pk ) *  Than  thara  exists  a  sequentially  alaost 
coaplata  set,  Q^,  of  order  a  of  real  n- tuples  such  that  for 
(«^0««20'  *  *  '*n0  *  *fla'  i#  #tiii  reactance  Hurwitz  and  that 
its  partial  degrees  in  Pa4.^  to  pk  are  the  saaa  as  for  g. 

Proof :  Due  to  Theorea  3.6.2  thars  exists  a  polynoaial  h  such 
that  r  •  h/g  is  a  reactance  function  in  irreducible  fora. 
Applying  Leaaa  3.5.4  to  the  rational  function  r  •  h/g,  we 
conclude  that  aay  be  chosen  in  such  a  way  that  the 
resulting  P'  *  h'/g*  is  a  reactance  function  in  irreducible 
fora,  with  g'  having  the  saae  partial  degrees  in  to  pk 
as  g.  The  rest  of  the  proof  follows  by  applying  itea  (i)  of 
the  second  part  of  Theorea  3.6.2  to  P‘ . 
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iaaittance  Hurwitz  polynoaials  as  given  in  Definition  3.2.7. 


Theorea  3.7.1:  If  P  -  h/g  is  a  positive  function  in 
irreducible  fora,  then  both  h  and  g  are  iaaittance  Hurwitz 
polynoaials . 

Proof :  g  is  widest  sense  Hurwitz  in  view  of  Theorea  3.5.1. 
Invoking  Theorea  3.3.7  it  follows  that  g  is  product  of  a 
scattering  Hurvits  polynoaial,  g^ ,  and  a  self-paracon jugate 
Hurwitz  polynoaial,  g^.  ly  virtue  of  Corollary  3. 3. 9.1, 
Leaaa  3.5.5a,  and  Corollary  3. 6. 4. 2,  g2  is  reactance  Hurwitz. 
The  saae  arguaent  holds  for  h. 

Corollary  3. 7. 1.1:  Pactors  of  iaaittance  Hurwitz  polynoaials 
are  iaaittance  Hurwitz.  Conversely,  products  of  iaaittance 
Hurwitz  polynoaials  that  do  not  have  any  coaaon  self- 
paracon  jugate  factors  are  also  iaaittance  Hurwitz. 

Proof :  Obviously  follows  froa  Definition  3.2.7  and  Corollary 
3  .  i  .  4 . 1 . 


Theorea  3.7.2;  (very  iaaittance  Hurwitz  polynoaial  is  the 
nuaerator  or  denoainator  of  a  positive  function  in 
irreducible  fora. 


Proof :  Since  the  reciprocal  of  a  positive  function  is  also  a 
positive  function,  it  is  enough  to  prove  the  theorea  for  the 
denoainator  polynoaial,  9.  Let  g  -  ab,  where  a  and  b  are, 
respectively,  the  scattering  Hurwitz  and  the  reactance 
Hurwitz  factors  of  g.  Then,  due  to  Theorea  3.6.2,  there 
exists  a  polynoaial  c  that  is  relatively  prise  with  b  and 
such  that  :  b  is  a  reactance  function.  L»t  us  define  F  -  * 
(a.,  a  1  *  c  bi.  T  can  also  be  written  in  the  fora  F  •  h  g, 


(i’- vAi*"  .1/ .y  ^  < 


vV, 


where  h  -  ab  ♦  a#b  +  ac.  Since  aft/a  is  an  all-pass  function, 
1  ♦  (a#/a)  is  a  positive  function  (cf.  Lena  2  of  [3]),  and 
the  same  is  true  of  r.  According  to  Corollary  3.6.4. 3,  a  and 
b  are  also  relatively  prime.  Furthermore,  in  view  of  the 
definition  of  a  scattering  Hurvitz  polynomial,  a  is 
relatively  prime  with  a*.  Hence,  since  h  -  a(b+c)>aAb  - 
{ a+a# ) b+ac ,  h  can  not  be  divisible  by  any  factor  of  a  or  of 
b.  The  proof  is  thus  completed  by  observing  that  any  common 
irreducible  factor  of  h  and  g  would  have  to  be  a  factor  of 


3.8.  Teats  for  Positive  function*  and  Related  Results 

Results,  which  are  useful  in  identifying  positive  functions 
based  on  their  behavior  on  the  distinguished  boundary  of  the 
doaain  of  holoaorphy  have  already  been  proved  essential  [11] 
in  probleas  related  to  aultivariable  network  synthesis  and 
are  elaborated  in  this  section. 

Theorea  3.9.1:  Assuae  that  the  rational  function  F  -  h/g 
satisfies  the  following  properties:  (i)  g  is  a  scattering 
Burwits  polynoaial  (ii)  deg^h  £  deg^g,  i  •  1  to  k; 

(ill)  Rer(jw)  >  0,  where  a  is  any  real  k-tuple  such  that  F(g) 
is  holoaorphic  at  g-j«.  Then  F  is  a  positive  function. 

Note  that  we  have  not  assuaed  the  polynoaials  h  and  g  to  be 
relatively  priae. 

Proof:  The  validity  of  the  theorea  for  k  •  1  is  classically 
known.  To  prove  the  result  in  the  general  case,  via 
induction  on  the  nuaber  of  variables,  we  assuae  that  the 
theorea  be  true  for  k-1  variables,  where  k  £  2. 

Let  us  freese  one  of  the  variables,  say  p^,  at  j«10  and 
define  the  rational  function  in  k-1  variables,  F^  -  h^/g^, 
where  h^tg' )  •  h(j*1Q,£')  and  g^g')  •  g(j«10,g').  Due  to 
Theorea  3.4.2,  there  exists  an  alaost  coaplete  set,  8^,  of 
real  nuabers  such  that  for  all  c  8^,  9^(£')  is  scattering 
Burwitx  and  deg^g  -  deg^^,  thus  deg^h^  <  deg^g^,  for  all  i  • 
2  to  k. 

Since  gx( j«' )  -  g(jw^Q,w' ).  the  condition  g^(jw')  *  0  replies 
g(jw1Q,jw' )  *  0.  Thus,  ReF1  ( j«' )  •  ReF( , j«' )  >  0  for 
gl( jw' )  »  0 . 

All  the  p  r e  r equi sites  for  the  validity  of  the  present  theorem 
are,  therefore,  satisfied  by  the  ( k-1 > -var table  rational 
function  F^.  Hence,  by  induction  hypothesis,  ReF^ig'i  £  0 
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for  R«2'  >  0.  Therefore,  if  «10  c  Qj  then  (3.7)  applies. 

ReF(  j«10»  £' )  >.  0  for  *11  R*£'  >  0.  (3.7) 

Let  us  now  freeze  the  variables  £'  at  an  arbitrary  point  £g 
with  Re£g  >  0  and  define  the  rational  function  Fq  *  h0/90  in 
the  variable  pj^  only,  where  hgtp^)  -  h(p^,£’g)  and  gQ  ( )  - 
g(p,,,g^).  Due  to  Theorem  3.4.6  with  m-k-1  the  polynomial  gQ 
is  scattering  Hurwitz,  thus  Hurwitz  in  the  classical  sense 
(cf.  Theorem  2  of  [3])  and  deg^gg  -  deg^  ^  deg^h  ^  deg^hg. 
Furthermore,  it  follows  from  (3.7)  that  ReFg(ju^)  >  0  for  co^ 
c  8^.  Hence  it  follows  from  a  classical  result  that  ReFg(p1) 
^  0  for  all  p^  with  Rep^  >  0.  Therefore,  we  conclude  that 
ReF (£)  £  0  for  Rep  >  0. 

The  following  comments  on  the  above  result  are  in  order.  In 
the  one-variable  case  the  proof  of  the  above  result  follows 
by  invoking  the  maximum  modulus  theorem  on  the  function 
exp(-F).  In  the  multidimensional  situation,  however,  a 
maximum  modulus  theorem  which  allows  for  special  types  of 
singularities  on  the  boundary  of  the  domain  of  holomorphy  is 
not  found  in  the  literature.  Reference  [2]  gives  a  version 
of  maxi mum- modulus  theorem  where  non-  essential  singularities 
of  the  second  kind  are  allowed  to  occur  on  the  boundary  of 
the  domain  of  holomorphy,  but  the  proof  is  restricted  to 
rational  functions  only.  The  non-  rationality  of  the 
function  exp  (-F)  makes  it  impossible  to  use  results  of  [2!, 
in  the  present  context,  thereby  calling  for  an  independent 
ptoof  of  Theorem  3.8.1. 

Note  that  Theorem  3.8.1  can  be  generalized  to  include  the 
possibility  that  the  domain  of  holomorphy  of  F  be  a  cartesian 
product  of  domains  other  than  half-planes  (e.g.,  discs)  in 
the  variables  p^.  i  -  1  to  k . 

The  following  partial  results  are  cf  seme  interest. 


Theorem  3.8.2a:  If  F  -  d/eg  is  a  positive  function  in  £ 

written  in  irreducible  fora  where  g  is  a  scattering  Hurwitz 
polynomial,  c  is  a  reactance  Hurwitz  polynomial  that  is  the 
product  of  reactance  Hurwitz  polynomials  in  one-variable 
only,  and  d  is  a  polynomial  such  that  deg^d  <  deg^cg),  i  -  1 

n 

to  k,  then  r  can  be  decomposed  as  r  ■  Z  Fv  ♦  (d^/g),  where 
each  rv  is  a  one-variable  reactance  fiXnition  and  (d^/g)  is  a 
positive  function  in  £,  with  deg^d^  <  deg^g  for  each  i  •  1  to 
k . 

Proof ;  write  c  -  cic2***cn'  where  cv  •  <  P^^vO * '  v  "  1  t0 
n,  each  p^  being  one  of  the  p^  to  pk  and  two  w^q  being 
necessarily  distinct  if  the  corresponding  p^  represents  the 
same  p^.  We  claim  that  the  rational  functions  R^  defined  in 
(3.8a)  are  constants.  We  show  this  for  v-1,  assuming  p^  - 
p1#  in  which  case  could  be  a  function  of  £' . 

lie  n 

VE'*  “  <  <P;-Iwv0,d/C9l  ;  C  “  d-c9  *  Kv/cv 

vml  ( 3 . 8a , b ) 

Por  each  £'  -  £'  in  He££  >  0,  )  is  the  residue  of  the 

positive  function  r'(pi)  -  P<Pj'Eo*  *t  the  pole  pi  -  and 

is  hence  positive  and  thus,  in  particular,  real  and  finite. 

Therefore,  R^(£')  is  real  and  holomorphic  in  Re£'  >  0. 

Consequently,  invoking  a  standard  result  from  the  theory  of 
functions  of  complex  variables,  it  follows  that  R^  is 
independent  of  £' ;  it  is  thus  a  positive  constant.  Consider 
next  the  polynomial  C(£)  as  defined  in  (3.8b).  Substituting 
for  the  R v  from  (3.8a)  in  (3.8b)  it  follows  that  C  is  zero 
for  p^  •  independently  of  the  values  of  the  other  p  . 

Hence,  all  c^  divide  C.  i.e.,  since  the  c^  are  distinct,  - 

C/c  is  a  polynomial  in  £.  Equation  f3.9a)  than  follows  by 

straightforward  algebraic  manipulation,  where  c '  : s  as 

de f i ned  i n  3  .  ?b  . 

n 

r  •  (C’,'CJ  ♦  !d,  g  ;  c '  •  1  i  R  c  :  c  3.9a,fc 

1  ,  V  V 

V- 1 


-75- 


J. 


wrm  wwwt  w  pj  w.  mw!  w  bwwb  w  w  evrni1  w;  ev.wr.  .r  .■.r^Fy 


Clearly,  c'/c  is  a  reactance  function.  Next,  we  note  that 
for  £  -  j«,  u  real,  Retd^/g)  ■  ReF  >  0  at  all  those  points 
where  c(jw)g(jw)  *  0  and  thus,  by  continuity,  where  g 1 ju>  * 
0.  Due  to  (3.9)  and  F  -  d/gc  the  inequalities  deg,d  < 
degi(gc)  imply  that  degid1  c  deg^  for  each  i  -  1  to  k . 
Hence,  invoking  Theorem  3.8.1  it  follows  that  d^g  is  a 
positive  function.  The  proof  is  thus  complete. 

Note  that  the  above  result  can  be  easily  extended  to 
the  case  with  simple  poles  at  p^  -  •  as  follows. 


Theorem  3.8.2b:  If  F  ■  d/g  is  a  rational  positive  function 
such  that  the  polynomial  g  is  scattering  Hurwitz  and  dec.d  > 
deg^g  for  some  of  the  id,  then  there  exist  nonnegative 
constants  such  that  F  can  be  written  in  the  fora 


k 

F  •  Fn*  r  K.p.  where  Fn 

u  i-1  1  i  u 


dg/g  is  a  positive  function  with. 


degidQ  <  degig.  i  -  1  to  k . 


Proof ;  Assume  first  deg^d  >  deg^g  for  i  •  1.  Let  A  £'  and 

*(£'  )  be  the  respective  leading  coefficients  of  d  and  g.  wr.en 

considered  as  polynomials  in  p.  with  the  coefficients  written 

»  » 

as  polynomials  in  £' .  For  any  £g  with  Re£,  >  £  the  function 
F1(pi)  •  F(pif£g!  is  a  positive  function  m  the  variable  p. 
only.  Furthermore,  since  due  to  Theorem  3.3.13,  A  £q  *  3. 
B(£q)  *  0  the  degrees  of  d  and  g  in  p,  remain  unaltered  due 

to  the  substitution  £’  •  .  Hence,  F.  p, 1  is  a  positive 

function  having  a  pole  at  infinity,  necessarily  simple,  with 
residue  K^(£q).  Therefore,  is  real  and  positive,  thus 

in  particular  finite  and  therefore  holomorphic  for  all  Sec'  ' 
0.  This  implies,  in  view  of  a  standard  result  m  tke  theory 
cf  functions  of  complex  variacles  that  <  .3  a  ::~3ta~' 
same  holds  for  t.ne  rema.mng  ..  .  ?  .  a  .  t  :  the  :  >*•  f.-o  t  *  r. 


constants  for  i 


1  to  k . 


Consider  next  the  rational  function  Fq  -  ^/g  defined 

k 

via  Fq  -  F-  I  ITp^.  Clearly,  ReFQ(jw)  -  ReF(jw)  where  w  is 

any  real  k-tuple  such  that  g(j«)  a  0.  Furthermore,  deg^Q  < 
deg^  for  i  •  1  to  k .  This  follows,  e.g.,  for  the  case  i  -  1 
by  writing  the  polynomials  d  and  g  appearing  in  F  -  d/g  in 
the  same  way  as  above  and  taking  into  account  that  - 
A(£')/B(£').  Invoking  Theorem  3. 8.1, it  then  follows  that  Fq 
is  a  positive  function. 


3.9.  Conclusions 

The  artifice  of  sequentially  almost  complete  and  sequentially 
infinite  sets,  which  proved  to  be  very  useful  in  the  present 
context,  have  been  introduced.  Properties  of  widest-sense 
Hurwitz  polynomials,  self-paraconjugate  Hurwitz  polynomials, 
strict  sense  Hurwitz  polynomials  and  scattering  Hurwitz 
polynomials  have  been  studied.  Several  properties  of 
multivariable  positive  rational  functions  have  been 
investigated  in  this  context.  Reactance  Hurwitz  polynomials 
and  immittance  Hurwitz  polynomials  have  been  introduced. 
They  fall  out  as  the  appropriate  polynomials  occurring  as  the 
numerators  and  denominators  of  (rational)  reactance  functions 
and  positive  functions  respectively.  The  hierarchical 
relationship  between  the  several  classes  of  multivariable 
Hurwitz  polynomial  thus  delineated  is  diagramatically  shown 
in  Pigure  3.1,  in  which  an  arrow  (single  or  double)  points  to 
subclasses  of  polynomials,  whereas  double  arrows  originate 
from  classes  formed  by  products  of  elements  of  classes  to 
which  they  point.  A  nontrivial  result,  which  proves  to  be 
very  useful  in  theoretical  tests  for  the  property  of 

positivity  of  holomorphic  functions  and  is  formulated  in 
terns  of  the  behavior  of  its  real  part  on  the  distinguished 
boundary  of  the  domain  of  holomorphy  has  been  derived, 
finally,  in  view  of  its  validity  in  the  one-variable  case,  it 
s««as  plausible  to  conjecture  that  given  any  positive 
function  with  a  self-paraconjugate  Hurwitz  factor  in  its 
denominator  it  is  always  possible  to  extract  a  reactance  from 
»t,  thus  leaving  a  positive  function  with  scattering  Hurwitz 
lenom*n«toc  only.  A  partial  result  in  this  direction  has  been 
r  .  ided  . 
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CHAPTER  4 


REALIZATION  OF  STRUCTURALLY  PASSIVE 
MULTIDIMENSIONAL  DIGITAL  FILTERS 


4.1. Introduction: 

Various  synthesis  schemes  such  as  the  Darlington  synthesis 
scheme  for  synthesizing  lossless  transfer  functions  as  a 
cascade  interconnection  of  most  elementary  lossless  building 
blocks  such  as  inductors,  capacitors,  gyrators  etc.  in  the 
continuous  time  domain  have  now  become  classical  in  the 
network  theoretic  literature.  The  corresponding  problem  in 
the  discrete  time  domain,  namely  that  of  synthesizing  a 
discrete  lossless  bounded  (or  positive)  transfer  function  as 
a  structurally  passive  interconnection  of  elementary  lossless 
building  blocks  was  first  resolved  via  transformation  from 
prototype  problems  in  the  continuous  time  domain,  and  the 
resulting  class  of  filter  structures  are  now  known  as  the 
wave  digital  filters  (1).  Recently,  however,  successful 
attempts  to  derive  these  and  similar  other  discrete  domain 
results  without  making  explicit  use  of  tools  of  classical 
network  theory  have  been  made.  Notable  among  these  are  the 
orthogonal  filters  [2],  and  the  class  of  filters  described  in 
[3],  [4]  and  in  related  other  publications. 


In  view  of  interest  in  the  synthesis  of  multidimensional 
(k-D)  structurally  passive  digital  filters,  the  problem  of 
synthesis  of  k-D  lossless  two-port  transfer  scattering  matrix 
via  the  bisection  of  a  prescribed  two-port  into  a  cascade 
connection  of  two  lossless  two-port  sections  of  smaller 
"degree"  has  been  addressed  in  the  continuous  time  domain  in 
15].  An  attempt  to  develop  a  self  consistent  theory  for  the 
synthesis  of  k-D  structurally  passive  digital  filters 
independent  of  the  continuous  time  methods  have  alceady  been 
initiated  in  (6]  by  discussing  the  discrete  domain  stability 
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properties  of  a  class  of  multidimensional  polynoaials.  The 
present  report  addresses  the  problea  of  synthesizing  a  k-D 
discrete  lossless  bounded  matrix  as  the  transfer  function  of 
a  structurally  passive  two-port  digital  filter  directly  in 
the  discrete  domain.  Our  approach  is  to  bisect  the  prescribed 
discrete  lossless  two-port  into  a  cascade  interconnection  of 
two  discrete  lossless  two-ports  as  shown  in  figure  4.1. 
Necessary  and  sufficient  conditions  as  to  the  feasibility  of 
the  bisection  is  obtained.  It  falls  out  that  in  the 
one-dimensional  (1-0)  case  the  aforementioned  bisection  is 
always  feasible.  Our  discussion  in  the  1-0  context  thus 
yields  yet  another  algorithm  for  the  structurally  passive 
synthesis  of  1-D  lossless  digital  filter  transfer  functions, 
previously  not  discussed  in  the  literature. 
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4.2.  Notation ,  Terminology  and  Problem  forauation : 


Ha  firat  axplain  tha  notation  to  ba  usad  in  the  rest  of  the 
paper  in  tha  following.  Notations  such  as  a,  b,  c  will  denote 
polynoaials:  a-a(z),  b-b(z),  c-c(z)  in  k-variables  z  - 
(*1'*2*  •  Notations  such  as  n^  or  deg^a  will  denote  tha 
partial  degree  of  a  in  the  variable  z^.  The  coapact 
notation: 

n  .  n  .  n  «  n  . 
z“a  *  zx  aA*2  *•••*)<  **  will  also  be  used. 

A  A  **1  . ” ^  ^  — 1  .  •  n 

Finally,  *«a  (zx  ,z2  , ...z^  ),  A  ■  a.z”a, 

where  *  denotes  coaplex  conjugation.  Corresponding  notations 
for  various  polynoaials  other  than  the  polynoaial  a  will  also 
be  used. 

A  k-D  discrete  lossless  two-port  is  characterized  [6]  by  an 
associated  transfer  function  aatrix  H  as  in  (4.1)  or  by  a 
transaission  matrix  T  as  in  (4.2). 


[Hiu-b/a,  (Hi12-d2z  a/a 

( 4 . la , b ) 

n. 

(h]21-c/a,  (hi22— dBz  a/a 

( 4.1c,d) 

(Tln"da/c,  (T]12-b/c, 

( 4.2a,b) 

n, 

!T]2X-dBz  a/c,(TJ22-a/c 

( 4 . 2c ,d ) 

where  a,b,c  are  polynomials  such  that  4  is  scattering  Schur 
(6],  deg^b<deg^a,  deg^c£deg^a  for  all  i-1  to  k,  d  is  a 
unimodular  complex  constant  i.e.,  |d|  -  1  and 

a3  -  bB  +  c5 ,  (4.3) 

Note  that  (4.1)  can  be  regarded  as  a  discrete  k-D 
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counterpart  of  Belevltch  canonical  torn  for  tha 
representation  of  lossless  bounded  two-port  scattering 
matrices,  well  known  in  classical  network  theory. 

In  sore  specific  terns  the  problen  dealt  with  in  the  present 
report  can  be  described  as  follows.  Given  T  as  in  (4.2),  two 
uninodular  coaplex  constants  d' ,  d"  with  d  -  d'd",  and  the 
polynonial  factorization  c  •  e'e",  along  with  two  sets  of 
integers  n' -( n£ ,n£ , . . .n£ )  and  n"-( n£ , n£ , . . .n£ )  such  that 
dag^c'^nj,  deg^c^nj  and  nai**n[*nj[  for  *11  !•  to  k,  eeek  a 
factorisation  T  •  T'T",  where  T'  and  T"  are  also  discrete 

lossless  two-port  transnission  natrices  with  associated 
polynonials  (a',b',c')  and  { a" ,b" , c" )  respectively.  In 
addition,  the  requirements  deg^a'£n^  and  deg^a"<n^  needs  to 
be  satisfied.  Thus,  both  T'  and  Ta  are  also  required  to  have 
representations  siailar  to  those  expressed  in  (4.2).  In 
particular,  the  polynonial  triples  (a',b',c')  and  ( a" ,b" , c" ) 
are  also  required  to  satisfy  the  condition  that  A',  A"  are 
scattering  Schur,  deg^b ' <deg^a ' ,  deg^b" <degia"  for  all  i-L  to 
k  and  (4.41)  holds  true.  The  discrete  lossless  two-ports  with 
associated  transnission  natrices  T'  and  T"  resulting  fron  the 
factorization  of  the  transnission  natrix  T  is  shown  in  figure 
4.1. 

a ' A '  -  b'B'ec'8*  (4.4a) 

a"I"  -  b-B"  ♦  c"S"  (4.4b) 

It  then  easily  follows  by  considering  representations  of  T' 
and  T"  such  as  that  expressed  in  (4.2)  for  T  that  the 
condition  T-T'T"  is  equivalent  to  the  conditions  expressed  in 
(4.5a)  and  (4.5b)  in  the  following. 

*  -  na- 

a  -  a ' a  ♦  d*  b'b"z  a  (4.5a) 
b  •  d'a'b"  *  b'A"  (4.5b) 

The  above  considerations  motivate  the  following  definition. 
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Definition  4.2.1:  The  pair  of  polynomial  two-tuples  (a',b') 


and  (a",b")  la  said  to  ba  a  solution  to  ths  algsbraic 
aquation  if  aquations  (4.4)  and  (4.5)  along  with  tha  degree 
restrictions  deg^'^n?  and  deg^a"$n?  for  i«l  to  k  are 
satisfied. 

Me  note  that  in  tha  above  definition  tha  degree  restrictions 
on  the  polynomials  a'  and  a"  are  expressed  as  weak 
Inequalities  rather  than  equalities  as  is  required  by  the 
solution  to  the  original  problem.  Also,  the  restrictions  that 
the  polynomials  A*  and  A"  be  scattering  Schur  polynomials  are 
not  imposed  at  all. 

Definition  4.2.2:  A  polynomial  triple  (a",b",b'l  is  said  to 
satisfy  the  fundamental  equation  if  (4.6)  along  with  (4.7) 
holds  true. 

d'.ab"  -  ba"  -  -b' c"c"z-"  (4.6) 
deg ia"Sni  and  (4.7) 

Note  that  equation  (4.6)  is  obtained  by  eliminating  the 
polynomial  a'  from  (4.5a,b)  and  (4.4b).  Obviously  then  any 
solution  of  the  algebraic  equation  also  satisfies  the 
fundamental  equation. 


\ 
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4.3.  Solution  to  the  Algebraic  aquation: 

Clearly,  any  solution  to  the  problea  of  factorization  of  T 
into  T' T"  is  also  a  solution  to  the  algebraic  equation.  The 
following  theorem  shows  that  the  scattering  Schur  properties 
of  S'  and  S"  and  the  degree  requirements  on  a'  and  a"  are 
automatically  satisfied  by  any  solution  to  the  fundamental 
equation,  and  therefore,  any  solution  to  the  algebraic 
equation  is  also  a  solution  to  the  problem  of  factorization 
of  T  into  T'T". 

Theorem  4.3.1:  If  the  pair  of  polynomial  two-tuples  (a',b') 
and  (aN,b")  constitute  a  solution  to  the  algebraic  equation 
then  the  polynomials  4'  and  S"  are  scattering  Schur  and 
deg^a'-nj,  deg^-nj  for  all  i  ■  1  to  k. 

Proof :  Consider  the  rational  function  defined  as: 

e-(S'S-)/a-( (S'a-)/S]zE  (4.8) 

where  E"<Pi'P2' • • •?]<) * 

and  Pi-nai-(na,i+naBi)  (4.9) 

Since  S  is  a  scattering  Schur  polynomial,  n(+n?-n  . ,  and 
factors  of  a  scattering  Schur  polynomial  are  also  scattering 
Schur,  the  denominator  polynomial  of  ♦  is  also  scattering 
Schur . 

Furthermore,  straightforward  algebraic  manipulation  of 
equations  (4.4b)  and  (4.5)  yield  the  following. 

♦•( a"/c" ) ( 8"/2" ) ( 1-d' { B/I ) ( b"/a" ) )  (4.10) 

Since  it  follows  from  (4.3)  that  |B/5|<1  and  |b"/a"|<l  for 
I  2  ^ | * 1  for  i-I  to  k,  an  examination  of  '4.10)  yields  that 
Rei»>0  for  I  z  i  I  - 1 ,  wherever  *  is  well  defined.  Thus,  by 


invoking  a  result  proved  in  (6)  it  follows  that  ♦  is  « 
discrete  positive  function.  Consequently,  the  numerator 
polynoaial  of  ♦  ,  in  irreducible  fora,  is  a  widest  sense  Schur 
polynoaial.  This,  however,  iaplies  that  nai“na*  for  a11 

i-1  to  k.  The  last  equality  along  with  the  facts  that 

ni*na'i'  nI-na"i  and  nai"nl*ni  t09®ther  i®ply  that  n?«na,  and 
ni"na"i * 

The  widest  sense  Schur  property  of  A'  has  already  been 

established.  Next,  if  for  soae  Zq  on  the  distinguished 
boundary  of  the  polydisc  l*il$l,  i-1  to  k  we  have  a'(zQ)-0 
then  froa  (4.4a)  it  follows  that  b'(Zg)-0,  which  in  turn  due 
to  (4.5a)  iaply  that  a(Zg)-0.  Consequently,  if  A',  and  thus 
a,  had  a  sequentially  alaost  coaplete  set  [6]  of  zeros  on 
the  distinguished  boundary  then  a  would  also  have  a 

sequentially  alaost  coaplete  set  of  zeros  there,  which  is 

iapossible  if  A  scattering  Schur.  Therefore,  A  cannot  have 
sequentially  alaost  coaplete  set  of  zeros  on  the 

distinguished  boundary.  The  scattering  Schur  property  of  A  is 
thus  established  in  view  of  results  in  [5].  Similar 

arguments  hold  for  A". 

A  lossless  two-port  is  said  to  be  an  allpass  if  the 

polynoaial  b  associated  with  it  is  identically  equal  to  zero. 

We  will  need  the  following  result  as  a  preparation  for  the 
rest  of  the  discussions  to  follow. 

Theorem  4.3.2:  Any  discrete  lossless  two-port  transmission 

matrix  T  can  be  factored  as  T-T.T  T  ,  where  T-,  T  ,  T  are 

r  o  r  ror 

also  discrete  lossless  two-port  transmission  matrices  such 
that  Tf  and  Tf  are  allpass  and  if  TQ  has  representation  in 
terms  of  polynomials  a,  b,  c  as  in  (4.2)  £hen  the  polynomial 
a  is  relatively  prime  with  b  as  well  as  Br-a. 


extraction  of  discrete  lossless  two-port  sections  fro*  the 
front  and  rear  end  of  the  prescribed  transmission  matrix. 
Thus,  without  loss  of  generality  it  will  be  assumed  in  all 
forthcomming  discussions  that  the  polynomial  a  is  relatively 


4.4.  Solution  to  the  fundamental  equation ; 

We  will  need  the  following  lemmas. 

Lemma  4.4.1:  If  the  polynomial  a  is  relatively  prime  with  b 

Ha 

at  well  as  Bz  then  neither  a  nor  &  can  have  a  factor  in 

“  « a 

common  with  the  polynomial  z-  c"5". 

Proof :  Since  A  is  scattering  Schur  &  cannot  have  z^  as  a 

factor  for  any  i.  Also  by  rewritting  equation  (4.3)  along 
with  c-c'c*  in  the  form  of  (4.11) 

aS-b(Bz"a)  «•  (z-'c'C'  )  <c"5"z-"  )  (4.11) 

it  can  be  seen  that  if  A  or  a  had  a  factor  in  common  with  the 

n  " 

polynomial  (c"C"z-  ^  then  it  would  also  have  a  factor  ^ 
common  with  the  b(Bz~a)  i.e.,  in  common  with  either  b  o-  Bz~ 
both  of  which  is  ^uled  out  by  the  fact  that  a  is  relatively 
prime  with  b  and  Bz_a. 

Lemma  4.4.2:  If  the  polynomial  triple  (  a" , 0" , 0' ) is  a  solution 
to  the  fundamental  equation  then  deg^0"<n?.  Furthermore, 
there  exists  a  polynomial  a'  given  by  (4.12)  such  that  the 

M  N  M  N 

n  n 

polynomial  triple  ( 0"z  ,a"z~  ,  -a'd')  is  also  a  solution  to 

the  fundamental  equation.  Also,  we  have  that  deg^o'^n'  for 
all  i«l  to  k. 

•  '-(  5"c'B'*z"-’Vb)/a  (4.12) 

Proof :  The  fact  that  deg^O’Xn^  follows  directly  from  the 

fundamental  equatio-  for  the  triple  •«".0".0*'.  Next,  by 


c  <o 


(4.13) 


n"  *  n"  n  . 

*_  (aa"-d'bB" )-( z~  c"5 ) ( a"2' c' +  z  —  B'b)z  a/a 

Since  the  left  hand  side  of  (4.13)  is  a  polynomial,  due  to 
lemma  4.4.1,  S  must  divide  ( a"6 ' c ' +z-—  B'b)z_a.  Thus,  a'  in 
(4.12)  is  a  polynomial.  The  fact  that  deg^a'<n?  then  follows 
by  considering  the  degree  restrictions  on  c',  c",  a",  b,  and 
0'  and  a. 

Lemma  4.4.3;  the  polynomial  a  is  relatively  prime  with  b 
as  well  as  5z  a,  and  (a£,0£,3£)  and  are  two 
polynomial  triples  satisfying  the  fundamental  equation  then 
the  rational  function  given  in  (4.14)  is  a  constant. 

(a£B5-(3"a^)/(Z^"c"2")  (4.14) 

Proof :  By  multiplying  the  fundamental  equations  for 
(a£,0£,0£)  and  (aj, 0^*02 5  respectively  by  and  (-a£)  and 
adding  the  resulting  equations  one  obtains  equation  (4.15). 

dMaJ0^-a20J)-(0[a5-aj[0^)(z2',c"2")/a  (4.15) 

Since  the  lefthand  side  of  (4.15)  is  a  polynomial,  by 
invoking  lemma  4.4.1  it  then  follows  that  a  must  divide  the 
polynomial  P-(  0.' af-aTSi  )  .  Since  deg  .  P<n ’  +n"-n  -deg.  a  for  all 
i-l  to  k  we  have  that  P/a  is  a  constant.  The  result  then 
follows  by  noting  that  the  expression  in  (4.14),  in  view  of 
(4.15),  is  equal  to  (Pd'/a). 

Lemma  4.4.4:  the  polynomial  a  is  relatively  prime  wit.n  fc 
as  well  as  Bz  a,  and  (a",  8",  S’)  is  a  polynomial  triple 
satisfying  the  fundamental  equation  then  the  expression  giver, 
in  (4.16)  is  a  constant. 

5"  :  - 


Proof :  Follows  from  lemma  4.4.2  and  lemma  4.4.3. 

Lemma  4.4.5:  If  the  polynomial  triple  (aM,3",0')  is  a 

solution  to  the  fundamental  equation  then  there  exists  an  a' 

n  "  n  " 

as  given  by  lemma  4.4.2  such  that  (pa"+qz-  B" ,  pB"+qz-  a", 
pB'-qda')  is  also  a  solution  to  the  fundamental  equation, 
where  p  and  q  are  arbitrary  complex  numbers. 


Proof :  Obviously  follows  from  lemma  4.4.2. 
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4.5.  Factorization  of  the  discrete  lossless  two-port 
transmission  matrix: 

Two  polynomial  triples  (a£,3£,3£)  and  each 

satisfying  the  fundamenatal  equation  will  be  said  to  be 
linearly  dependent  if  there  exists  constants  p  and  q  not 
simultaneously  zero  such  that  pa£+qot2»pBJ+  q82*p8£+q32*0 * 

Also,  a  solution  (  a" ,  fi" ,  (S' }  to  fundamental  equation  will  be 
said  to  be  nonsingular  if  a"a"^8"B" . 

The  following  two  theorems  constitute  the  major  results  of 
this  report. 

Theorem  4.5.1:  Assuming  £hat  the  polynomial  a  is  relatively 
prime  with  b  as  well  as  Bz“a,  the  problem  of  factorization  of 
discrete  lossless  two-port  transmission  matrix  T  admits  a 
solution  if  and  only  if  there  exists  a  nonsingular  solution 
(a",  8"  ,8')  to  the  fundamental  equation. 

Proof :  Necessity  is  obvious.  If  (a", 3", 6')  is  a  nonsingular 

solution  to  the  fundamental  equation  then  due  to  lemma  4.4.5, 
n"  n" 

a"-p«"+qz-  3",  b"»p8"+qz-  a",  b'«p3'-qda'  is  a  solution  to 

the  fundamental  equation.  straightforward  algebraic 
manipulation  then  yields  that 

(a"S"-b"B")/c"S"-( | p | 2- 1 q I  2 ) K  (4.17) 
where  K-( a"5"-3"B" )/c"5"  (4.18) 

Since  due  to  lemma  4.4.3  and  nonsingularity  of  (a",8",8'),  K 
is  a  nonzero  constant,  by  proper  choice  of  p  and  q  in  the 
right  hand  side  of  (4.17)  it  is  possible  to  have 
(a"S"-b"B")-c"5". 


Furthermore,  there  exists  a'  such  that  ,  £"z-  ,  -a'd), 

by  virtue  of  lemma  4.4.2,  satisfies  the  fundamental 


equation.  It  can  then  be  verified  via  routine  algebraic 
manipulation  that  the  pair  of  two-tuples  (a',b')  and  (aN,b”) 
satifies  the  algebraic  equation,  and  thus,  due  to  theorem 
4.3.1,  is  a  solution  to  the  problem  of  factorization  of  T. 

Theorem  4.5.2:  Assuming  that  the  polynomial  a  is  relatively 
prime  with  b  as  well  as  with  Bz“a,  the  problem  of 
factorization  of  discrete  lossless  two-port  transmission 
matrix  T  admits  a  solution  if  and  only  if  there  exists  two 
linearly  independent  polynomial  triples  (a?,0?,0f),  i-1,2 
each  of  which  satisfy  the  fundamental  equation. 

Proof:  Necessity  is  obvious.  If  one  of  the  solutions 

i-1,2  is  nonsingular  then  sufficiency  follows 
from  theorem  4.5.1.  If  both  solutions  are  singular  then  the 
triple  (a",b",b')  obtained  as:  a"»pa£+q0&,  b"-p0£+q02j, 

b'-p0£+q0£,  where  p  and  q  are  complex  numbers,  satisfies  the 
fundamental  equation.  Algebraic  manipulation  then  yields  that 

(a"3"-b"B" )/c"2"  •  L  ♦  L  (4.19) 

L  -  p*q(«5«£-05B*)/c"e"  (4.20) 

By  invoking  lemmas  4.4.2  and  4.4.3  it  then  follows  that  L  in 
(4.20)  is  a  constant,  and  thus,  t-L* .  Furthermore,  by 
following  arguments  similar  to  those  in  [5]  it  can  be  proved 

via  the  use  of  results  in  (7]  that  L*0  if  p#0  and  q#0.  (The 

details  of  this  derivation  is  left  out  of  here  for  the  sake 

of  brevity).  Consequently,  by  proper  choice  of  p  and  q  in 

(4.19)  and  (4.20)  it  is  possible  to  have  a"3"-b"B"-c"5" .  The 
rest  of  the  proof  follows  by  imitating  the  last  paragraph  in 
the  proof  of  theorem  4.5.1. 

The  fundamental  equation  (4.6),  when  considered  as  a  set  of 
linear  simultaneous  equations  involving  the  coefficients  of 
the  polynomials  a",  b",  b',  along  with  the  upper  bounds  on 

their  degrees,  turns  out  to  be  ove rdete trained  in  general. 


More  explicitly,  we  note  that  the  unknown  polynomials  a",  b" 
and  b'  contain  a  total  of  u  unknown  coefficients,  whereas  the 
total  number  of  linear  simultaneous  equations  can  easily  be 
found  to  be  equal  to  e,  u  and  e  being  as  given  in  (4.21)  and 
(4.22)  below. 

k  k 

u  -  2  H  (n?+l)  ♦  n  (n.'+l)  (4.21) 

i-1  1  i-1  1 

k 

e  -  B  (2n?  ♦  nf  +  1)  (4.22) 

i-1  1  1 

Since  for  k>l  we  have  e>u  in  a  generic  situation  a  solution 
to  the  problem  may  not  exist. 

Furthermore,  in  order  for  the  digital  network  so  synthesized 
to  be  'computable'  it  may  not  contain  delay  free  loops 
arising  from  cascading  of  two  elementary  sections.  In  spite 
of  the  fact  that  it  is  known  (8]  that  this  problem  can  always 
be  circumvented,  at  least  in  the  one-dimensional  case,  by 
incorporating  digital  equivalents  of  unit  elements  it  is  of 
interest  to  note  that  by  properly  utilizing  the  flexibility 
in  the  choice  of  p  and  q  in  (4.20)  it  is  always  possible  to 
avoid  the  occurrence  of  such  delay  free  loops  in  the  filter 
structure.  This  point  is  further  elaborated  in  the  following 
section. 


4.6.  One-dimensional  synthesis  as  a  special  case: 

In  the  one-dimensional  case  i.e.,  if  k  -  1,  a  closer 
examination  of  (4.21)  and  (4.22)  reveals  that  we  have  u-e-2, 
and,  therefore,  there  are  two  more  unknown  coefficients  than 
the  number  of  linear  equations  in  the  set  of  linear 
simultaneous  equations  which  determine  the  solution  to  the 
fundamental  equation.  Thus,  there  are  at  least  two  linearly 
independent  solutions  of  the  fundamental  equation,  and  in 
view  of  theorem  4.5.2,  the  problem  of  factorization  of  T 
always  admits  of  a  solution.  Consequently,  structurally 
passive  synthesis  for  T  is  achieved  by  performing  a  sequence 
of  further  factorizations  of  T'  and  T"  into  discrete  lossless 
transmission  matrices  of  progressively  lower  complexity, 
until  a  stage  is  reached  when  each  of  the  resulting 
transmission  matrices  cannot  be  factorized  any  further.  This 
latter  situation  corresponds  to  the  case  that  each  of  the 
two-ports  resulting  from  the  decomposition  satisfy  deg  a  -  1, 
deg  c  i  1  and  deg  b  £  1.  However,  if  the  specified  two-port 
transmission  matrix  T  has  real  coeif f icients  for  its 
numerator  and  denominator  polynomials  and  realization 
involving  only  real  multipliers  are  sought  then  the 
constituent  two-ports  may  also  be  of  the  type  deg  a  -  2,  deg 
c  ■  2,  and  deg  b  <  2.  Two-port  sections  of  the  above  types 
will  be  called  elementary  sections  and  can  in  turn  be 
realized  in  structures  possibly  other  than  the  cascade 
structure  by  exploiting  synthesis  techniques  as  discussed, 
for  example,  in  [4]. 

To  address  the  issue  of  absence  of  delay  free  loops  at  the 
junction  of  the  two-ports  associated  with  T'  and  T"  it  may 
be  noted  that  the  only  restriction  governing  the  choice  of 
the  numbers  p  and  q  is  that  the  right  hand  side  of  (4.17)  or 
(4.19).  be  equal  to  one.  This  flexibility  in  the  choice  of  p 
and  q  may  thus  be  exploited  to  make  b"(0)  -  0,  which  ensures 


the  absence  of  the  delay  free  loops  of  the  type  mentioned 
above.  Furthermore,  since  A  is  Sehur  and  thus  &(0)*0,  it 
follows  from  (4.5b)  that  if  b(0)>0  and  b"(0)«0  then  b'(0)«0. 
This  fact  gaurantees  that  the  original  two-port  can  be 
decomposed  into  cascade  interconnection  of  elementary 
two-ports  in  such  a  way  that  the  b-polynomial  associated  with 
each  of  the  constituent  two-ports,  except  possibly  the  one  at 
the  extreme  left,  is  equal  to  zero  for  z«0.  Absence  of  delay 
free  loops  from  each  junction  is  thus  gauranteed. 

Realizations  for  elementary  sections  with  deg  a-1,  deg  c  <  1, 
deg  b  $  1  and  b(0)>0  as  interconnections  of  Gray-narkel 
sections  and  delays  are  shown  in  figures  4.4  and  4.5. 
Gray-Markel  sections  of  two  different  kinds  used  in  these 
figures  are  shown  in  figures  4.2  and  4.3.  An  elementary 
section  with  deg  c  ■  2,  deg  a  ■  2,  deg  b  $  2  and  b(0)-0  is 
shown  in  figure  4.6.  Thus,  an  arbitrary  lossless  two-port  can 
indeed  be  synthesized  as  a  cascade  interconnection  of  these 
elementary  sections  only.  It  turns  out  that  elementary 
sections  just  referred  to  are  exactly  the  same  as  those 
discussed  in  the  literature  [11,(2], [3]. 
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4.7  Conclusions: 


A  simple  algorithm  involving  the  examination  of  rank  of  a  set 
of  linear  simultaneous  eqations  for  studying  the 
synthesizability  of  an  arbitrary  multidimensional  lossless 
two-port  in  a  cascade  structure  has  been  derived  via 
factorization  of  the  associated  transmission  matrix  T.  It 
turns  out  that  under  a  generic  situation  synthesis  in  a 
cascade  structure  may  not  be  feasible.  In  the  special  case  of 
one-dimension  the  algorithm  provides  a  new  method  of 
realizing  structurally  passive  filters  directly  in  the 
digital  domain.  The  problem  of  multidimensional  synthesis, 
though  not  necessarily  in  cascade  structure,  can  also  be 
addressed  via  factorizations  of  hybrid  matrix  or  the  tranfer 
function  matrix  associated  with  the  lossless  two-port.  The 
class  of  multidimensional  lossless  two-ports  thus 
synthesizable  along  with  the  class  identified  in  the 
present  study  would  thus  broaden  the  whole  class  of 
synthesizable  multidimensional  structurally  passive  lossless 
two-ports.  It  may  be  remarked  that  even  though  in  the  k>l  case 
synthesis  may  not  be  feasible  for  an  arbitrary  discrete 
lossless  T,  the  possibility  of  synthesis  for  special  classes 
of  discrete  lossless  T  is  by  no  means  ruled  out.  This  is 
especially  true  in  view  of  synthesizability  of  certain 
classes  [9],  [10]  of  two-dimensional  continuous  time  systems 
arising  in  studies  of  lumped-distributed  netwoks.  The  class 
of  multidimensional  discrete  lossless  two-port  transmission 
matrices  T,  which  admits  of  such  synthesis  remains,  however, 
to  be  identified. 
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CHAPTER  5 


MULTIDIMENSIONAL  INTERPOLATION  AND  DECIMATION 
SCHEMES  FOR  FAST  IMPLEMEMNTATION 


5.1.  Introduction 


The  processes  of  sampling  rate  increase  and  sampling  rate 
reduction  commonly  referred  to  as  interpolation  and 
decimation,  are  required  whenever  it  is  necessary  to 
change  from  one  sampling  rate  to  another.  The  fact  that 
many  commonly  encountered  one-dimensional  signal 
processing  tasks  such  as  speech  processing,  single  side 
band  frequency  multiplexing  require  the  processes  of 
sampling  rate  change  is  now  well  known  [3].  Similarly  a 
large  number  of  multidimensional  signal  processing  [5] 
tasks  also  require  the  operation  of  digital  interpolation 
and  decimation.  Such  application  areas  include 
transmission  of  television  pictures  [13],  antenna 
beamforming  [5],  target  tracking  [1],  astronomical  data 
processing  [17]  geophysical  signal  processing  [16],  and 
medical  tomography  [9].  We  point  out  exactly  how  the 
specific  problem  of  multidimensional  sampling  rate 
alteration  enters  into  some  of  these  applications.  In 
radio  astronomical  observations  it  is  often  desirable  to 
estimate  the  radio  brightness  of  the  sky  at  intermediate 
points  from  observations  made  by  directing  the  measuring 
antennas  at  a  regular  array  of  points  in  the  sky  [17];  in 
the  problem  of  transmission  of  television  pictures 
efficient  coding  schemes  require  that  the  time-varying 
image  signal  be  known  between  two  successive  picture 
frames  [131;  whereas  in  X-ray  computed  tomography  the 
problem  of  interpolation  manifests  itself  when  a  higher 
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resolution  in  reconstruction  is  necessitated  [9].  Similar 
other  examples  of  need  of  multidimensional  interpolation 
as  well  as  decimation  schemes  can  be  drawn  from  these  and 
other  application  areas  mentioned  earlier.  A  feature 
common  to  all  of  these  multidimensional  processing  tasks 
is  the  enormous  amount  of  computational  requirement.  This 
fact  becomes  especially  prohibitive  to  practical 
implementation  if  real  time  or  adaptive  applications  are 
called  for. 

On  the  other  hand,  recent  advent  of  VLSI  technology  has  the 
potential  to  make  such  computation  intensive 
multidimensional  signal  processing  tasks  by  increasing  the 
throughput  rate  via  utilization  of  new  concepts  such  as 
parallelism,  pipelining,  concurrency  modularity  of 
implementation  etc.  [8].  The  more  recent  optical 
technologies  [2]  provide  yet  another  potential  means  for 
highspeed  implementation  of  many  multidimensional  signal 
processing  algorithms.  In  this  context,  the  need  for 
reconsidering  existing  signal  processing  algorithms  as 
well  as  that  of  designing  algorithms  for  previously 
intractable  problems  have  already  been  recognized  in 
general  [11],  and  both  new  computational  schemes  and  their 
implementations  in  hardware  for  solving  specific 
multidimensional  signal  processing  tasks  are  now  beginning 
to  emerge  [6]  [12].  It  is  in  this  perspective  that  the 
general  problem  of  designing  an  algorithm  for 
interpolation  and  decimation  of  a  broad  class  of 
multidimensional  signals  is  investigated  in  the  present 
report. 

For  the  type  of  applications  under  consideration,  it  is 
important  to  understand  the  processes  of  interpolation  and 
decimation  from  the  point  of  view  of  digital  signal 


processing  rather  than  from  numerical  analysis  standpoint. 
For  example,  linear  interpolation  is  not  satisfactory  in 
most  digital  signal  processing  tasks.  In  classical 
numerical  analysis,  the  inadequacies  of  polynomial 
interpolation  schemes  such  as  the  multidimensional 
Lagrange  type  interpolation  method  lead  to  the  use  of 
higher  order  polynomials,  the  inappropriateness  of  which 
has  been  pointed  out  in  [3]  in  the  one-dimensional 
context.  In  the  present  study  we  consider  the  problems  of 
interpolation  and  decimation  based  on  the  frequency  domain 
description  of  the  multidimensional  signal.  It  then  turns 
out  that,  as  has  already  been  discussed  in  the  literature 
for  one-dimensional  signals  [3],  the  problem  of 
interpolation  and  decimation  of  multidimensional  signals 
can  also  be  interpreted  as  linear  filtering  operations  in 
the  frequency  domain.  More  importantly,  the  filtering 
scheme  leads  to  computational  structures,  which  is  highly 
modular  and  derive  full  advantages  of  parallel  and 
pipeline  implementation.  Performance  of  the  computational 
algorithm  so  designed  is  also  examined  by  experimenting 
with  both  real  and  synthetic  two-dimensional  signals. 

In  section  5.2  the  notation,  terminology  and  the  sampling 
scheme  to  be  used  for  the  rest  of  the  report  is 
introduced.  The  problems  of  interpolation  and  decimation 
are  formulated,  and  the  filtering  schemes  leading  to  their 
solution  are  discussed  in  section  5.3.  The  fact  that  such 
computational  schemes  can  be  implemented  by  exploiting  the 
concepts  of  both  parallelism  and  pipelineability  forms  the 
contents  of  section  5.4.  Some  filter  design  examples  of 
interpolators  and  decimators  and  their  performances  on 
image  data  are  presented  in  section  5.5  and  in  section  5.6 
conclusions  are  drawn. 


5.2.  Multidimensional  Periodic  Samplin 


Of  the  several  ways  to  generalize  (5]  one-dimensional 
(1-0)  periodic  sampling  schemes  to  multidimensions  (N-D) 
(N£2),  the  most  straightforward,  although  not  the  most 
efficient  [15],  is  periodic  sampling  in  rectangular 
cartesian  coordinates,  which  we  will  simply  call 
rectangular  sampling. 

In  what  follows  underlined  characters  will  be  used  to 
denote  column  vectors  and  the  notation  will  be  used  to 
denote  the  transpose  of  a  vector.  For  example, 
T'«(T^,T-,...  ) ,  where  T  is  a  column  vector.  Similar 

notations  will  be  used  for  variables  such  as  t,  2,  «,  w, 
n,  k  etc.  the  later  physical  meanings,  of  which  will  be 
made  clear  when  the  context  arises. 

If  x  (t),  is  a  multidimensional  continuous  signal,  the 
discrete  signal  x(n)  obtained  from  it  by  rectangular 
sampling  is  given  by: 

x(n)-xa(niT1,n2T2*  *  *  *nNTN)  (5.2.1) 

where  T^,T2,...TN  are  positive  constants  known  as  the 
sampling  intervals  or  periods  in  the  respective  sampling 
directions . 

The  N-D  Fourier  transform  X_(2)  of  the  continuous  signal 

a  *■ 

x  (t),  and  its  inverse  are  given  in  (5.2.2a,b). 

-  -j(a.t')  2-  j(fl.t') 

Xa(2)-r  xa(t)e  dt;  xa(t)-(l/4a  )/  Xa(2)e  dt 

(5.2.2a,b) 

The  discrete  sequence  x(n)  obtained  by  sampling  x  (t)  at 

a  ™ 

spatial  locations  t0  N  can  then  be  obtained 


as  in  (5.2.3a),  whereas  the  Fourier  transform  of  x(n)  is 
given  by  X(w)  in  (5.2.3b). 


,  •  -j(2.t' )  j ( w.n' )  . 

x(n)»( l/4n^ ) J  Xa(2)e  ~  d2;  x(«)-Zx(n)e 

-  ( 5.2. 3a, b) 

We  assume  that  the  signal  x.(t)  is  bandlimited.  More 

O  * 

specifically,  we  assume  that  the  Fourier  transform  X_(2) 
have  a  support  which  is  contained  in  the  hypercube 
| 2^ | <W| <•  for  i-1  to  N  in  the  N-dimensional  frequency 
space  i.e.,  Xft(2)-0  for  (2^|>w^,  i-1  to  N.  It  then  follows 
from  the  multidimensional  version  of  the  Nyquist  sampling 
theorem  [5]  that  the  continuous  signal  x_(t)  can  be 
recovered  from  the  discretized  multidimensional  signal 
x(n)  according  as  equation  (5.2.4),  in  which  ta<r/Wa. 


N 

xm(t)-E  x(n)  H 
a  “  n  ~  i-1 


sin( r/T. ) ( t.-niTi ) )/( r/t. ) ( ti-n.Ti ) 1 

(5.2.4) 


Equations  (5.2.1)  and  (5.2.4)  taken  together,  form  the 
basis  of  the  multidimensional  sampling  theorem  in 
rectangular  cartesian  coordinates. 


104- 


5.3. 


Multidimensional  Sampling  Rate  Conversion 


The  process  of  sampling  rate  conversion  is  one  of 
converting  the  sequence  x(n)  obtained  from  the  sampling 
of  the  bandlimited  signal  x#(t)  with  periods  ,  to  another 
sequence  y(m)  obtained  from  sampling  xft(t)  with  periods  T£ 
for  each  i  -  1  to  N.  Conceptually  at  least,  the  most 
straightforward  way  to  perform  this  conversion  is  to 
reconstruct  x  (t)  (or  the  low-pass  filtered  version  of  it) 
from  the  samples  of  x(n)  and  then  resample  xa(t)  (assuming 
that  it  is  sufficiently  bandlimited  for  the  new  sampling 
rate)  with  periods  T[,  i-1  to  N  to  give  y(m).  For  any  m, 
the  value  of  y(m)  can  be  then  obtained  as: 


y(m)»xa(t)  for  t^«m^T£,  i“l  to  N  (5.3.1) 

By  substituting  (5.3.1)  into  (5.2.4)  and  renaming  the 
variables  n  as  k,  we  have: 

N 


y(m)  -  Z  x(k )  0  sin(x((miT^/Ti)-ki))/ii((miT^/Ti)-ki) 
^  1-1  (5.3.2) 


5.3.1  Interpolation  off  Multidimensional  Signals 

If  the  sampling  rate  in  the  i-th  dimension  is  increased  by 
an  integer  factor  ,  then  the  new  sampling  period  T£  for 
i>l  to  N  are  given  by 

T^/Ti-1/Li  (5.3.3) 

This  process  of  increasing  the  sampling  rate 
(interpolation)  of  a  signal  x(.)  by  implies  that  we 
must  interpolate  (L^-l)  new  sample  values  between  each 
pair  of  sample  values  of  x(.).  By  substituting  (5.3.3) 


into  ( 5.3.2) ,  (5. 3.4a)  is  obtained,  in  which  h^(s^)  for 

integer  values  of  s^,  i-i  to  N  are  given  in  (5.3.4b). 


y(a)- 


N 

E  x(k)  n  h. (m.-k.L. ) ; 
k  ”  i-1  1  1  1  1 


h^s^-sint  KSj/L^/t  KSj/L^) 

(5.3.4a,b) 


An  alternative  formulation  of  equation  (5.3.4a),  as  given 
in  (5.3.6),  can  be  obtained  via  the  introduction  the 
change  of  variables: 

ki"l®i/I»iJ  -n^,  for  i-1  to  N  (5.3.5) 


where  LuJ  denotes  the  integer  less  than  or  equal  to  u. 


y(m)- 


N 

£  x( Im/L-nJ )  n  h, (m.-La./L. J .L.+n.L. ) 

j  i  1  1  1  1 

n 


N 

E  n  h, ( n. L. +a. ®L, )x( In/L-nJ ) 
n  i-1  1  11  1  1 


(5.3.6) 


where  denotes  the  value  of  m^  modulo  for  i-  1  to 
N  and  the  notation  Lm/L-nJ  is  taken  to  mean  the  N- tuple  of 
integers  s-(  s^Sj, . .  .sN)  with  si*  J”ni »  i-1  to  N. 
Equation  (5.3.6)  expresses  the  output  y( . )  in  terms  of  the 
input  x(.)  and  the  set  of  one-dimensional  sequences  h^(.) 
i-1  to  N,  as  given  in  (5.3.4b).  Thus,  in  a  compact 
notation  y(a)  can  be  written  as  in  (5.3.7a),  where  gm(n) 
is  as  expressed  in  (5.3.7b)  for  all  N-tuple  of  integeTs  m 
and  n. 


y(m)-E  g 
n  ' 


,(n)x(  Im/L-nJ ) ; 


gaU)- 


N 

n 

i-1 


hi(niLi+Bi®Li } 


(5.3.7a,b) 


Note  that  ga(n)  is  periodic  in  m^  with  period  for  each 
i-1  to  N.  Furthermore,  by  referring  to  (5.3.7b)  it  follows 
that  9a(n)  is  a  product  separable  function,  and  each  of 
its  factors  are  periodic  in  nr  with  respective  periods 
i .  e . , 
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(5.3.8a,b) 


9.(n)-  J  9m  (Hi)  ?  ^(n^-h^n^+m^) 


where  (5.3.8b)  holds  for  each  i»l  to  N.  Thus,  by  using 
(5.3.8a),  (5.3.7a)  can  be  written  as  (5.3.9). 


y(m)“  I  g  (nN)...  I  g  (n, )x( Im/L-nJ )  (5.3.9) 

nN - "l -  1 

An  alternative  formulation  of  equation  (5.3.9),  which 
yields  the  implementation  of  the  equation  in  digital 
filtering  terms  is  given  in  equation  (5.3.10). 


y  ( # • • • * v . vM ) 


-  1  V(ni,y 

ni—  1 


(i-1) 


( ® ^ » •  •  •  ^ »  L®4 /l»i J ) 

(5.3.10) 


for  l£i£k,  where  y  (.)  -  x(.) 


Note^N^hat  (5.3.10)  together  with  (5.3.9)  yields  that  y( . ) 
-  y  (.).  Thus,  output  y(m),  as  given  in  (5.3.9),  can 
also  be  computed  via  ^  £l^e  recursive  construction  of  the 
intermediate  signals  y  (.),  for  1  <  i  <N.  Furthermore, 
for  each  i  in  1  $  i  $  N,  (5.3.10)  can  be  interpreted  as 
the  input-output  equation  of  a  one-dimensional  spatially 
varying  filter  operat^jg^ jOn  the  i-th  dimension  of  the 
intermediate  signal  y  (.),  whose  impulse  response  is 
periodically  space  varying  with  a  period  equal  to  .  A 
closer  examination  of  (5.3.10)  reveals  that  (5.3.10)  can 
be  digitally  implemented  in  N  stages  as  shown  in  Figure 
5.1,  where  the  i-th  stage  represents  a  set  of  one¬ 
dimensional  filters  o^rj^ing  on  the  i-th  dimension  of 
intermediate  signal  y  (.),  one  for  each  value  of  the 
indices  m^  >  •  •  .m^_^ *  +  i *  •  • • v|<  »  t0  produce  the  next 


intermediate  signal  y  ( . ) .  In  the  special  case,  when  N-2 
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i.«.,  foe  two-dimensional  image  signals,  for  example, 
(5.3.10)  can  be  written  as  follows. 

y(1)(m.,v2)-  I  g  (n.)x(  Lm./L.J-n.,-^)  (5.3.11) 

n^--»  1 

y(m1,m2)-y<2) (m1,m2) 

■  I  gB  (n2)y<3,)  (m2,  Lm2/L2J-n2)  (5.3.12) 

n2— •  2  “ 

Equation  (5.3.11)  represents  a  one-dimensional 
interpolator  which  interpolates  the  rows  of  the  given 
input  image  x(.).  This  row  interpolation  is  performed  on 
each  row  of  the  input  signal  x(.)  and  the  results  are 
stored  in  an  intermediate  image  signal  y^^(m^,v2).  The 
second  step  is  an  implementation  of  (5.3.12),  which 
amounts  to  performing  the  operation  of  interpolation  on 
the  columns  of  the  intermediate  image  y*^(. ).  The  entire 
process,  therefore,  can  be  implemented  in  two  stages  of 
1-0  interpolations.  It  must  be  noted  that  by  rewriting  in 
(5.3.9)  the  order  in  which  the  summations  over  different 
indices  are  considered,  it  is  also  possible  to  perform  the 
column  interpolations  first  and  the  row  interpolations 
next. 

5.3.2  Decimation  of  Multidimensional  Signals 

The  process  of  reducing  the  sampling  rate  (decimation)  of 
x(n)  by  an  integer  factor  in  the  i-th  dimension,  is 
considered  next.  If  TJ/T^-H^/1  for  i*  1  to  N  then  the  new 
sampling  rate  is  given  by  r£«l/Tr -l/m^Tj-rj/ft^ . 

In  order  to  lower  the  sampling  rate  and  to  avoid  aliasing 
at  this  lower  rate,  it  is  necessary  to  filter  the  signal 
x(n)  with  a  digital  low-pass  filter  whose  unit  impulse 


response  is  denoted  by  h(n) .  The  sampling  rate  reduction 
is  then  achieved  by  forming  the  sequence  y(m)  by  saving 
for  each  i  in  1  <  i  <  N  only  every  J^-th  sample  in  the 
i-th  dimension  of  the  filtered  output. 

If  h(.)  denotes  the  N-D  impulse  response  of  the  ideal 
low-pass  filter  then  it  follows  that  the  signal  at  the 
output  of  the  low-pass  filter  is  given  by  : 


w(n)-  I  h(k)x(n-k)  (5.3.13) 

k  ~ 

Furthermore,  if  y(m)  is  obtained  by  considering  every 
M^-th  sample  in  the  i-th  dimension  of  the  signal  w(n)  then 
y(m)  is  given  by  (5.3.14),  where  the  notation  x(H.m-k)  is 
taken  to  mean  xfn^m^-k^ ,N2m^  kj, . . .M^m^-k^) . 

y(m)-[w(n)l  „  »£  h(k)x(M.m-k)  (5.3.14) 

-  ni-miMi  k  -  - 

If  the  frequency  response  of  the  ideal  low-pass  filter 
having  impulse  response  h(n)  is  given  by  (5.3.15)  then 
h(n)  is  product  separable  (5)  and  can  be  written  as  in 
(5.3.16) . 

1  for  |ujJ<x/M^,  for  i«  1  to  N 
H(w)-<  (5.3.15) 

0,  otherwise 


N 

h(n)-  n  h^tn^);  h^fn^  )-sin(  xn^  )/*n^  for  i>  1  to  N 

l"1  (5. 3. 16a, b) 

Making  use  of  (5.3.16)  in  (5.3.14)  it  follows  that  the 
output  signal  from  the  entire  decimator  y(m)  can  be 
written  as  in  (5.3.17). 


y(m)' 


h'(kk) 


N 


£ 

k. 


h' ( k^ )x( M.m-k ) 


(5.3.17) 
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Furthermore,  as  for  the  interpolator  we  note  that  the 
output  signal  y  ( m)  can  be  computed  via  the  recursively 
defined  intermediate  signals  y^  ( . ) ,  1  <  i  <  N  as  given 
in  (5.3.18)  below,  where  y^(.)  -  x(.).  We  then  have 

y(m)«y(N)(m). 

y(i)<*i'.. **i/Vi+i#  * • •vn)“ 

I  hMk^y*1  1)(m1,...mi_1,Mimi-ki,vi+1,...vN)  (5.3.18) 

*i 

Clearly,  for  each  i  equation  (5.3.18)  represents  the 
operation  of  performing  1-D  decimation  in  the  i-th 
direction  of  the  intermediate  signal  y(i-l)(.),  as  a 
result  of  which  the  next  intermediate  signal  y^N.)  is 
obtained,  and  the  computation  of  y^(.)  in  (5.3.18)  for 
values  of  i-  1,2,...N  correspond  to  the  implementation  of 
N  stages  of  such  1-D  decimators  in  succession.  For 
two-dimensional  signals,  for  example,  i.e.,  if  N-2  the 
entire  operation  can  be  interpreted  as  first  decimating 
the  rows  of  the  discrete  signal  and  then  decimating  the 
columns  of  the  resulting  signal  obtained  from  the  output 
of  the  first  stage. 


5.4.  High  Speed  Implementation  of  Interpolators  and 
Decimators 

We  note  that  (5.3.10)  and  (5.3.18)  are  the  basic  equations 
for  implementation  for  N-D  interpolators  and  decimators, 
which  can  be  implemented  in  N  decoupled  stages.  Each  such 
stage,  in  fact,  consists  of  a  set  of  1-D  interpolation  or 
decimation  filters.  In  this  section  the  fast  computational 
scheme  associated  with  (5.3.10)  and  (5.3.16)  that  results 
in  this  highspeed  implementation  will  be  discussed.  It 
will  be  shown  that  the  set  of  1-D  interpolators  as  well  as 
the  decimators  just  referred  to  can  be  implemented  by 
using  certain  types  of  1-D  interpolators  or  decimators 
known  to  be  the  polyphase  filters  [3]  as  our  basic  module 
of  implementation.  Only  nonrecursive  implementations  of 
this  basic  module  will  be  sought  in  the  present  study. 

We  first  consider  the  i-th  stage  of  implementation  of  the 
interpolator.  Similar  considerations  apply  to  each  such 
stage.  The  impulse  responses  gB  ( n ^ )  of  the  1-D  filters 

mentioned  in  the  previous  paragraph,  by  virtue  of 

equations  (5.3.8b),  are  periodically  shift  varying  in  m^ 

with  a  period  L. .  Due  to  this,  for  any  fixed  set  of  values 

1  (  i ) 

of  m^, . .  +  the  computation  of  y  (.)  can  be 

carried  out  via  the  use  of  different  shift-invariant 

filters  in  such  a  way  that  each  filter  provides  every 

L^-th  sample  of  y^  ( . )  in  the  i-th  dimension. 

Consequently,  for  a  given  value  of  i  as  shown  in  Figure 

5.3,  the  entire  filtering  operation  represented  by 

(5.3.10)  can  be  implemented  as  a  parallel  interconnection 

of  different  shift  invariant  filters  having  y^*”^ (  . ) 

at  its  input,  followed  by  the  process  of  sampling  rate 

expansion  by  a  factor  of  i.e.,  insertion  of  (L^-l) 


zeros  between  two  consecutive  samples  in  the  direction  m^. 
Finally,  for  each  j-1  to  (L^-l)  the  output  from  j-th  such 
shift  invariant  filter  is  shifted  j  space  units,  and  the 
resulting  signals  are  added  to  obtain  y^(.). 


Similar  implementational  considerations  also  apply  to  each 
of  the  N  stages  of  the  decimator.  By  making  the 
substitution  ,ciBriNi4>pi  in  (5.3.18)  it  follows  after  some 
algebraic  manipulations  that  (5.4.1)  to  m  (5.4.3)  hold 
true. 


y(i)-(m1, . . • • *VN) 


Mj-1 

-  £ 

Pi 


,( i-1 ) 


1  n  r  1  .pe,(ri)ye,  <®i,...«i_1,*i-rl,vi+1,...vN) 
’i"°  ri*“"  1  .  1  (5.4.1) 

where  y^"1*  (mi# . .  •  •  V 

-  y(i  1)(mi,...mi_1,miMi-pi,vi+1,...vN)  (5.4.2) 


and 


( rj^J-h'  (miHi+pi) 


for  pj^-0, 1, . . .  (Mi-1 ) 


(5.4.3) 


Thus,  for  fixed  integer  values  of  m^,...m^_^,  *  •  •  •  ,'>N 
(5.3.18)  or  equivalently  (5.4.1)  to  (5.4.3)  can  be 
implemented  as  a  parallel  connection  1-0  filters,  one 
for  each  value  of  p^,  and  having  impulse  responses  as 
given  in  (5.4.3).  Note  that  the  input  y^-1*(.)  to  the 
Pj^-th  such  filter  is  obtained  from  y*i_1^(*)  by  shifting 
the  samples  by  an  amount  p^  and  then  by  considering  only 
every  M^-th  sample.  The  resulting  scheme  for 
implementation  is  shown  in  Figure  5.4. 


It  is  important  to  notice  that  for  interpolation  as  well 
as  for  decimation  the  computation  of  y  ( . )  in  (5.3.10) 
and  in  (5.4.1)  for  different  values  of  m^,...m 


i-1' 


Vi+1'***VN  are  completely  independent  of  one  another,  and 
thus  can  be  ■  carried  out  concurrently.  For  each  value  of 
the  set  of  (N-l)  tuples  just  mentioned  it  would  thus  be 
necessary  to  have  an  identical  copy  of  the  filter 
described  in  Figures  5.3  or  5.4  the  required  number  of 
such  copies  is  obviously  determined  by  the  size  of  the 

<i-l) 

support  of  the  signal  y  ( . ) . 

Furthermore,  in  spite  of  the  fact  that  as  shown  in  Figure 
5.1  and  5.2  the  entire  scheme  is  to  be  implemented  in  N 
decoupled  stages,  a  thorough  examination  of  the 
computational  scheme  reveals  that  it  is  at  least  in 
principle,  possible  to  begin  partial  computation  of 

(i+1)  (i+2) 

intermediate  signals  y  (.),...  etc.  even 
before  the  computation  of  y  (.)  is  completed,  thus 
potentially  resulting  in  further  speedup  in  arithmetic.  To 
exemplify  this  situation  in  the  2-D  case,  it  may  be 
noticed  that  since  the  copies  of  row  interpolators  all 
operate  in  parallel  -  each  on  one  row  of  the  input  image  - 
if  each  row  interpolator  is  made  to  sequentially  process 
the  rows  from  the  same  edge  of  the  picture  frame  then  as 
computati^)  proceeds,  the  columns  of  the  intermediate 
signal  y  (.)  begin  to  make  themselves  available  to  the 
column  interpolator ^^f  the  succeeding  stage  even  before 
the  computation  of  y  (.)  is  completed. 

In  what  follows  the  implementation  of  the  1-D  shift 
invariant  filters  g_  (n.)  for  the  interpolator  or  the 
pp  ( r ^ )  for  the  decimators  will  be  sought  in  direct  form 
noArecursive  (FIR)  structures.  Other  structures  can  be 
also  used,  but  our  choice  is  motivated  by  the  particular 
nature  of  the  interpolation  and  decimation  filters,  which 
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■akes  the  design  of  FIR  structure  especially  simple  as 
well  as  by  the  fact  that  such  structures  can  be 
conveniently  implemented  in  systolic  (pipeline) 
architecture  [2]  [10].  Figure  5.5  shows  an  example  of  such 
an  1-D  nonrecursive  structure,  which  occurs  in  each  branch 
of  the  polyphase  filters  for  the  case  when  the  filter 
order  is  equal  to  NQ.  Since  the  design  of  such  1-0  filters 
is  extensively  documented  [4]  in  the  literature  we  shall 
not  undertake  the  discussion  of  this  issue  here,  but  only 
highlight  the  fact  that  by  multidimensional  sampling  rate 
alterations  can  be  performed  at  a  higher  speed  by 
exploiting  in  parallel  several  copies  of  the  type  of  1-0 
filter  structures  discussed  in  existing  signal  processing 
literature.  Reduced  computation  time  may  thus  be  achieved 
at  the  expense  of  increased  amount  of  hardware 
requirement.  Furthermore,  it  may  be  noticed  that  for  the 
implementation  schemes  under  consideration  the  entire 
interpolator  or  decimator  have  a  high  degree  of 
modularity. 


5.5.  Design  examples  in  two-dimensions . 


It  follows  from  equations  (5.3.8b)  and  (5.4.3)  that  the 

impulse  responses  gm  (n^)  and  p^  (n^)  which  constitute  the 

respective  branches  the  1-D  polyphase  filter  structures 

for  the  interpolator  and  the  decimator  are  space  shifted 

versions  of  impulse  response  of  an  ideal  low-pass  filter. 

In  the  following  design  examples  the  window  method  of 

designing  FIR  filters  was  used  for  realization  of  the 

impulse  responses  just  mentioned.  The  Kaiser  window  was 

used  for  each  of  the  following  examples,  in  which  &-ripple 

in  the  passband  and  stopband  from  ideal  response; 

u>  -highest  frequency  of  interest  in  the  input  signal; 
c 

u  -passband  edge  frequency;  <•>  -stopband  edge  frequency; 
p  s  * 

and  N q —  required  filter  order.  The  details  of  design  can 
be  found  in  [14]. 


Interpolator:  L-2. 

The  following  choices  are  made  5-0.1,  w  -.65rt,  «  -w  /2L, 

c  pc 

»#-r/L.  Then  Ng-5.  The  resulting  polyphase  filter  was  used 
in  each  of  the  two  stages  of  implementation.  The 
performance  of  the  interpolator  was  tested  for  three 

different  test  signals:  (SI)  x(n. ,n, )-sin( r )/r,  where 

2  2  1  * 
r-/(u  +v  )  with  u,  v-  0,+l,  +2,... +25  as  shown  in  Figure 

5.6a;  (S2)  the  "Jet"  image  of  size  64  x  64  pixels,  with  4 

bits/pixel  resolution  as  shown  in  Figure  5.6b;  and  (S3) 

the  "Saturn”  image  of  size  64  x  64  pixels,  with  resolution 

4  bits/pixel  as  shown  in  Figure  5.6c.  The  interpolated 

signals  are  shown  in  Figures  5.7a,  5.7b  and  5.7c 

respectively. 
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L-3.  Again  8-0.1,  a>c«.65it,  o>p-<i>c/2L,  ws«n/0.8L,  NQ-5,  and 
Figures  5.8a,  5.8b,  and  5.8c  show  the  result  of  the 
interpolation  scheme  corresponding  to  input  signals  SI,  S2 
and  S3  respectively. 


Decimation:  M-2. 

The  following  parameters  were  chosen  for  each  branch  of 
the  polyphase  filters  of  the  basic  module  for  each  stage. 
8-0.7,  uc-.65x,  »p-o»c/1.3M,  Consequently  NQ-9.  The 

composite  filtering  scheme  was  tested  on  the  following  set 
of  test  signals.  (Tl)  xfn^njl-sint r)/r,  where  r-/(u^+v^ ) , 
u,  v-  0,  +1,  +2,... +75  as  shown  in  Figure  5.9a;  (T2)  the 
128x128  pixel,  4  bits/pixel  "Jet"  data  as  displayed  in 
Figure  5.9b  and  (T3)  the  128x128  pixel,  4  bits/pixel 
"Saturn"  data  as  displayed  in  Figure  5.9c.  The  respective 
decimated  signals  are  shown  in  Figures  5.10a,  5.10b,  and 
5. 10c. 

H-3.  Here  8-0.1,  «c-.65n,  Wp-«C/2L,  «s-x/.8L.  Consequently 
Nq-5.  The  input  signals  are:  (Ul)  same  as  in  (Tl)  above; 
(U2)  the  192x192  pixel,  4  bits/pixel  "Jet"  data  of  Figure 
5.8b  and  (U3)  the  192x192  pixel,  4  bits/pixel  "Saturn" 
data  of  Figure  5.8c.  The  corresponding  decimated  signals 
are  as  shown  in  Figures  5.10a,  5.10b  and  5.10c 
respectively. 

All  progams  were  written  in  a  DEC10  computer  in  sequential 
mode  as  opposed  to  parallel/  pipeline  modes  suggested  in 
the  present  study. 
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5.6.  Conclusion 


Motivated  by  practical  applications,  a  highspeed  computing 
scheme  for  the  interpolation  and  decimation  of  a  broad 
class  of  multidimensional  signals,  which  can  potentially 
derive  advantages  from  parallel,  pipelined  and  modular 
implementation  has  been  proposed.  Only  interpolation  or 
decimation  by  integer  factors  are  discussed,  but  sampling 
rate  alterations  by  non-integer  rational  factors  can  also 
be  performed  via  the  techniques  discussed  by  cascading  in 
two  successive  stages  interpolators  and  decimators  of  the 
type  discussed  here.  The  computational  scheme  is  based  on 
the  frequency  domain  description  of  the  signal,  and  uses 
one-dimensional  interpolation  or  decimation  filters  as  a 
basic  module  for  implementation.  Only  nonrecursive 
implementation  of  this  basic  module  has  been  considered  in 
the  present  study.  Apart  from  the  ease  of  design  and 
convenience  of  pipelineabili ty ,  the  choice  of  such 
structure  may  have  advantages  in  applications  such  as 
image  processing,  where  linear  phase  is  a  highly  desirable 
characteristic  of  the  processing  scheme.  The  signal  has 
been  assumed  to  be  sampled  according  to  the  rectangular 
cartesian  sampling  scheme.  Other  sampling  geometries  as 
multidimensional  (N>2)  generalizations  of  the  hexagonal 
sampling  scheme  [15]  prove  to  be  more  economical  in  terms 
of  the  number  of  samples  per  unit  volume  required  to 
represent  a  bandlimited  signal.  Interpolation  and 
decimation  schemes  for  multidimensional  signals  sampled  in 
such  geometries  and  their  implementation  in  currently 
emerging  highspeed  architectures  remain,  however,  to  be 
investigated . 
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CHAPTER  6 


CONCLUSIONS 

Sinple  and  rigorous  proofs  of  results  on  tests  for  the 
property  that  a  multivariable  polynomial  be  devoid  of  zeros 
in  the  closed  unit  polydisc  are  given.  The  proof  technique 
rests  on  a  complete  formulation  of  the  fact  that  the  zeros  of 
a  polynomial  are  continuous  functions  of  its  coefficients.  It 
is  shown  that  all  other  stability  related  results  can  be 
derived  in  this  manner. 

New  classes  of  multivariable  polynomials  arising  in  studies 
of  passive  multidimensional  systems  have  been  identified  and 
their  properties  have  been  studied.  In  particular, 
polynomials  occurring  as  the  numerators  and  denominators  of 
multivariable  reactance  functions  and  positive  functions,  are 
characterized.  Related  properties  of  these  and  other  classes 
of  multivariable  Hurwitz  polynomials  are  also  studied.  A 
nontrivial  test  for  the  property  of  positivity  of  rational 
functions,  holomorphic  in  a  domain,  in  terms  of  their 
behavior  on  the  distinguished  boundary  is  formulated. 

The  problem  of  structurally  passive  synthesis  of 
multidimensional  digital  filters  as  a  cascade  interconnection 
of  more  elementary  building  blocks  has  been  addressed  via  the 
factorization  of  the  associated  discrete  lossless  two-port 
transmission  matrix.  Necessary  and  sufficient  conditions  for 
the  factorization  to  be  feasible  are  obtained.  In  particular, 
it  is  shown  that  in  one-dimension  the  factorization  can 
always  be  performed,  and  as  a  consequence,  known  filter 
structures  fall  out  as  special  cases  of  the  results 
developed.  Thus,  an  alternative  algorithm  for  synthesizing 
one  dimensional  structurally  passive  digital  filters  is  also 
obtained. 
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The  problem  of  sampling  rate  alteration  of  deterministic 
multidimensional  signals  is  addressed  on  the  basis  of 
frequency  domain  description  of  the  signal.  It  is  shown  that 
the  problem  can  be  formulated  in  such  a  way  that  solutions 
can  be  obtained  via  filtering  techniques  known  for 
one-dimensional  signals.  Fast  non-recursive  implementations 
of  such  interpolation  and  decimation  schemes  are 
investigated.  The  resulting  algorithms  can  be  potentially  be 
implemented  in  a  combination  of  parallel  and  pipelined 
architecture.  Experiments  with  digitized  images  are  also 
reported  to  demonstrate  the  performance  of  the  designed 
interpolation  and  decimation  schemes. 

The  fundamental  results  developed  in  the  present  report  open 
up  ways  of  investigation  into  a  large  number  of  problems  of 
both  theoretical  and  practical  importance  in  the  area  of 
multidimensional  signal  processing.  Efficient  test  procedures 
for  the  various  classes  of  polynomials  identified  in  chapter 
3,  namely  the  scattering  Hurwitz,  reactance  Hurwitz  and  the 
immittance  Hurwitz  polynomials  etc.  are  lacking  and  needs  to 
be  developed.  Since  polynomials  of  this  type,  particularly 
the  scattering  Hurwitz  polynomials,  enter  into  the 
description  of  passive  systems  in  a  fundamental  manner,  this 
should  prove  to  be  an  important  step  in  designing  various 
types  of  multidimensional  filters.  In  this  context,  a 
detailed  study  into  the  properties  of  discrete  domain 
counterparts  of  the  various  multidimensional  polynomials 
discussed  in  chapter  3  and  their  testing  procedures  also 
remain  to  be  carried  out.  In  the  area  of  synthesis  of 
structurally  passive  multidimensional  digital  filter  design 
the  problem  has  been  addressed  only  in  the  context  of 
synthesis  in  cascade  type  structure  via  the  factorization  of 
the  transmission  matrix  associated  with  a  multidimensional 
lossless  two-port.  Other  possibilities  of  investigating 
synthesizability  of  lossless  multidimensional  two-ports  in 
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structures  other  than  the  cascade  structure  also  exist,  e.g., 
via  the  factorization  of  the  hybrid  matrix  or  the  transfer 
function  matrix  itself.  The  need  for  this  investigation  is 
strongly  felt  in  view  of  the  result  established  in  chapter  4 
that  in  a  generic  situation  multidimensional  lossless 
two-ports  may  not  be  synthesizable  as  cascade  interconnection 
of  most  elementary  passive  building  blocks.  Moreover,  exactly 
how  these  synthesis  schemes  can  be  utilized  in  special  cases 
of  practical  interest  when  the  frequency  response  of  the 
filter  is  required  to  have  certain  symmetries,  for  example, 
spherical  symmetry  (image  processing  applications)  or  planar 
or  conical  symmetry  (direction  finding  applications)  remains 
to  be  investigated.  Attention  has  only  been  restricted  to  the 
quarter  plane  type  recursive  schemes  so  far.  Other  recursive 
schemes  such  as  the  symmetric  or  the  asymmetric  half  plane  or 
the  fully  recursive  half  plane  recursive  schemes  and 
multidimensional  generalizations  thereof  also  needs  to  be 
considered. 

The  close  relationship  between  passive  filtering  and 
modelling  of  stationary  or  nearly  stationary  stochastic 
processes  is  well  known  for  one-dimensional  signals.  The 
results  of  the  present  investigation  can  thus  be  potentially 
utilized  towards  resolving  problems  in  the  domain  of 
modelling  of  random  fields  e.g.,  in  (spectral)  estimation, 
and  prediction  problems  associated  with  multidimensional 
signals.  Multidimensional  extensions  of  various  time/space 
varying  adaptive  filtering  schemes  can  also  prove  to  be  an 
important  topic  of  future  research  in  this  context. 

Since  generic  multidimensional  signal  processing  tasks  are 
severly  computation  intensive  the  feasibility  of  implemention 
of  the  algorithms  resulting  from  the  above  studies  in 
high-speed  architechtures  is  also  an  important  area  of 
investigation.  For  example,  in  the  specific  problem  of  fast 
image  interpolation  and  decimation  dealt  with  in  chapter  5 
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details  of  issuss  relating  to  implementation  in  systolic  VLSI 
and/or  optical  architechtures  remain  to  be  studied  and  can 
form  a  topic  of  future  research. 
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Figure  3.1 

Hierarchical  relationship  between  various  classes 
of  multidimensional  stable  polynomials. 
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Figure  4.4 

Elementary  section  of  type  1. 


e2 


Figure  4.5 

Elementary  section  of  type  2 
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Figure  5.8  Interpolated  (Si),  (S2) ,  (S3).  L=3 
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Figure  5.10  Decimated  (Ul) ,  (U2) ,  (U3) .  M-3. 


MISSION 

cf 

Rom  Air  Development  Center 

RAPC  plant  and  execute*  research,  development,  test 
ana  setected  acquisition  programs  In  support  of 
Command ,  Contact,  Communication*  and  Intelligence 
(C* 1}  activities.  Technical  and  engineering 
support  mlthln  areas  of  competence  Is  provided  to 
CSV  Program  Offices  [PCs }  and  other  ESV  elements 
to  perform  effective  acquisition  a*  CM  systems. 

The  areas  oj  technical  competence  include 
communications,  command  and  control,  battle 
management.  Information  processing,  surveillance 
■sensors.  Intelligence  data  collection  and  handling, 
solid  state  sciences ,  electromagnetics ,  and 
propagation,  and  electronic,  maintainability, 
and  compatibility. 


